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ABSTRACT 


In  the  past,  the  mobility  of  atmospheric  ions  and  the  num¬ 
ber  density  in  different  mobility  ranges  (ion  spectrumV have 
been  measured  in  the  terrestrial  atmosphere  in  the  heights 
from  gound  to  about  5  km.  Electric  conductivities  of  the 
atmosphere  have  been  measured  up  to  about  30  km  and  ion 
number  densities  to  about  75  km.  To  extend  the  possibility 
of  ion  spectrum  measurements  up  to  the  same  height,  a  new 
method  has  been  developed  and  tested  in  the  laboratory.  Its 
application  in  the  free  atmosphere  is  being  prepared.  The 
instrument  may  be  used  in  planetarian  atmospheres. 

The  method  consists  of  a  GERDIEN-type  "differental  ion 
counter  of  the  second  order,  "  which  provides  a  predeter¬ 
mined  location  of  ion  intake  and  a  number  of  separate  re¬ 
ceiving  electrodes  for  the  ion  impact.  The  driving  voltage 
is  ac;  its  amplitude  is  increasing  downstream  in  the  chamber. 

A  low-pressure  wind  tunnel  for  continuous  operation,  most¬ 
ly  in  the  subsonic  range,  has  been  developed  and  construc¬ 
ted  for  the  laboratory  experiments.  This  tunnel  and  its  pos¬ 
sible  applications  for  other  purposes  are  discussed. 

An  outline  of  the  theory  of  the  different  GERDIEN  type  ion 
counters  is  communicated.  The  laboratory  experiments 
are  described,  and  ion  spectrums  obtained  in  the  low-pres¬ 
sure  wind  tunnel  are  presented. 
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PARTI 

THEORY 


1.  1  ATMOSPHERIC  IONS 


1.1.1  GENERATION,'  LIFE  TIME,  AND  DESTRUCTION  OF  IONS 

The  existence  of  charged  particles  in  the  atmosphere  --  aerosol,  groups  of  mole¬ 
cules,  molecules,  atoms,  subatomar  particles  --  is  causedby  severalprocesses, 
the  relative  significance  of  which  varies  with  altitude  and  in  some  cases  with  time. 

Three  main  groups  of  ion  generation  in  a  volume  under  consideration  may  be  dis¬ 
tinguished:  (1)  the  entrance  of  charged  particles  from  outside  the  volume  or 
even  from  outside  the  atmosphere  (such  as  natural  or  anthropogeneous  aerosol 
particles,  blown  up  snow  or  sand,  smoke,  and  the  like,  or  penetrating  subatomar 
particles,  electrons,  protons,  alpha-rays,  mesons);  (2)  the  simultaneous  gener¬ 
ation  of  two  (or  more)  charged  particles  of  opposite  sign  by  some  superatomar 
"breaking-up  process"  (such  as  breaking  of  ice-crystals,  splashing  of  raindrops, 
and  so  forth,  within  the  atmosphere);  and  (3)  the  simultaneous  generation  of  a  free 
electron  and  a  positive  ion  by  a  "subatomar  breaking-up  process"  (the  usual 
"ionization"  by  removing  an  electron  from  an  atom  or  a  group  of  atoms). 

The  average  life  times  of  these  ions  differ  by  many  orders  of  magnitude.  Gen¬ 
eration,  variation,  and  destruction  of  atmospheric  ions  is  a  dynamic  process. 

The  number  of  ions  of  a  particular  kind  to  be  found  depends  on  the  time  constant 
for  that  particula:  ion  species,  and  this  time  constant  depends  on  atmospheric 
density  as  well  as  on  the  number  of  the  different  partners  with  which  the  ion  under 
consideration  may  react. 

Free  electrons  close  to  the  ground  are  attached  within  some  10 ”8  second  to  a  neu¬ 
tral  particle,  thus  forming  a  negative  ion.  In  the  tropopause,  the  time  constant 
(i.  e.  ,  the  time  in  which  l/e  of  all  freshly  generated  free  electrons  disappears) 
comes  to  about  10“^  second,  in  the  stratopause  it  reaches  about  10”2  second,  and 
in  the  mesopause  tenths  or  hundredths  of  seconds.  In  all  these  ranges,  attachment 
is  the  dominant  process  for  electron  removal,  but  the  combination  with  positive 
ions  gains  relative  significance  with  increasing  altitude.  Positive  and  negative 
"small  ions",,  i.  e.  ,  the  ions  which  carry  the  atmospheric  electric  current,  have 
life  time  constants  in  the  order  of  tenths  to  several  hundredths  of  seconds  close 
to  ground,  depending  on  the  contamination  of  the  atmosphere  by  aerosol  particles. 
In  greater  heights,  the  average  life  time  of  negative  ions  first  increases,  then 
flattens  off,  and  finally  goes  drastically  down,  reaching  values  of  a  few  seconds 
in  the  altitude  of  the  mesopause.  The  life  time  constant  for  pos  .ive  ions  gener¬ 
ally  increases  with  height.  Charged  aersol  particles  have,  in  general,  a  much 
longer  life  than  smaller  ions. 
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A  typical  pattern  may  be  presented  by  showing  the  scheme  of  WAIT  (1934)  for  the 
life  of  small  ions  in  the  troposphere  (figure  1). 


Figure  1  AN  EXAMPLE  FOR  THE  CREATION  OF  A  PAIR  OF  "SMALL  IONS"  OR“CLUSTER  IONS” 

IN  THE  LOWER  ATMOSPHERE  (AFTER  WAIT  1934): 

(A  neutral  molecule  (a)  is  met  by  radiation,  ond  one  positive  molecular  ion  (b)  and  one  fr*’e 
electron  (c)  are  formed.  The  free  electron  attaches  to  another  neutral  molecule  and  forms  a 
negative  molecular  ion  (d).  Attachment  of  more  neutral  molecules  brings  cluster  ions  (e)  into 
existence,  which  later  may  combine  or  attach  to  aerosol  porticles.) 

/  »« 

The  average  total  duration  tj  through  1 4  is  only  about  10"°  second  (ISRAEL,  (1957), 
p.  20).  When  one  moves  upwards  in  the  atmosphere,  the  time  durations  tj  through 
t4  are  enlarged,  and  chances  that  c  and  b  or  d  and  b  (figure  1)  meet  and  thus  can¬ 
cel  each  other  are  increased. 

The  life  of  a  particular  species  of  ions  is  terminated  by  several  processes,  which 
again  depend  in  their  relative  significance  on  height.  Three  main  possibilities  may 
be  distinguished  here:  two  ions  with  opposite  electric  charges  may  disappear  sim¬ 
ultaneously  with  their  electric  charges  by  combination;  an  ion  may  switch  over 
from  one  "group  of  ions"  to  another  such  as  electrons  to  molecular,  molecular  to 
small,  small  to  large  ions,  etc,,  by  attachment  or  detachment;  and  an  ion  may  re¬ 
main  within  the  same  "group",  but  still  change  its  nature.  Attachment  and  detach¬ 
ment  processes  may  involve  one  charged  and  one  uncharged  particle  (free  electron 
to  neutral  molecule,  or  liberation  of  an  electron  from  an  uncharged  aerosol  par¬ 
ticle,  and  other  cases)  or  two  charged  particles  of  the  same  sign  (in  this  case  at 
least  one  rather  coarse  aerosol  particle  which  is  able  to  hold  more  than  one  ele¬ 
mentary  charge  is  involved).  Change  of  the  nature  of  an  ion  may  occur  in  many 
different  ways:  a  charged  raindrop  may  partially  evaporate,  a  charged  ice  crystal 
may  melt,  a  clustered  ion  may  lose  or  attract  one  neutral  molecule,  a  charged 
atom  may  combine  with  a  neutral  atom,  charge  transfer  may  occur,  such  as 
0*  *  02  >  0  ♦  02  and  many  others.  The  fact  that  we  have  differentiated  between  the 
two  last  possibilities  is  slightly  arbitrary,  depending  on  the  denomination  "group 
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of  ions"  as  used  above.  It  is  justified  by  the  habit  of  quoting  differential  equations 
for  the  variations  in  number  densities  without  distinguishing  between  different  ions 
of  the  same  "group".  However,  for  our  later  mobility  discussions  such  differen- 
tiation  is  essential. 

1.1.2  INDICATION  TO  THE  DIFFERENTIAL  EQUATIONS 
FOR  ION  EQUILIBRIUM 

Analytical  expressions  for  the  variation  of  number  densities  of ‘ions,  dn/dt,  have 
been  proposed  in  great  numbers.  The  attempt  to  write  down  a  complete  set  of  dif¬ 
ferential  equations  covering  all  ion  generation,  ion  trano^er,  and  ion  decay  pro¬ 
cesses  for  all  heights  and  times  would  constitute  a  challenging  task  for  a  great  com 
puter  program  and  is  not  made  here.  The  simplest  equation  of  them  all,  already 
used  by  THOMSON  and  RUTHERFORD  (1896)  and  RUTHERFORD  (1897),  is  still 
the  most  important  one.  It  describes  the  ionic  state  in  an  environment  in  which 
only  ions  of  one  kind,  but  with  two  signs,  exist: 

dn/dt  =  q  -  a ; 
with 

n  the  number  density  of  ions  of  one  sign  (equal  to  that  of  ions  with  the 


q  the  number  of  ions  of  one  sign  generated  per  unit  volume  and  unit  time 
a  the  so-called  "recombination"  coefficient 
For  ion  equilibrium,  dn/dt  must  be  zero. 


v.  SCHWEIDLER  (1919,  19  24)  extended  the  discussion  by  including  large  ions 
with  one  elementary  charge  and  uncharged  aerosol  particles,  applying  an 
additional  "recombination  coefficient,  "  which  in  fact  is  a  composed  coefficient 
for  true  combination  between  ions  of  different  sign  and  for  attachment  of  small 
ions  to  uncharged  particles.  On  this  basis,  much  research  work  has  been 
done  to  obtain  quantitative  values  for  the  different  coefficients  and  the  discus¬ 
sions  have  been  extended  to  intermediate  ions  (e.  g.,  WAIT  1935),  but  only  ions  with 
one  elementary  charge  have  been  considered.  Using  five  different  coefficients  for  ion 
production,  four  different  combination  coefficients,  and  two  attachment  coefficients 


*  The  most  important  contributions:  1916:  KENNEDY;  1918,  1919,  1924:  v.  SCHWEIDLER;  1925:  NOLAN  et  al;  1927: 
HESS,  NOLAN  et  al;  1928:  HESS,  1929:  NOLAN  et  al.;  1933:  WHIPPLE,  WRIGHT;  1934:  HOGG;  1935:  SCRASE,  WAIT; 
1936:  WRIGHT;  1940:  FLEMING,  FISH  &  SHERMAN;  1941:  v.  SCHWEIDLER,  THELLIER;  1944:  NOLAN  et  al.;  1946: 
WAIT;  1949!  BRICARD;  and  the  different  works  of  ISRAEL,  reported  in  1957;  some  more  literature  is  given  by  GISH 
and  SHERMAN  1940  and  PARKINSON  1943. 


ISRAEL  (1957'p.  .189)  gives  a  set  of  five  differential  equations,  in  which  small  ions 
and  large  ions  and  neutral  aerosol  particles  are  combined,  but  none  has  more  than 
one  elementary  charge.  Detachment  coefficients  are  included  in  the  ion  production 
coefficients.  No  free  electrons  are  discussed.  For  the  inclusion  of  multiply-charged 
aerosol  particles,  PEDERSEN  (1964)  applies  the  plasma  probe  theory  by  LANGMUIR 
et  al  (1923  to  1939,  in  particular:  MOTT -SMITH  and  LANGUIR,  1926,  and  LANGUIR 
and  COMPTON,  1931)  to  a  non-charged  or  charged  aerosol  particle,  and  he  derives 
a  set  of  equations  in  which  the  number  of  charges  on  an  aerosol  particle  is  included. 
Since  he  is  concerned  about  the  height  range  between  40  and  80  km,  his  equations 
also  reflect  the  inclusion  of  attachment  coefficients,  detachment  coefficients,  and 
combination  coefficients,  which  play  a  role  in  these  heights.  These  have  not  been 
considered  by  the  earlier  workers  mentioned  in  the  footnote  above,  but  have  been 
added  in  the  last  decades.**  In  his  equations  PEDERSEN  assumes  that  the  aerosol 
particles  have  either  no  charge  or  a  number  of  Z  elementary  charges;  thus,  he  sim¬ 
plifies  somewhat  the  given  distribution  of  charges  on  the  particles. 

This  fact,  together  with  the  already  indicated  situation  that  in  all  these  proposed 
aeronomic  reaction  equations  the  possible  differences  between  ions  of  one  "kind" 
are  not  considered,  has  the  consequence  that  most  of  the  applied  coefficients  are 
related  to  a  certain  mixture  of  ions,  and  not  to  one  species  alone. 

11.3  BRIEF  DISCUSSION  OF  THE  CHEMICAL  NATURE  OF  IONS 

On  the  other  hand,  chemical  reaction  equations,  which  may  play  a  role  in  the  atmos¬ 
phere,  havv.  been  listed  in  great  numbers.  The  book  by  NAWROCKI  and  PAPA 
(1961,  chapter  3)  provides  a  selection.  Number  densities  are  not  much  considered 
here,  although  indications  on  the  energy  transfers  and  on  probability  are  given. 
According  to  COLE  and  PIERCE  (1965),  the  most  important  positive  ions  are  prob¬ 
ably  NO+  and  02+  but  N2+  and  complex  ions,  such  as  03+  and  N4+,  may  be  present 
in  significant  quantities  while  negative  ions  are  most  probably  02~  and,  at  greater 
heights,  0 In  lower  altitudes,  in  particular,  other  ions  have  been  considered 
as  well,  e.  g.  ,  H20+  .  There  may  sti*ll  be  found  other  ions  in  small  numbers,  such 
as  H“  ,  OH  “ ,  N02~,  Oj-  ,  and  others,  and  large  ions  in  addition  to  these. 
PEDERSEN  (1964)  discusses  the  existence  of  ions  in  the  mesosphere  which  move 
but  slowly  under  the  influence  of  the  electric  fieled  (heavy  ions).  NARCISI  and 
BAILEY  (1965),  measuring  with  quadrupole  massenfilter,  moving  with  a  velocity 
greater  than  that  ol  sound,  and  after  pumping  down  to  reduce  collisions  within  their 
instrument,  observe,  at  altitudes  between  64  and  82.  5  km,  ions  of  the  species 
OH*  ,  N2  ,  NO+  ,  02  ,  but  also  ions  of  mass  19  and  37  interpreted  to  be  HjO*- ,  and 
possibly  H5O2  .  That  ions  of  the  last  mentioned  species  may  exist  has  been 

shown  in  the  laboratory  in  the  range  from  40  to  80  N/m^  (0.  3 . 6  torr)  by 

KNEWSTUBB  and  TICKNER  (1963). 


**1958:  SMITH,  IHJRCH  and  BRANSCOMB;  1959:  CHAN1N,  PHELPS,  and  BIOND1;  1961:  GUNTON  et  al.,  HOLT  ct  a)., 
Phelps  and  PACK;  1963:  B10NDI  (quoted  after  PEDERSEN). 
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NARCISI  and  BAILEY  (1965)  also  observed,  using  another  measuring  mode,  many 
positive  ions  of  atomic  mass  greater  than  45  in  altitudes  between  65  and  75  km,  while 
the  defined  metallic  ions  such  as  Fe+  ,  NV,  Ca+, and  have  but  small  number  densities 

below  the  mesopause.  For  the  transition  from  the  measured  current  to  number  den¬ 
sities,  again  only  one  charge  per  particle  is  assumed. 

For  all  measurements  made  from  rockets,  the  possibility  must  be  investigated  that 
species  of  ions  found  in  the  measuring  results  have  been  injected  into  the  atmosphere 
by  the  rocket  itself,  and  thus  do  not  constitute  genuine  atmospheric  ions.* 

1.1.4  CONCLUSION 

Considering  the  discussion  of  this  chapter  on  atmospheric  ions,  one  comes  to  the 
conclusion  that  a  "measuring  stick"  for  ions,  which  attaches  one  particular  number 
to  each  ion  species,  is  needed.  The  natural  measuring  stick  for  this  purpose  is  mo¬ 
bility,  as  defined  and  discussed  in  the  next  chapter.  For  most  cases,  mobility  pro¬ 
vides  a  unanimous  "label,  "but  some  ambiguity  may  appear  in  particular  cases.  This 
is  rather  probable  in  the  case  of  large  ions  carrying  more  than  one  elementary  charge. 
However,  even  in  these  cases,  mobility  is  the  most  promising  tool  for  the  separation 
of  different  groups  of  ions  from  each  other  and  for  the  separation  of  different  ions 
within  the  same  "group.  " 


Supplement  made  when  proofreading:  r.  recent  paper  by  LINDZEN  and  GOODY  (J .  Atm.  Sci.  22  (1965)  341)  provides  an 
impressive  evidence  for  the  complexity  of  the  region  between  30  and  70  km. 
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1.  2  IONIC  MOBILITY 


1.  2.  1  DEFINITION  OF  THE  PARAMETER  MOBILITY 

A  generally  applied  formula  for  the  relative  motion  of  the  two  types  of  particles 
("molecules")  in  a  binary  gas  mixture  follows: 


Vel 


nj  +  ri2 
nl  n2 
m  n  —  m 


Dj2  grad  n. 


D 


12 


n  i  m  j  +  n  2  ^2 

kBT 

nl  +  n2 

D12 

nj  mj  +  n2  m2 

kBT 

(nj  +  n2)2 

grad  p 


(1) 


mj  m2  (a j  -  a2) 


n  n  grad  T 

0<)  'j*  ° 

[Sydney  CHAPMAN  and  COWLING,  I960,  eq.  (14.  1/1),  with  substitutions  accord¬ 
ing  to  eqs.  (1.21/1),  (2.5/2),  (2.5/3),  (2.5/4),  (2.5/6),  (2.5/12),  and  (8.  3/5), 
and  written  in  our  own  symbols],  where: 

is  the  average  relative  velocity  between  the  two  constituents  of  the 
gas  [m/sec] 

n j  and  n0  are  the  number  densities  of  the  two  constituents  [m“^] 

Djj  is  the  coefficient  of  mutual  diffusion  [m2/sec] 


ml 


and  m2  are  the  masses  of  the  molecules  of  the  two  constituents  [kg] 

-23, 


k  is  the  BOLTZMANN  constant,  1.  38  x  10"  [Newton-meter  per  degree 


or  Watt-seconds  per  degree 


T  is  the  temperature  [degrees  KELVIN  j 
p  is  the  hydrostatic  pressure  [Newton/m^  ] 


a]  and  a2  are  the  accelerations  of  the  molecules  of  the  two  constituents, 
effected  by  external  forces  [m/sec2  1 

Dt  is  the  coefficient  of  thermal  diffusion  [nr /sec]. 
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For  the  significance  of  the  second  and  the  fourth  term  of  eq.  (1),  refer  to 
CHAPMAN  and  COWLING  (I960),  p.  143,  lines  10  to  12  of  the  second,  and  all 
of  the  third,  paragraph.  If  we  consider  local  conditions  in  the  atmosphere,  con¬ 
fined  to  dimensions,  in  which  grad  p  and  grad  T  may  be  neglected,  the  second 
and  the  fourth  term  in  eq.  (1)  vanish. 


The  first  term  represents  the  "gaskinetic"  diffusion,  which  occurs  if  the  two 
constituents  are  not  evenly  distributed  throughout  the  gas  (grad  0),  and  the 
third  term  represents  the  motions  of  one  of  the  constituents  through  the  other 
one,  if  external  forces  act  on  this  one  constituent,  and  not,  or  not  in  the  same 
amount,  on  the  other  one.  We  denote  with  the  index  1  the  ions,  and  with  index 
2  neutral  molecules,  and  we  restrict  our  discussion  to  only  one  kind  of  external 
force,  an  electric  field  (i.  e.  ,&2  =  0).  Then  we  may  write  for  the  motion  of  ions 
(of  one  kind)  through  a  gas,  constituted  by  only  one  type  of  neutral  molecules, 
the  following  equation: 


nj  +  f)2  nj  +  n  2  ^12 

-  Di-j  grad  n,  +  -  - 

rij  n2  n jm j  +  kgT 


nil  3 1  n^2 


(2) 


Replacing  the  external  force  acting  on  the  ion  with  the  electric  charge  e. 


mjaj  =  e 


(3) 


we  get: 


vrel  = 


‘‘I  +  n2  ni  +  n2 

-  Dt2  grad  n,  +  - 

nln2  ml 


eD^ 

kBT 


(4) 


where  E  is  the  electric  field  strength.  Thus,  the  average  relative  velocity  be¬ 
tween  the  two  constituents  of  the  binary  gas  in  this  case  (grad  T=grad  p  =  0) 
is  composed  of  two  components: 


vrel  =  <vrel>D  +  (vrel>E 


with 


„  nl  +  n2 

( vrel =  ”  TT —  D12  «rad  nl 


D 


n,  n 


1  2 


and 


<VI>E 


n  i  +  n2 


mi 

nj  — —  +  02 

m2 


e  P12 
k,,T 


E  . 


(5) 

(6) 


(7) 
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We  restrict  out  discussions  in  this  section  to  eq.  (7).  The  term  with  the  dif¬ 
fusion  coefficient  may  be  substituted  by  the  parameter  "mobility"  as  follows: 


e  d12  (8) 

k  a  -  . 


Expression  (8)  is  often  referred  to  as  the  EINSTEIN  relation.  In  MKSA-units, 
the  mobility  is  given  in  m2  V”1  sec-1  (for  the  derivation,  remember  l  Nm  =  1  Wsec). 
With  this  parameter,  eq.  (7)  comes  to: 

- ,  nl  f  n2 

(vrel)£  - - k  E  ;  (9) 

mj 

"1  —  +  n2 
m  2 

and,  if 

ml 

nj  <<  n2  and  nj  -  <<  n2  ,  (10) 

m2 


(vrel>E 


=  k  E  . 


(ID 


This  expression  leads  to  a  definition  of  the  parameter  "mobility":  In  a  slightly 
ionized  gas  (condition  (10)),  the  mobility  is  the  average  velocity  (drift  velocity) 
of  an  ion  under  the  influence  of  an  electric  field  of  unity  strength. 


If  we  measure  the  average  velocity,  or  the  field  strength  with  an  accuracy  of  1 
percent  the  term  "slightly  ionized"  means  that 

nj  +  nj 

- 1 

m, 

nl—  +  n2 

n^2 

In  the  free  atmosphere,  the  lightest  molecular  or  atomic  ions  may  be  hydrogen, 
and  the  heaviest  radon;  thus  mj/m2  would  vary  from  about  3/30  to  about  8.  If  we 
include  cluster  ions  (up  to  20  neutral  molecules  fixed  to  one  ionized  molecule), 
the  upper  limit  may  shift  to  20  or  more,  and  if  we  include  large  and  ultralarge 
ions,  it  may  go  up  to  about  1000.  This  means,  that  eq.  (12)  holds  for  molecular 
icns,  cluster  ions,  and  ultralarge  ions  up  to: 


<  0.01 


(12) 


Ion 

n;  less  than 

Molecular  ion 

1.  4  x  10' 3 

-4 

Cluster  ion 

5.  2  x  10 

Ultralarge  ion 

10“  5 

of  n 


21- 
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If  we  knov.  fhat  in  the  gas  under  consideration  only  ions  with  m j  =  m2  exist,  we 
may  permit  nigher  ionization  ratios.  However,  we  must  be  cautious  so  as  not 
to  hurt  other  limitations,  e.  g. ,  the  one  mentioned  in  the  next  paragraph  con¬ 
cerning  the  mutual  influence  between  ions.  For  electrons,  a  ratio  of  even  less 
than  10“J  applies.  Furthermore,  definition  (11)  is  applicable  only,  where  the 
CHAPMAN-COWLING  equation  (1)  is  valid.  That  means,  at  first,  that  we  must 
consider  the  average  behavior  of  the  ion  in  a  space  large  enough  to  allow  the  de¬ 
rivation  of  mean  values  from  a  great  number  of  particles  participating  in  a  great 
number  of  collisions,  and  do  not  study  the  fact  that  between  collisions  the  ions 
have  an  accelerated  motion.  Secondly,  the  number  of  ions  must  be  "so  small 
that  their  mutual  influence  on  each  other's  mean  velocities  is  negligible" 
(CHAPMAN  and  COWLING  (I960),  p.  321,  valid  not  only  for  more  than  one  species 
of  ions). 


Although  eq.  (1)  is  derived  only  for  a  binary  gas,  it  is  an  established crusto m 
to  speak  of  a  mobility  also  in  the  case  of  gases  composed  of  many  different  mol¬ 
ecules,  e.  g. ,  in  the  atmosphere.  In  this  case,  n2  represents  the  number  density 
of  all  neutral  particles,  and  m2  is  the  weighted  average  mass  of  them,  in  the  ter¬ 


restrial  atmosphere  up  to  more  than  hundreds  of  kilometers,  approximately 
29  x  1.  7  x  10"^7  kg. 


1.  2.  2  GENERAL  MACROSCOPIC  DISCUSSION  ON  MOBILITY 


In  this  section,  we  treat  the  gas  in  which  the  ion  moves,  as  a  continuum,  while 
in  the  following  section  (microscopic  discussion)  the  molecular  structure  of  the 
gas  is  considered.  Be  definition,  mobility  is  a  "macroscopic"  parameter. 

If  in  a  gas  containing  ions  there  is  a  difference  of  electric  potential  across  some 
distance,  there  will  be  an  electric  "conduction"  current;  and  if  there  is  an  electric 
current,  there  will  be  an  electric  field.  In  the  gas  mixture  under  consideration 
(slightly  ionized,  short  mean  free  paths,  many  collisions),  the  parameter  con¬ 
necting  electric  fields  and  "conduction^currents,  "  conductivity,  is  composed  of 
ion  number  density,  electric  charge  of  the  ions,  and  ionic  mobility.  Be  defini¬ 
tion  (eq.  (11)),  mobility  is  the  ionic  velocity  in  a  unity  field.  Accordingly,  the 
most  direct  method  to  measure  mobility  is  to  measure  the  velocity  of  the  ion, 
i.  e. ,  to  measure  a  current.  In  fact,  this  is  the  usual  method. 

There  are  two  principal  difficulties  related  to  this  measurement. 

At  first,  with  the  usual  ion  (or  electron)  velocities,  the  current  represented  by 
the  motion  of  a  single  elementary  charge  is  too  small  to  be  measured  with  pre¬ 
sent  methods.  Thus,  we  measure  the  sum  <->f  the  motions  of  a  great  number  of 
ions  which  may  be  different  kinds,  having  different  mobilities  (and  different  signs). 
Such  a  measurement  would  yield  an  "average  mobility.  "  For  a  more  informa¬ 
tive  measurement,  the  ions  must  be  separated  into  different  ranges  of  mobility 
prior  to  the  measurement  of  the  current.  If  this  is  done  successfully,  we  obtain 
a  mobility  spectrum,  i,  e.,  an  information  on  the  number  of  ions  in  each  range 
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of  mobility,  or 


dn 

dk 


f(k). 


(14) 


The  separation  must  be  done  prior  to  the  measurement,  because  in  the  usual 
methods  the  measurement  destroys  the  ion.  Also,  in  the  case  of  the  measure¬ 
ment  of  the  ion  spectrum,  the  current  measured  for  each  mobility  range  must 
be  represented  by  a  great  number  of  ions  moving,  and  moving  through  a  great 
number  of  collisions. 

This  is  not  only  required  by  the  fact  that  we  do  not  have  methods  to  measure  so 
small  a  current  as  the  one  representing  the  motion  of  just  one  or  a  few  ions. 

More  essentially,  the  integration  over  many  ions  and  many  collisions  is  required 
by  the  facts  that  ions  of  any  mobility  are  constantly  generated  and  destroyed,  and 
that  the  motion  of  an  ion  between  two  collisions  is  an  accelerated  one.  Thus, 
the  measuring  instrument  must  comprehend  a  sufficiently  large  volume,  and 
the  measurement  must  be  extended  over  a  sufficiently  long  period  of  time.  Both 
volume  and  time,  depend  on  the  gas  density.  They  must  be  expanded,  when 
measurements  in  the  upper  layers  of  the  atmosphere  are  intended,  and  the  at¬ 
tempt  becomes  unrealistic  if  atmospheric  pressure  becomes  smaller  than  about 
1  N/m“  (or  7  x  10”^  torr). 

Second,  any  measurement  of  the  current  provides  information  on  the  total  cur¬ 
rent  only,  and  does  not  distinguish  between  the  "conduction  current"  mentioned 
above  and  other  possible  types  of  current,  represented  by  the  motion  of  the  same 
group  of  ions  which  represent  the  conduction  current.  To  evaluate  these  possi¬ 
bilities,  we  should  discuss  the  different  kinds  or  shares  of  current  which  may 
exist  in  the  gas  or  atmosphere  under  consideration,  outside  and  inside  the  meas¬ 
uring  instrument.  Returning  to  eq.  (1),  we  want  to  state  the  introduction  of  elec¬ 
tric  forces  as  the  only  ones  acting  from  "outside"  eq.  (3)  is  a  restriction.  There 
may  be,  in  place  of  them,  or  in  additipn  to  them,  other  external  forces  such  as 
gravity  forces  or  the  inertial  forces  of  particles  which  have  been  accelerated  ex¬ 
tremely  strong  at  an  earlier  time  and  outside  of  the  volume  under  consideration 
(penetrating  particles).  Furthermore,  there  may  be  a  motion  of  the  binary  gas 
as  a  whole,  depending  on  the  selected  coordinate  system,  and  this  will  be  a  mo¬ 
tion  of  the  particles  indicated  with  the  index  1  as  well.  If  these  particles  (in 
which  we  include  the  penetrating  particles  also)  carry  an  electric  charge  all  these 
motions  will  represent  electric  currents.  If  some  or  all  of  the  mentioned  ef¬ 
fects  (seven  in  total)  are  in  about  the  same  order  of  magnitude,  and  act  on  the 
same  ion,  a  complicated  pattern  of  current  lines  may  be  the  result.  Usually, 
more  than  one  species  of  ions  exist,  and  the  pattern  is  even  more  complicated. 

Disregarding  again  grad  p  and  grad  T,  we  may  distinguish  between  the  following 
shares  of  the  electric  current: 
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"Conduction  current,  "  represented  by  eq.  (7),  applied  to  all  ions 
existing  in  the  gas  mixture 

"Convection  current"  (in  the  narrower  sense  of  this  denomination),  repre¬ 
sented  by  the  motion  oi  the  gas  mixture  as  a  whole,  including  the  ions 
of  different  kinds 

"Gravitational  current,  "  included  in  the  third  term  of  eq.  (1)  if  mj  ^  m-,,  and 
applicable  to  the  different  kinds  of  ions 

"Penetrating  current,  "  if  such  ions  exist 

"Diffusion  generated  electric  current,  "  represented  by  eq.  (6). 

Sometimes,  the  last  four  shares  are  combined  in  the  denomination  "convection 
current,  "  used  in  a  wider  sense  and  including  all  motions  of  ions  accelerated 
by  non-electric  forces.* 

Consequently,  to  derive  mobility  from  the  measured  value  of  the  current,  we 
must  ascertain,  that  only  the  "conduction  current"  is  considered.  This  r.  _y  be 
proven  by  additional  experiments  or  by  theoretical  computations  if  the  related 
parameters  are  known  sufficiently.  Fortunately,  this  task  is  often  an  easy  one. 
Convection  currents  in  the  narrower  sense  may  be  detected  easily,  if  not  in 
magnitude  so  in  direction,  and  if  the  conduction  current  in  an  instrument  can  be 
made  to  be  directed  perpendicular  to  these  convection  currents,  the  latter  ones 
are  easily  cancelled  out.  Gravitational  currents  and  penetrating  currents  do 
occur,  but  they  occur  only  occasionally,  and  often  are  too  small  to  be  considered, 
because  either  the  charges  transported  are  small  or  because  the  velocity  of  this 
transport  is  small.  In  the  troposphere,  strong  gravitational  currents  are  re¬ 
presented  by  the  falling  of  charged  rain  drops,  etc. ,  while  in  the  mesosphere 
and  stratosphere,  penetrating  electrons  and  protons  or  alpha  rays  may  constitute 
a  current.  While  charged  rain  or  snovj  is  easily  detected,  the  detection  of  pene¬ 
trating  charged  particles  may  be  difficult,  and  this  possibility  should  be  kept  in 
mind. 

Thus,  in  addition  to  the  conduction  current,  there  remains  mostly  the  "diffusion¬ 
generated  electric  current,  "  which  must  be  considered  for  all  these  measure¬ 
ments  as  a  possible  and,  in  principle,  always  present  disturbance. 

In  a  discussion  on  mobility  carried  out  by  applying  macroscopic  view  points,  a 
word  must  be  said  on  the  behavior  of  this  parameter  when  applied  to  different 
ions  and  in  different  environments. 


The  displacement  current,  although  in  atmospheric  processes  often  stronger  than  the  ether  shares,  is  not  discussed  here 
since  it  does  not  represent  a  motion  of  a  charge. 
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At  first,  it  must  be  stated,  that  the  mobility  of  the  same  kind  of  ion  in  the  same 
kind  of  gas  is  directly  dependent  on  the  inverse  value  of  the  density  of  the  gas:* 


kP’  T  JJp  _T_  (16) 

\T0  =  p‘  T0 

Second,  the  nature  of  the  atmospheric  ions  is  not  well  known.  It  is  still  a  con- 
t-'wersial  item,  whether  the  atmospheric  ions  in  the  troposphere  exist  of  ion 
clusters  or  of  molecular  ions.  Also,  the  chemical  nature  of  these  ions  is  a  mat¬ 
ter  of  discussion.  Furthermore,  it  is  expected,  that  the  nature  of  ions  changes 
with  height.  If  there  are  cluster  ions  close  to  ground,  more  and  more  non-clus- 
tered  ions  may  appear  if  height  increases,  and  finally  even  more  and  more  atomar 
ions  will  replace  the  molecular  ions.  The  transition  region  will  be  different  for 
chemically  different  ions. 

Since  generation  of  ions,  ion  equilibrium  maintainance,  and  decay  of  ions  is  a 
never  ceasing  dynamic  process,  the  time  constants  for  the  different  processes 
play  the  main  role  here,  as  has  been  discussed  in  section  1.  1,  above. 

Third,  mobility  is  defined,  according  to  eq.  (11),  as  being  independent  of  the 
field  strength.  This  condition  is  mostly  fulfilled  in  the  atmosphere.  However, 
under  certain  circumstances,  it  may  be  hurt  in  an  instrument.  If  the  electric 
field  increases,  or  the  pressure  decreases,  or  both,  such  that 


E 

—  >  200 
P 


Vm 

=  1.5  - 

N 


the  definition  of  eq.  (11)  is  no  more  applicable.  The  velocity  of  the  ion  does  not 
more  increase  linearly  with  the  field  strength;  instead  it  becomes  more  close  to 
a  relationship  such  as: 

v  =  f  (\/T)  .  H8) 

This  fact  will  be  discussed  later  on  (section  1.  7.  3). 

Fourth,  it  may  be  possible  that  ions  exist  which  can  be  destroyed  by  the  appli¬ 
cation  of  a  small  amount  of  energy.  We  know  from  laboratory  experiments  and 
theoretical  considerations  (e.  g, ,  E.  A.  MASON  and  SCHAMP,  1958;  DALGARNO, 
McDOWELL  and  WILLIAMS,  1958)  that  ions  with  bindings  energies  as  low  as 


♦  #  5/2 

This  'imple  relationship  is  often  used  to  "reduce"  measured  mobility  values  to  the  pressure  of  "ground  level",  or  10  N/m  . 

‘'uch  a  reduction  facilitates  the  comparison  of  different  kinds  of  ions,  which  may  exist  in  different  heights,  but  care  must 
be  taken  not  to  confi  sc  the  absolute  values. 
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0.  014  eV  may  occur  (Cs+  -  He),  and  we  expect  to  find  in  the  higher  atmosphere 
ions  with  binding  energies  in  the  order  of  1  eV.  Since  any  GERDIEN  measure¬ 
ment  applies  stronger  fields  to  the  atmosphere  than  normally  existing  in  undis¬ 
turbed  conditions,  the  breaking- up  of  such  ions  is  a  possible  disturbance*  and 
must  be  considered.  We  shall  submit  a  quantitative  discussion  on  this  problem 
for  the  application  in  the  mesosphere  later  on  (section  1.  7.  4). 

1.  2.  3  GENERAL  MICROSCOPIC  DISCUSSION  ON  MOBILITY 


Mobility  is  sometimes  referred  to  as  being  a  macroscopic  expression  for  the 
integral--taken  over  some  time  and  over  some  space--of  a  variety  of  frictional 
forces  which  turn  the  accelerated  differential  motion  of  the  ion  in  an  unaccel¬ 
erated  integral  motion  with  a  definite  velocity.  This  image  is  a  reminder  of  the 
motion  of  a  solid  particle  through  a  fluid  with  high  viscosity,  under  the  influence 
of  an  external  force,  such  as  gravity.  In  fact,  as  has  been  shown  by  ISRAEL 
(1957,.  p.  26  and  288)  the  hydrodynamic  law  of  STOKE S,.  as  modified  by 
CUNNINGHAM  and  MILLIKAN,  and  valid  downwards  to  particle  sizes  in  the  order 
of  the  mean  distance  between  the  particles: 


67TT/T j  Vj 

1  +  l/r1  [0.864+  0.290  exp  (-1.25  rj/L)] 


(19) 


smoothly  changes  over  into  the  formula  by  LENARD  (1913,  1919,  1920)  which 
is  supposed  to  be  valid  for  small  particles  upwards  to  the  order  of  the  mean  free 
path: 

p  _  77  P  (fi  +  f2)2  ^2  vj  (20) 

Q(n»j,  m2) 

where 


F  is  the  force  acting  on  the  particles  with  the  index  1 
Tj  is  the  dynamic  viscosity  of  the  gas  (atmosphere) 

rj  and  tj  are  the  radii  of  the  particles  with  indices  1  and  2  (ions  and  gas 
molecules) 

vj  is  the  stationary  terminal  velocity  of  the  particles  with  index  1 
l  is  the  mean  free  path 
P  is  the  density  of  the  gas  (atmosphere) 


It  may  be  noted  here  that  the  usually  observed  non-fulfillment  of  Ohm's  Law  in  the  measured  quantities  of  the  atmos¬ 
pheric  electric  field,  current  density,  and  conductivity  may  find  another  explanation  (additional  to  DOLEZALEK,  1960) 
in  the  possibility  of  changing  ionic  mobility  when  measuring  conductivity. 
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Q  is  a  factor  representing  the  masses  of  the  particles  with  the  indices  1  and 
2  (ions  and  gas  molecules) 

is  the  mean  velocity  of  the  gas  molecules. 

If  we  determine  a  "mobility"  also  for  the  case  of  the  STOKES-CUNNINGHAM- 
MILLIKAN  Law  as  well  as  for  the  LENARD  Law,  and  substitute  for  F  the  force 
acting  on  an  electric  charge  in  an  electric  field  of  unity  strength,  we  come  to 

1  +  1  hx  [0.864  +  0.29  exp  (-1.25  q/f)  ]  (21) 

k  =  -  e 

6  rr  rj  rj 

and 


3  1  +  m* 

8(1—  m*)  2m* 

e,  (22) 

2  \/  m*  3  +  m* 

*t* 

n  (1  +  m*)  1  +  m* 

- - - 

it  p  (rj  +  r2)“  v2 

where  we  have  now  introduced  the  explicit  value  for  the  Q  of  eq.  (20),  with 

m*  =  m  j/  (mj  +  TTI2)  (23) 


and 


mj  and  m2  being  the  masses  of  an  ion  and  of  a  gas  molecule, 

ej  is  the  electric  charge,  for  "small"  or  molecular  or  atomic 
ions,  being  1.6x  10"  19  coulombs. 

Since  at  present  we  want  to  discuss  the  fundamental  relations,  we  may  apply 
(21)  for  small  l  Aj  values,  so  that 

,  1  '  (24) 

k  =  -  e, 

6mit1 

and,  with  the  usual  formula  for  the  dynamic  viscosity,  written  with  our  sym¬ 
bols: 


1  m2  v2  (25) 
3  4n\f?.(  rj  +  r2)2 

(see,  for  this,  e.  g.  ,  SEARS  (1959),  eqs.  (1 3- 1 3)  with  p.  256  last  line) 


we  arrive  at  the  following  equation  for  the  "mobility"  in  STOKES'  Law: 


k 


y/T 


‘2  v2 


:i ; 


(26) 
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Introducing: 


where  we  have  taken  into  account  that  rj  >  >  ^  . 


we  may  also  write  this  as: 


On  the  other  hand,  if  we  make  in  LENARD's  Law  (22)  rrq  >  >  m->  and  rj  >>  ^  > 
and  apply  eq.  (27),  we  have: 


Not  discussing  the  fact  that  density  p  is  not  included  in  eq.  (28),  since  dynamic 
viscosity  is  independent  of  gas  density,  and  that  radii  and  masses  of  "ion"  and 
gas  molecules  have  different  functional  relationships  to  mobility  in  eqs.  (28)  and 
(29),  we  see,  that  in  both  cases  mobility  is  a  function  of  the  electric  charge  of 
the  ion,  the  size  and/or  mass  or  ion  and  gas  molecules,  and  the  temperature  of 
the  gas: 

k  =  f  (r,  m,  e,  T,  p)  ;  (30) 

where  p  cancels  for  very  large  charged  particles  (e,  g,  ,  electrified  raindrops). 

The  problem  of  giving  an  anaTytical  expression  for  eq.  (30)  is  extremely  com¬ 
plicated  and  far  from  being  solved. 

• 

The  attempts  begin  with  an  extended  discussion  between  LENARD  (1913,  1919, 
1920),  who  starts  from  the  concept  of  a  solid  impact,  and  LANGEVIN  (1902, 
1905,  and  later),  who  derives  one  set  of  his  equations  from  the  center-of-force 
theory,  which  is  expressively  doubted  by  LENARD.  In  the  center-of-force 
theory,  the  electric  charge  of  the  ion  and  its  mirror  charge  on  the  neutral  gas 
molecule  are  supposed  to  alter  significantly  the  free  path  of  the  ion,  and  here 
not  only  the  static  fields  play  a  role,  but  the  dipole  and  quadrupole  and  higher 
fields  also  may  be  discussed.  These  discussions  considered  gas  pressures  of 
about  1  atmosphere,  and  they  will  be  altered  if  we  consider  gas  densities  as 
they  occur  in  the  stratosphere  and  mesophere.  Eq.  (22)  is  the  basic  LENARD 
equation. 

A  basic  LANGEVIN  equation,  which  does  not  yet  consider  the  electric  forces, 
written  with  our  symbols  and  transferred  to  MKSA-system,  follows: 

C  / 

k  =  5.5  X  10>  — —=z —  ej  . 

mj  V2  v  m* 
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(31) 


Based  ontheorectical  considerations  and  laboratory  experiments,  a  series  of 
improvements  of  this  equation  has  been  proposed.  By  introducing  the  electri¬ 
cal  forces  (center-of-force-theory),  the  mobility  becomes  independent  of  the 
ionic  charge,  and  a  function  of  density,  molecular  mass  of  the  gas,  the  dielec¬ 
tric  constant  and  the  reduced  mass-ratio  m*.  For  this  LANGEVIN-LOEB  center- 
of-force  equation,  see  LOEB  (1955).  In  the  same  book  other  derivations  are 
treated,  among  them  that  of  HASSE  and  COOK  (1931),  in  which  the  ionic  charge 
again  appears,  and  the  radius  of  the  ion  constitutes  a  factor. 


The  radius,  or  more  exactly  the  distance  between  the  centers  of  ion  and  mole¬ 
cule  at  impact,  is  a  complicated  parameter,  since  sphere  shapes  for  the  parti¬ 
cles  involved  can  not  be  assumed.  In  fact,  since  mobility  measurements  can 
be  made  rather  accurately,  these  equations  have  been  used  to  determine  these 

center  distances. 

l 


Another  approach  can  be  made  by  starting  from  the  equation  for  conductivity  as 
used  in  ionospheric  physics  (see,  e.  g.  ,  RAWER  (1953)): 


nl  el 

A  =  - 

ni|  v 


(32) 


from  which,  with  the  usual  definition  of  conductivity: 

A  =  rq  etk  (33) 

we  get  for  the  mobility: 

el  1  (34) 

m  j  v 

where  i/is  the  collision  frequency  between  ions  and  gas  molecules.  Thus,  if  the 
collision  frequency  were  known  from  ojher  measurements,  and  we  can  measure 
nj  and  the  ionic  mass  mj  would  be  easily  derived.  Unfortunately,  the  colli¬ 
sion  frequency  is  not  independent  of  the  ion  mass  nij  ,  For  rather  coarse  con¬ 
siderations,  the  application  of  eq.  (32)  will  yield  satisfactory  results,  but  for  the 
determination  of  the  ion  mass  from  a  mobility  spectrum  a  more  accurate  solu¬ 
tion  should  be  approached. 

After  the  fundamental  book  of  LOEB  was  published  (1955),  more  investigations 
on  ionic  mobilities  were  done  by  calculation  and  in  the  laboratory.  These  in¬ 
vestigations  were  often  restricted  to  simple  cases,  if  we  compare  the  situation 
in  the  atmosphere.  Metallic  ions  in  inert  gases,  or  ions  in  their  atomic  parent 
gas  are  examples  for  this  type  of  investigation  (see,  e.  g.  ,  DALGARNO, 
McDOWELL  and  WILLIAMS,  1958;  DALGARNO  1958;  E.  A.  MASON  and  SCHAMP, 
1958).  Since  the  number  of  possible  types  of  ions  is  either  one  or  small,  and 
mostly  known,  and  since  the  neutral  gas  is  not  a  mixture  in  most  of  these  in¬ 
vestigations,  the  application  of  such  results  to  the  atmosphere  poses  some  still 
unsolved  problems. 
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On  the  other  hand,  in  the  region  between  about  60  and  70  km,  there  may  be  a 
possibility  of  measuring  mobility  and  ionic  mass  {without  pumping  down)  at  the 
same  time  and  the  same  locality.  For  mobility  measurements  a  great  number 
of  collisions  within  the  instrument  is  required,  for  mass  spectrometry  of  ions, 
the  number  of  collisions  within  the  instrument  must  be  very  small.  In  princi¬ 
ple,  it  should  be  possible  to  build  a  rather  great  mobility  instrument  and  a  very 
small  mass  spectrometer  (using  very  small  distances  and  high  frequencies  of 
the  control  voltages,  such  as  proposed  by  EIBER  (1965),  and  to  make  simultane¬ 
ous  measurements  at  the  same  locality.  Such  a  combined  measurement  would 
yield  borderline  values  for  the  upper  region,  and  some  extrapolation  to  mobility 
values  in  lower  regions  certainly  is  permissible. 

In  order  to  give  an  indication  on  a  mass-mobility  relation,  we  have  plotted  in 
figure  2  the  dependency  of  the  factor  m*  on  ionic  mass  for  a  neutral  gas  of  29 
amu,  These  curves  hold  for  all  equations,  in  which  mobility  is  a  linear  func¬ 
tion  of  the  square  root  of  m*,  and  in  which  the  parameter  does  not  depend  strong¬ 
ly  on  the  kind  of  ion. 


Figure  2  MASS-MOBILITY  RELATION  ACCORDING  TO  THE  INDICATED  FORMULA  FOR  ol; 
(Abscissa:  molecular  mass  in  amu;  upper  curve:  same  as  lower  curve,  but  with  lineor 
ordinate.) 


1.2.4  CONC LESION 

Knowledge  on  mobilities  of  atmospheric  ions  may  be  used  for  a  variety  of  pur¬ 
poses.  Electric  conductivity  is  determined  in  an  accurate  manner  by  the  num¬ 
bers  of  ions  in  different  mobility  ranges.  A  further  important  application  of 
mobility  is  the  determination  of  the  mean  free  path,  e.  g. ,  the  LANGEVIN- 
BRICARD  mean  free  path  for  ions  (BRICARD,  1964).  Since  the  mean  free  path 
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plays  an  essential  role  in  the  equations  for  recombination  and  attachment  coef- 
ficients  (cf. ,  e.g.,  SIKSNA,  1964),  knowledge  on  mobility  is  in  many  cases  still 
the  most  desired  information. 

Under  equilibrium  conditions,  if  the  combination  coefficients  are  derived  in  a 
quiet  situation,  information  on  the  ionization  rate  in  disturbed  situations  may  be 
obatined  for  different  heights. 

Information  on  diffusion  and  diffusion  coefficients  is  another  result  to  be  ex¬ 
pected  from  the  measurement  of  the  ion  spectrum. 

Of  course,  the  total  ion  number  densities  is  an  immediate  consequence  of  the 
ion  spectrum.  Variations  in  this  number,  which  occur  upon  descent  through  the 
atmosphere,  may  indicate  changeover  points  between  diffex-ent  kinds  or  groups 
of  ions.  One  of  such  changeover  points  may  concern  the  variation  from  mole¬ 
cular  to  "cluster"  ions,  and  this  knowledge  may  yield  a  contribution  to  the  study 
of  the  ion  nature  in  lower  heights. 

For  given  combinations  of  the  involved  parameters,  the  different  theories  con¬ 
cerning  the  mobility-mass  relation  may  be  compared,  and  it  may  be  determined 
to  which  degree  of  resolution  and  accuracy  they  may  yield  the  same  result. 

However,  if  such  accuracy  is  not  sufficient  in  certain  cases,  and  even  if  it  is 
not  known  which  theory  should  be  applied  and  how  it  should  be  applied,  it  may 
still  be  possible  to  derive  a  much  more  accurate  result  by  a  combination  of  dif- 
erent  approaches.  These  would  involve  the  results  of  laboratory  investigations 
on  ionic  mobilities  in  cases  of  predetermined  gases  and  ions. 

Finally,  the  possibility  of  carrying  oui  simultaneous  measurements  of  ionic  mo¬ 
bility  and  ionic  mass  by  a  common  carrier,  which  seems  to  be  within  reach  in 
the  next  few  years  for  a  limited  range  of  altitude,  and  the  further  possibility  of 
extrapolating  these  results  downward  j to  a  certain  extent,  and  upward  for  a  com¬ 
parison  with  other  measurements  executed  in  higher  ranges,  may  help  consider¬ 
ably  to  close  the  gap  of  knowledge  in  these  mass-mobility  relations  and  to  es¬ 
tablish  a  solid  basis  for  statements  on  the  number  densities,  chemical  natures, 
masses,  and  mobilities  of  atmospheric  ions. 

it  may  be  said,  that  mobility  is  the  most  essential  measuring  parameter  con¬ 
cerning  ‘ons  in  the  upper  stratosphere  and  mesophere  of  the  Earth  or  in  similar 
environments  somewhere  else. 
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1.3  MEASUREMENT  OF  IONIC  MOBILITY 


1. 3.  1  THE  MEASURING  PRINCIPLE  OF  THE  PERPENDICULAR  FORCES 

If  a  gas  containing  ions  is  drawn  through  an  electrical  condenser  of  capacity 
C  and  with  a  voltage  difference  U  between  its  plates,  a  current  I  which  is  a 
measure  of  the  time  constant  of  the  gas  according  to  DOLEZALEK  (1956)  may  be 
measured  at  either  one  of  the  plates: 

A  I/U  (35) 

c  C 

where 

A  is  the  conductivity  of  the  gas 

f  is  the  capacitivity  of  the  gas  (permittivity,  dielectric  constant). 

Since  I/U  is  the  conductance  of  the  condenser,  this  equation  seems  almost 
to  be  an  identity.  In  fact,  it  is  the  basic  equation  for  the  method  of  the  perpend¬ 
icular  forces  for  the  measurement  of  conductivity  and  mobility.  This  method 
is  truly  a  method  to  measure  the  time  in  which  an  atmospheric  ion  travels  a 
certain  given  and  known  distance  .ider  the  influence  of  an  electric  field  (one 
force).  By  the  second  force,  th»  "wind  speed"  in  the  condenser,  this  time 
measurement  is  transferred  into  measurement  of  a  length. 

If  the  current  quoted  in  eq.  (35)  is  measured  at  the  negative  electrode  of  the 
condenser,  the  conductivity  A  in  the  time  constant  e/X  is  the  positive  polar 
conductivity,  and  vice  versa,  and  the  total  conductivity  is  the  sum  of  both: 

A  =  A+  +  A~  ;  .  (36) 

and  in  all  cases,  only  ions  with  a  mobility  smaller  than  a  critical  value  are 
measured  in  their  full  number  (a  more  detailed  discussion  of  this  will  be  given 
in  section  1 . 3. 2,  ,  below). 

The  first  quotation  of  this  method  is  found  in  RUTHERFORD'S  paper  of  1899. 
ZELENY  and  LENARD,  both  in  1900,  and  KA'HLER  in  1903  applied  about  the 
same  principle,  before  GERDIEN  turned  it  (1905)  into  the  most  common  tool 
for  the  measurement  of  the  polar  electric  conductivities  in  the  atmosphere. 

Since  then,  the  condenser  as  used  here  --  mostly  a  cylindrical  one  --  is  often 
referred  to  as  the  "GERDIEN-condenser,  "  "GERDIEN- chamber,  "  "GERDIEN-  . 
counter,  "  or  simply  the  "GERDIEN.  " 
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1. 3.  2  THE  "CRITICAL  MOBILITY"OF  A  GERDIEN  CONDENSER* 

The  "critical  mobility"  of  a  GERDIEN  condenser  is  not  an  atmospheric  para¬ 
meter;  it  is  an  instrumental  feature  of  a  certain  GERDIEN,  depending  on  the  type 
and  geometrical  configuration  of  it,  and  on  the  applied  "driving"  voltage  and  the 
applied  wind  speed  in  the  condenser. 

To  define  the  "critical  mobility,  "  we  consider  a  plate  condenser  through  which 
air  is  sucked,  and  we  investigate  the  motion  of  the  ions  under  the  applied 
electric  field  (driving  voltage  divided  by  electrode  distance)  and  the  velocity 
of  the  air. 

From  figure  3  we  derive  directly: 


x  s  L 


*-x 
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Figure  3  DERIVATION  OF  CRITICAL  MOBILITY  OF  A  GERDIEN 


movement  of  ions  in  horizontal  direction: 
dx 

—  _  v 
dt 


(37) 


X  Vt  •  }Q  ; 

and  in  vertical  direction: 


(38) 


-  k  E  <39> 

dt 


>'  -  k  K  t  ,  y0  ; 


(40) 


where  !•  is  the  electric  field,  k  the  mobility,  v  the  velocity  of  the  air,  and 
xq  and  y0  are  integration  constants. 


The  combination  of  cqs.  (38)  and  (40)  gives 


foil  nvniR  (Ik  derivation  as  it  is  given  by  ISRAEL  (10S7),  pp  108  ff. 
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y  -  y  o 


kE 


{41) 


X  -  X0 


This  is  the  equation  for  the  ion  paths  lto  6  shown  in  figure  3.  All  ions  entering 
the  condenser  will  reach  the  lower  electrode,  if 


y  =  0  for  x  =  0  I 

y  =  d  for  x  =  L  j  (42) 

where  d  is  the  distance  between  the  condenser  plates,  and  L  the  length. 
With  eq.  (41),  we  get  the  "limiting1’  or  "critical"  mobility  as  follows: 


k 


c 


(43) 


or,  by  introducing  the  capacitance  of  the  condenser  C, 


c  "  cu 


(44) 


where  U  =  Ed  ,  and  M0  the  volume  of  air  which  is  sucked  per  second  through 
the  condenser. 

All  ions  which  have  a  mobility  equal  to  or  greater  than  kc  ,  are  reaching  the 
lower  electrode  and  thus  measured.  Ions  with  a  smaller  mobility  may  pass 
the  condenser  without  being  captured,  if  they  enter  at  a  value  of  y  which  is 
small  enough.  In  fact,  of  ions  with  a  mobility  k  <  kc,  just  the  ratio  k/kc 
is  captured.  As  can  be  seen  from  figure  3,  all  ions  which  enter  at  or  below 
the  solid  line  with  the  number  4  are  captured.  They  are  entering  at  an  y- 
value  p-d,  where  p  is  a  number  smaller  than  unity.  Our  conditions  eq.  (42) 
are  now  modified  to: 

y  =  pd  for  x  =  0  j 

y  =  d  for  x  =  L  j  (45) 

and  the  critical  mobility  k'  for  the  condenser  with  the  distance  d  (1  -  p)  be¬ 
tween  the  electrodes  comes  to: 

^  =  d(l-p)  •  (46) 
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Comparing  eqs.  (43)  and  (46),  we  see  that: 


k^/kc  =  1  -  p  . 


(47) 


Since  it  can  easily  be  shown  that  the  volume  of  air  passing  through  this  smaller 
condenser,  M'  ,  is  related  to  the  value  M  according  to: 

M'/M0  =  1  -  p,  (48) 

and  since  this  derivation  for  k'  may  be  repeated  for  any  othe-  of  k  <kc  , 

our  statement  is  verified  that  ions  with  mobility  k  <  kc  are  meaa  .....  d  with  the 
ratio  k/kc  of  their  number. 

1.  3.  3  APPLICATION  OF  THE  GERDIEN  TO  THE  MEASUREMENT  OF  THE 
ION  SPECTRUM 


The  cylindrical  condenser  as  proposed  by  GERDIEN  (1905),  i.  e. ,  with  one 
"driving"  and  one  "receiving"  electrodes  (the  latter  one  being  the  electrode  to 
which  the  particular  species  of  ions  which  is  to  be  measured  is  driven  by  the 
applied  voltage),  and  with  ions  entering  across  the  full  opening  of  the  cylinder, 
can  be  used  for  the  measurement  of  the  ion  spectrum.  To  show  this,  we  follow 
the  derivation  of  ISRAEL  (1931). 


The  mobility  spectrum  is,  by  definition: 


dn 

7k 


f(k); 


(49) 


where  dn  is  the  number  of  ions/m^  in  the  mobility  range  between  k  and  k  +  dk ; 
k  i  s  the  mobility,  and  f(k)  is  a  function  of  the  mobility. 

• 

If  ions  of  different  mobilities  are  intercepted  at  the  "receiving  electrode"  of  the 
aspiration  condenser,  the  current  represented  by  these  ions,  is  generally  given 
by 


I  -  M0 


k  f(k)  dk  +  Mq  e 


I  f(k)dk; 


(50) 


where  M0  is  the  amount  of  air  passing  through  the  aspiration  condenser  (m  /sec) 
e  is  the  charge  per  ion  (i.  6  x  10" *9  cb);  andkcis  the  "critical  mobility"  of  the 
condenser  (see  section  1.3.2)  , 
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If  mobility  is  entered  in  (m/sec)  /  (V/m),  the  current  in  eq,  (50)  becomes  in  A. 
If  we  want  to  derive  the  mobility  spectrum  from  the  GERDIEN-characteristic, 
i.e..from  the  current  versus  voltage  curve  of  the  GERDIEN  condenser,  we  are 
looking  for  a  relation  between  f(k)  and  the  parameters  "current"  I  and  "voltage" 
U  of  the  GERDIEN  condenser.  In  eq.  (50)  we  have  a  relation  between  I  and  kc. 
Another  relation,  that  between  kc  and  U.may  be  derived  as.  follows. 

The  critical  mobility  is  a  constant  of  the  measuring  apparatus,  the  GERDIEN 
condenser.  It  is  defined  as  (see  eq.  (44) 


k 


C 


(51) 


where  C  is  the  effective  capacity  of  the  condenser  (i.e.  ,  the  capacity  without 
consideration  of  the  capacity  of  the  connected  wiring,  supports,  etc.  ),  and  ( 
is  the  capacitivity  of  the  air  (5.  86  x  F/m).  By  differentiation  we  get: 


dkc 

dU 


(52) 


Combining  eq.  (50)  after  differentiation,  and  eq.  (52)  we  get 


for  the  second  differentiation,  we  apply 

d2I  _d_  /dl\  ^ 

=  dTc  \m)  dU 

and  get: 

d2I  efM0 

-  =  -  -  f(k)  . 

dU2  CU3 


(53) 


(54) 


(55)* 


*  This  derivation  was  first  demonstrated  hy  H.  BENNDORF.  For  the  differentiations  the  rule  of  LEIBNIZ  has  been 
applied,  namely,  foi  the  differentiation  of  eq.  (50)  in  the  form: 

d  /*y  ^  d  F  (x,  y) 

I  F(x,y)dx  =  I  “  F  F(y,y); 

dy  I  J  d  y 

Jv  =  A  •'v  =  Cl 


and  for  the  differentiation  of  eq.  (53),  with  eq.  (54),  in  the  following  form: 


o 

Ty 


F  (x)  dx  =  F  (y)  . 


0 
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We  may  replace  here  kc  by  k  ,  if  we  are  doing  the  measurement  in  such  a  way 
that  the  particular  k  under  consideration  is  always  just  the  kc  and  this  is  done 
by  proper  adaptation  of  U.  Thus,  we  arrive  at: 


dn  C  U3  d2 1 

f(k)  =  —  =  -  - r  — - 

dk  erM2  dU2 


(56) 


or, 


specified  for  a  cylinderical  condenser: 

dn  2LU3  d2 1 

dk  /rev2  (R2  —  r2)2  In  R/r  dU2 


(57) 


with  L  as  the  length,  and  R  and  r  as  the  radii  cf  the  outer  and  the  inner  electrodes, 
respectively,  and  In  being  the  logarithm  to  the  basis  2.  71828,  . , 

It  is  shown  by  eqs.  (56)  and  (57),  that  the  ion  spectrum  can  by  measured  by  the 
GERDIEN  condenser,  if  the  applied  driving  voltage  U  is  varied  ("swept")  while  the 
current  from  the  receiving  electrode  is  measured,  so  that  the  secon  derivative 
may  be  calculated  of  taken  from  the  I  over  U  curve. 

1.  3.4.  GENERAL  OPERATING  CONDITIONS  OF  THE  GERDIEN  CONDENSER 

A  certain  number  oi  conditions  must  be  fulfilled  if  the  measurements  made  by 
a  GERDIEN  condenser  are  to  give  correct  values. 

The  life  time  of  the  ions  to  be  measured  must  be  so  long,  that  the  ion  number 
density  does  not  change  due  to  combination  or  attachment  to  other  particles 
while  the  ions  are  traveling  inside  the  GERDIEN  condenser.  The  number  of 
ions,  as  compared  with  the  number  of  neutral  gas  molecules,  and  the  average 
distances  between  the  ions  must  be  neglected  as  compared  with  the  electric 
field  applied  io  the  chamber.  There  must  be  no  space  charges  of  a  magnitude 
so  as  to  compete  with  the  applied  fiehi. 

The  number  density  of  ions  inside  the  chamber  must  be  the  same  as  it  is  out¬ 
side  the  chamber.  That  means  that  no  ions  which  approach  the  chamber  in  a 
cylindrical  volume  of  air  of  the  diameter  of  the  orifice  of  the  chamber  are  inter¬ 
cepted  at  the  entrance  or  moved  to  a  location  outside  the  chamber;  in  other 
words,  the  so-called  "edge -effect"  must  be  avoided,  The  edge  effect  usually 
is  created  by  electric  stray  fields  at  the  front  end  of  the  condenser,  generated 
by  the  applied  driving  voltage  or  by  potential  differences  between  the  in¬ 
strument  and  its  environment.  This  effect  has  been  studied  a  great  deal  since 
1905. 

The  density  of  the  air  inside  the  chamber  must  be  the  same  as  outside  of  it. 
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That  means,  that  the  chamber  must  not  provide  a  remarkable  air  flow  resistance, 
and  there  must  be  no  deviation  in  the  direction  of  the  relative  air  flow  in  front 
of  the  chamber. 


The  air  flow  in  the  chamber  must  be  a  laminar  one,  i.  e.  ,  it  must  not  be 
turbulent.  Also  this  effect  has  been  studied  by  a  number  of  authors. 

Since  edge  effectand  turbulence  disturbances,  in  the  general  meaning,  do  not 
play  a  role  in  the  case  of  the  particular  "modified  GERDIEN  chamber"  to  be 
particularly  discussed  in  this  report,  we  refer  to  the  literature  as  quoted  in 
the  textbooks  on  atmospheric  electricity,  and  do  not  discuss  these  effects  here. 

The  possibility  of  disturbances  given  by  any  motion  of  ions  due  to  other  forces 
than  electrical  and  convection  with  the  air  motion  in  the  chamber  is  already 
indicated  in  eq,  (1)  (section  1.2.1);  in  particular,  diffusion  is  to  be  mentioned 
here,  Also  this  is  discussed  extensively  in  the  literature  ,  and  we  shall  con¬ 
sider  this  influence  later  on. 

At  the  same  place  (i.  e. ,  in  chapter  1.7  ),  some  more  limitations  and  disturb¬ 
ances,  which  occur  with  measurements  of  ion  number  densities,  conductivities 
and  mobilities  in  the  upper  stratosphere  and  in  the  mesosphere,  will  be  invest¬ 
igated. 

Almost  none  of  the  conditions  quoted  here  can  be  fulfilled  with  a  perfect  accuracy. 
Thus,  the  influences  of  the  different  deviations  must  be  weighted  to  determine 
the  accuracy  of  a  measurement. 

i.  3.  5  VARIATIONS  IN  THE  DC-OPERATED  GERDIEN  CHAMBER 

In  section  1.  3.  3,  we  have  shown  that  the  ion  spectrum  may  be  derived  from  a 
"usual"  GEP.DIEN  chamber,  if  the  driving  voltage  U  is  varied  ("swept")  and 
the  measured  current  is  differentiated*twice  to  this  voltage.  "Usual"  means 
that  the  cylindrical  condenser  has  only  the  two  electrodes,  and  both  are  not 
subdivided,  and  that  ions  are  taken  in  through  the  full  opening  of  the  cylinder. 

Since  a  double  differentiation  generally  reduces  the  accuracy  obtainable  with  a 
measurement,  means  of  developing  instruments,  which  would  do  with  only  one 
or  even  without  any  differentiation  have  been  investigated. 

A  classification  of  these  varied  GERDIEN  chambers  has  been  proposed  by 
TAMMET  (1960).  He  calls  the  "usual"  GERDIEN  chamber  an  "integrated 
chamber".  Next  he  quotes  two  different  kinds  of  "differential  chambers  of  the 
first  order".  In  both  cases,  the  mobility  spectrum  is  derived  by  sweeping  the 
driving  voltage  and  differentiating  the  measured  current  but  once.  Finally, 
there  is  the  "differential  GERDIEN  chamber  of  the  second  order",  with  which 
the  measured  current  gives  the  mobility  spectrum  directly.  In  figure  4,  the 
principle  of  these  four  varieties  is  shown  in  the  example  of  the  plate  cond^r-cer. 
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The  same  principles  of  variation  apply  to  cylindrical  (and  other)  condensers: 
while  the  integrated  chamber  allows  the  ions  to  enter  through  the  full  crosss 
section,  and  has  them  intercepted  all  at  the  same  "receiving  electrode,  11  and 
while  the  differential  chamber  of  the  second  order  lets  the  ions  pass  only 
through  a  nari  -\v  slit  and  makes  the  interception  at  a  subdivided  receiving  e- 
lectrode,  the  two  kinds  of  differential  chambers  of  the  first  order  apply  either 
one  of  these  two  additional  features. 

The  analytical  expression  for  the  mobility  spectrum  for  the  integrated  chamber 
has  been  derived  above,  in  section  1.  3.  3.  Specified  for  the  plate  condenser, 
it  is  entered  in  figure  4. 

To  derive  the  mobility  spectrum  formula  for  the  dc-operated  differential  cham¬ 
ber  of  the  second  order,  we  start  from  figure  5,  for  both  the  plate  and  the  cy¬ 
lindrical  condenser  form. 


DRIVING  ELECTRODE 


DRIVING  ELECTRODE 
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F*9u.e  5  SCHEMATIC  FIGURE  FOR  THE  DERIVATION  OF  THE  FORMULA  FOR  THE  MOBILITY 
SPECTRUM  FOR  THE  DIFFERENTIAL  GERDIEN  CHAMBER  OF  THE  SECOND  ORDER 
(o  -  plofe  condenser,  b  cylindrical  cona'enser.) 


Since,  in  a  differential  ion  counter  of  the  second  order,  ions  are  entering  only 
at  a  predetermined  slit,  the  calculation  of  the  mobility  spectrum  is  rather 
simple.  The  current  I  is  given  by  the  number  density  of  ions  times  their  elec¬ 
tric  charge  times  the  product  (Ay  bv)  which  is  the  numbers  of  unit  volumes  en¬ 
tering  per  time  unit  through  the  slit  \y  is  the  width  of  the  slit  or  intake  open¬ 
ing,  b  is  the  width  of  the  condenser,  and  v  is  the  air  velocity): 
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I  =  n  Ay  b  v; 


or,  written  for  n 


e  Ay  b 


Because  of  the  obvious  relation: 


dx  /  dt 
dy  j  dt 


k  E  k  U 


eq.  (61)  can  be  written  as  follows: 


n  =  -  I; 

e  Ay  b  k  E  Xj 


and  then  be  differentiated  to  get  the  mobility  spectrum: 


VO  d  b  e  Ay 

(not  noting  the  sign) 
which  is,  with 

v  b  d  = 


v  b  Ay  =  Mq 

the  equation  for  the  mobility  spectrum  as  entered  in  figure  4: 


dk  yj  eM0  Mg 

(67) 

By  the  application  of  any  one  of  the  many  possibilities  for  substitution  given 
with  eq,  (62),  a  number  of  modified  expressions  for  eq.  (67)  may  be  derived, 
and  this  may  be  helpful  for  an  interpretation  of  this  mobility  spectrum  formula. 
Considering: 


Mq  =  v  Ay  2i7(yj  +  r) 


(68) 


(compare  figure  5)  and  replacing  eq.  (61)  by 


ve  Ay  2/r  (y  j  +  r  ) 


and  eq.  (62)  by 
2 

yr  i  2  uk 

Xj  v  In  R/r 

we  derive,  in  a  similar  way,  the  formula  for  the  cylindrical  condenser: 


I  ?  e  Ay  y,  (y,  *r)  In  R/r 

I  1  1-1  (71) 

which  has  been  written  here  in  a  form  comparable  to  eq,  (67),  and  which  may 
be  modified  in  various  forms  by  substitutions  such  as  given  in  eq.  (70). 

In  figure  4,  the  corresponding  formulas  for  the  differential  GERDIEN  of  the 
first  order  are  quoted,  namely,  for  type  a  (i.  e.  ,  full  intake,  subdivided  re¬ 
ceiving  electrode): 


and  for  type  b  (i.  e.  ,  subdivided  intake,  full  receiving  electrode: 


dk  y,  e  Mg  M^  dll 

(73) 

The  detailed  derivation  of  these  equations  is  not  given  here. 

Historically,  the  differential  ion  chamber  of  the  first  order  has  already  been 
applied  by  BLACKWOOD  in  1920,  and  the  differential  chamber  of  the  second 
order  was  used  by  ERIKSON  in  1921  (ERIKSON  1921,  1922,  1924,  1929 
MAHONEY  1929,  and  Seville  CHAPMAN  1937). 

Recent  reports  on  the  measurement  of  ion  spectra  have  been  published  by 
MISAKI  especially  for  large  ions  close  to  ground  (1950,  1961,  1964);  TSVANG 
(1956);  E,  WHIPPLE  (I960);  and,  for  ions  in  the  troposphere,  by  HOPPEL 
(1963)  and  HOPPEL  and  KRAAKEVIK  (1964).  Disturbances  occurring  in 
GERDIEN  chambers  have  been  discussed  recently  (aft  jr  many  others  since 
1905)  by  KOMAROV  (I960  a  and  b)  and  TAMMET  (1962).  The  paper  by  TAMMET 
(1960)  in  which  the  chambers  are  classified,  has  already  been  mentioned. 


Surveys  on  methods  to  measure  ionic  spectra  are  given  by  LOEB  (1949,  1955) 
and  by  KNOLL,  EICHMEIER,  and  SCHON  (1  964,  pp.  198-205). 


-31- 


1.  3.  6  RESOLUTION  OF  THE  DC -OPERATED  GERDIEN  CHAMBERS 


An  important  problem  is  given  by  the  fact  that  the  resolution  of  all  types  of 
GERDIEN  chambers  decreases  with  increasing  mobility.  By  differentiation 
of  eq.  (62),  we  get 

v  yj  d 

~u 

(74) 

for  the  resolution.  The  driving  voltage  u  cannot  be  decreased  beyond  a  level 
which  is  given  by  the  unavoidable  noise.  This  noise  is  due  mostly  to  the  VOLTA 
effect.  Since  the  condenser  has  to  be  exposed  to  the  free  atmosphere,  pollu¬ 
tions  of  the  surfaces  of  the  electrodes  and  other  differences  between  the  elec¬ 
trode^  can  hardly  be  avoided,  and  the  VOLTA  effect  plays  a  role  (compare 
ISRAEL  and  DOLEZALEK,  1957).  The  smallest  applicable  voltage  will  be  be¬ 
tween  1  and  5  V,  and  it  is  advisable  to  stay  well  over  that. 

The  wind  velocity  v  should  not  be  increased  beyond  the  speed  of  sound,  be¬ 
cause  otherwise  the  shock  wave  may  bring  disturbances.  The  geometrical 
dimensions  in  radial  direction  yj  and  d  have  certain  limits,  too.  Finally, 
the  length  Ax  must  be  kept  great  enough  to  allow  for  the  disturbances  due  to 
diffusion. 

Obviously,  this  is  the  difficulty  which  explains  why  no  direct  measurements 
of  conductivity,  let  alone  mobility,  have  been  presented  for  altitudes  higher 
than  abcut  30  km.  If  we  want  to  extend  the  measurements  into  higher  ranges, 
either  v  or  U,  or  both  parameters  in  eq.  (74)  should  be  replaced  by  another 
one.  In  other  words,  a  new  measuring  principle  must  be  found. 
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1.  4  THE  MODIFIED  GERDIEN  CHAMBER 


l.  4.  1  THE  MEASURING  PRINCIPLE 

It  has  been  said  above  (section  1.3.  1),  that  the  measuring  method  of  the  perpen¬ 
dicular  forces  is  essentially  a  measurement  of  the  time  in  which  an  ion  travels 
over  a  given  distance.  Since,  for  great  mobilities,  the  resolutions  become  too 
small  to  allow  a  successful  measurement  with  instruments  of  practical  size  and 
applicable  voltages  and  air  speeds,  it  seems  reasonable  to  try  to  enlarge  that 
"given  distance"  without  changing  the  diameter  of  the  chamber.  In  order  to  ver¬ 
ify  this,  one  or  more  of  the  parameters  in  eq.  (74)  must  be  replaced  by  another 
one. 

There  are  a  number  of  possibilities  for  this  approach.  One  of  them  is  the  re¬ 
placement  of  U  by  a  parameter  u/  «  where  <o  is  the  circular  frequency  of  an 
ac  voltage  -  in  other  words,  to  replace  the  dc  driving  voltage  discussed  so  far 
in  this  report  by  an  ac  driving  voltage.  It  will  turn  out  that  by  introducing  the 
ac  voltage,  the  wind  velocity  does  not  influencing  the  resolution  anymore.see 
below. 

Here. a  word  should  be  said  with  respect  to  our  application  of  the  denomination 

ac.  In  the  past,  GERDIEN  chambers  have  been  proposed  which  operate  with  an 
ac  voltage,  but  in  these  cases,  the  frequency  was  so  small  that  for  the  movement 
of  a  particular  ion  on'  a  dc  field  was  effective.  If  one  of  the  GERDIEN  chambers 
quoted  in  section  1.  3.  5  should  be  used  for  ions  of  both  signs,  and  with  a  sweep¬ 
ing  driving  voltage,  the  solution  of  the  application  of  a  low  frequency  ac  is  a 
natural  one.  However,  in  all  these  cases,  the  length  of  the  ion  motion  is  not 
essentially  changed,  and  the  equations  of  section  1.  3.  5  and  i.  3.  6  basically 
apply. 

If  we  want  to  increase  the  length  of  the  ion  paths  within  the  condenser  by  the  ap¬ 
plication  of  an  ac  driving  voltage,  the ’applied  frequency  must  be  so  great  that 
the  ions  follow  a  sinusoidal  path.  This  principle  was  proposed  by  DOLEZALEK 
in  1962  and  discussed  by  DOLEZALEK  (1963)  and  by  DOLEZALEK  and  OSTER 
(1963,  1964,  1965). 

Figure  6  shows  the  principle  and  figure  7  demonstrates  the  difference  in  reso¬ 
lution  for  dc  and  ac  chambers.  The  lower  part  of  figure  6  shows  the  influence 
of  a  superimposed  dc  voltage  of  a  magnitude  comparable  to  the  amplitude  of  the 
ac  driving  voltage,  and  it  indicates  how  this  influence  can  be  minimized.  As 
can  be  seen,  the  ions  are  forced  to  follow  a  sinusoidal  pattern  inside  the  cham¬ 
ber,  the  amplitude  of  which  is  constantly  increasing  while  the  ions  travel  down¬ 
stream.  The  increase  in  the  amplitude  of  the  driving  voltage  along  the  chamber's 
axis  is  done  by  a  rather  fine  subdivision  of  the  driving  electrode,  and  by  a  net¬ 
work  of  relatively  small  resistances  every  partial  driving  electrode  is  supplied 
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Figure  6  ION  MOVEMENT  IN  A  PLATE  CONDENSER  WITH  AN  AC  VOLTAGE  INCREASING 
FROM  LEFT  TO  RIGHT,  AIR  FLOW  FROM  LEFT  TO  RIGHT 
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(A)  D.C.  DRIVING  VOLTAGE! 
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Figure  7  ION  PATHS  IN  DIFFERENTIAL  GERDIEN  CHAMBERS  OF  THE 

SECOND  ORDER 

(Two  ions  with  mobilities  k^  and  k2  enter  from  the  ion  gate  into  the 
chamber  at  the  same  time.) 


with  a  slightly  greater  ac  amplitude  than  its  upstream  neighbor.  This  increase 
must  be  done  by  small  resistors  to  avoid  phase  delays,  and  it  must  be  done 
smoothly  to  avoid  any  phase  shifts  for  the  ions  traveling  downstream. 

The  actual  ion  paths  for  a  given  cylindrical  condenser  have  been  calculated: 


/  2  k  U 

/'  =  / — ~  —  sin  (ojt  t  <3)  •  Cj 

o  In  R  r 

Here,  is  the  radial  location  of  the  ion  at  the  time  t,  <1>  is  the  phase  of  the 
sinusoidal  driving  voltage,  and  Cl  is  an  integration  constant.  The  significance 
of  the  other  letters  used  is  as  before  (see  figures  4  and  5);  however,  U  is  now 
"amplitude"  of  the  driving  voltage,  i.  e.  ,  the  zero -peak  potential  difference  and 
o>  is  the  circular  frequency  of  it. 

Figure  8  shows  an  example  for  these  ion  paths,  demonstrated  for  ions  of  one 
mobility  only,  but  entering  at  different  times  with  respect  to  the  zero  point  of 
the  driving  voltage.  It  can  be  seen  from  this  figure  that,  for  example,  an  ion 
entering  at  c  -  *>  S  T  would  not  reach  a  receiving  electrode  placed  at  r  =  ]. 5 
cm  during  the  first  sine  wave  it  follows  in  the  chamber,  but  because  of  the  in¬ 
creasing  U,  it  may  intercept  during  one  of  the  following  waves.  However,  only 
ions  with  a  smaller  mobility  are  supposed  to  intercept  more  downstream.  For 
this  reason,  only  ions  which  arrive  at  t  =  o  will  be  allowed  to  enter.  This 
poses  a  new  technical  problem  for  the  realization  of  the  method. 


Figure  3  ION  PATHS  IN  A  CYLINDRICAL  CONDENSER  WITH  AN  A.C.  DRIVING  VOLTAGE. 
IONS  OF  ONE  MOB  I  LHY  ONLY.  IONS  ENTERING  AT  THE  SAME  ORIFICE  (ION  GATE)  AT 
DIFFERENT  MOMENTS,  NAMELY  AT  OT,  I  AT,  1  2T,  5  6T,  3  4T.  T:  DURATION  OF  ONE 
PERIOD  OF  THE  SINE  WAVE. 
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By  positioning  a  "collecting  electrode  "  -  in  the  case  of  figure  8  -  at  about  2.  6 
cm  and  with  a  length  reaching  downstream  to  about  T/8,  all  ions  with  a  phase 
other  than  zero  would  be  caught  there.  This  solution  of  the  problem  is  applica¬ 
ble  in  the  case  of  known  wind  speed,  since  the  abscissa  in  figure  8  is  time  and 
the  length  of  the  collecting  electrode  must  be  chosen  according  to  velocity.  The 
collecting  electrode  must  not  be  longer  than  corresponding  to  T/8  because  of  the 
increase  of  driving  voltage  along  the  axis. 

Another  method  of  solving  this  problem  is  given  by  the  possibility  of  making  the 
ion  orifice  to  an  ion  gate,  which  is  closed  all  the  time  with  the  exception  of  a 
short  duration  at  t  =  O  in  figure  8.  This  can  be  achieved  by  a  particular  gen¬ 
erator  triggered  by  the  driving  voltage.  Figure  9  shows  its  principle.  At  the 
upstream  end  of  the  modified  GERDIEN  chamber,  there  are  three  sets  of  rings, 
the  purpose  of  which  will  be  discussed  in  the  next  figure.  The  ion  gate  is  lo¬ 
cated  between  rings  13  and  14.  Between  these  rings,  a  dc  voltage  which  is  great 
enough  to  capture  all  ions  of  the  expected  mobility  range  is  maintained.  Only 
during  the  time  t  this  dc  voltage  is  removed  and  ions  will  pass.  Since  the 
rings  03,  04,  14,  23,  and  24  are  always  on  the  same  potential,  ions  will  not  meet 
an  electric  field  while  passing  between  03  and  04  and  while  passing  between  23 
and  24.  (All  actual  rings  and  their  numbers  in  this  example  refer  to  the  labor¬ 
atory  model  of  the 'hnodified  GERDIEN  chamber'1,  the  so-called  "Model  A," 
which  will  be  described  in  Part  2  of  this  repoit). 

The  purpose  of  the  rings  11  to  18,  as  indicated  in  the  upper  left  part  of  figure  9, 
is  to  allow  the  air  to  pass  between  them,  but  to  remove  all  ions  out  of  it.  To  do 
this,  dc  voltages  are  maintained  between  all  these  rings.  This  same  set  of  rings 
is  to  be  seen  in  figure  10  as  well;  in  the  same  figure  we  also  see  the  two  other 
sets  of  rings,  parts  of  which  are  already  shown  in  the  upper  left  part  of  figure 
9.  The  purpose  of  these  two  sets  of  rings  is  explained  below. 

The  first  set  of  rings,  numbers  00  to  07,  are  to  cancel  the  so-called 
edge  effect.  Because  of  the  strong  dcVoltages  to  be  maintained  between  the 
pairs  of  rings  of  the  set  10  to  17,  stray  fields  would  occur  going  across  the  in¬ 
let  to  the  pair  13-14.  The  effect  of  these  stray  fields  would  be  to  reject,  more 
or  less,  ions  of  one  sign,  and  to  promote  ions  of  the  other.  Since  this  would 
falsify  the  number  densities,  the  rings  00  to  07  had  to  be  added.  They  are  all 
on  the  same  potential  as  ring  14,  Thus,  during  the  time  r  (figure  9)  when  there 
is  no  potential  difference  between  13  and  14,  ions  entering  between  03  and  04  will 
feel  no  electric  field  while  passing.  The  same  had  to  be  done  at  the  downstream 
side  of  rings  10  to  17.  However,  this  side  faces  the  interior  of  the  chamber  and 
there  is  a  field.  For  this  reason,  the  rings  of  the  set  20  to  27  are  all  carried 
on  different  potentials.  Each  ring  gets  the  potential  which  a  point  with  the  same 
radius  as  the  ring  would  have  in  an  infinitely  long  cylindrical  condenser.  This 
is  done  with  one  exception:  rings  23  and  24  both  have  the  same  potential,  namely, 
that  of  the  point  with  the  radius  ( r2 j  +  r 24)  / 2  in  an  infinity  long  cylinder.  The 

consequence  of  this  is  that  an  ion  entering  between  03  and  04  will  "see"  no  elec- 
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Figure  9  OPERATION  OF  ION  GATE  DEMONSTRATED  WITH  THE  LABORATORY 
MODEL  OF  MODIFIED  GERDIE.N  CHAMBER  ("MODEL  A") 

Upper  left:  three  pairs  of  the  sets  of  rings  at  the  upstream  end  of  modified 
GERD1EN  chamber  with  the  ring  numbers  as  indicated  in  figure  10. 

Upper  right:  Schematic  view  of  the  rings  as  seen  when  looking  parallel  to  the 
axis  of  the  chamber. 

Lower  left:  Variation  of  voltage  versus  ground  at  rings  13  and  14  with  time;  T 
is  the  period  of  the  driving  voltage. 

Lower  right:  The  three  conditions  for  the  operation  of  the  ion  gate:  (I)  is  the 
condition  for  applying  formulas  of  the  plate  condenser;  (II)  is  the  condition 
for  capturing  all  ions  of  mobility  k  during  the  closing  time  of  the  gate;  and 
(III)  is  the  condition  for  the  passage  of  ions  of  mobility  k  during  the  opening 
time  of  the  gate.  N  is  the  total  number  of  ions  of  mobility  k  which  will  pass 
during  one  gate  opening  •  ,  when  the  ion  number  density  prior  to  entering  the 
gate  is  n. 
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trie  field  prior  to  entering  the  chamber  proper,  and  this  is  what  must  be  reques¬ 
ted.  To  avoid  any  fields  outside  the  chamber  for  the  laboratory  experiments 
done  in  the  wind  tunnel,  the  wind  tunnel  is  lined  by  an  electrically  insulated 
sleeve,  and  this  is  carried  along  at  the  same  potential  as  the  rings  00  to  07  have. 
The  beginning  of  this  sleeve  is  shown  at  the  left  hand  end  of  figure  10. 

The  rings  20  to  27,  furthermore,  serve  to  equalize  the  electric  field  in  the  neigh¬ 
boring  part  of  the  chamber,  which  otherwise  would  show  the  wall  known  effect 
at  the  end  of  a  cylindrical  condenser.  For  the  same  reason,  another  set  of  rings, 
30  to  37,  is  located  at  the  downstream  end  of  the  chamber  (figure  10). 

The  potentials  of  all  these  rings  are  functions  of  the  driving  voltage.  This  is  to 
be  seen  from  the  fact  that  the  "undisturbed  potentials"  at  the  locations  of  rings 
20  to  27  (and  30  to  37)  are  determined  by  the  potential  of  the  adjacent  driving 
electrode  (which  is  No.  60  and  No,  75  in  the  case  of  figure  10,  where  the  inner 
electrodes  represent  the  driving  electrodes,  and  the  outer  ones  are  used  as  re¬ 
ceiving  electrodes).  Since  this  potential  is  a  sinusoidal  one,  all  other  potentials 
are  as  well,  and  the  additional  dc  voltages  as  applied  to  rings  11  to  17  are  "riding" 
on  this  sinusoidal  v/ave. 

The  selection  of  the  outer  electrodes  as  receiving  electrodes  was  determined  by 
the  necessity  to  keep  the  input  capacity  of  the  electrometers  connected  to  them 
as  small  as  possible. 

In  the  following  sections,  some  more  theoretical  aspects  of  this  chamber  will 
be  discussed.  In  Part  2  of  this  report,  the  actual  design  of  this  model  will  be 
described  and  the  laboratory  experiments  and  their  results  will  be  reported. 

For  the  actual  measurements  in  the  free  atmosphere,  the  modified  GERDIEN 
chamber  will  descend  from  peak  latitude,  while  a  parachute  or  decelerating  body 
will  keep  the  falling  speed  below  the  velocity  of  sound. 

1.  4.  2  THE  MOBILITY  SPECTRUM  OF  THE  MODIFIED  GERDIEN  CHAMBER 

The  chamber,  as  shown  in  figure  10,  is  designed  as  a  differential  GERDIEN 
chamber  of  the  second  order.  However,  if  all  driving  electrodes  are  connected 
to  each  other,  and  if  all  receiving  electrodes  are  also  connected  to  each  other, 
it  may  be  applied  as  a  differential  chamber  of  the  first  order  (type  b)  as  well. 

The  analytic  expressions  for  the  mobility  spectra  are  not  derived  here.  They  are 
for  the  first  order  mode: 

(2-)2 

dn  T/r  \  co/  dl 

-  «  -  -  -  (76) 

dk  y  j  d  e  d  (U/w) 
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and  for  the  second  order  mode: 

dn  T/t  2 U/eu  {  ,  (77) 

dk  yj  d  eMg 

in  both  cases  for  plate  condensers.  The  same  letters  as  before  are  used.  Both 
modes  are  demonstrated  in  figure  11. 

Again,  it  can  be  seen,  that  for  the  first  order  mode  the  first  derivative  must  be 
used,  while  the  spectrum  is  directly  given  in  the  second  order  mode.  The  cur¬ 
rent  measured  for  the  same  ion  number  density  per  mobility  range  is  the  great¬ 
er  the  more  ions  are  sucked  into  the  chamber  per  second  (Mg),  but  the  parame¬ 
ter  M0  (the  total  amount  of  air  passing  through  the  chamber)  is  no  more  contain¬ 
ed  in  these  equations,  as  it  was  in  the  dc  case  (eqs.  (67)  and  (73)).  That  means 
that  the  air  flow  velocity  in  the  chamber  does  not  influence  the  resolution  any¬ 
more,  This  also  becomes  evident  from  figure  7(C):  the  location  for  the  inter¬ 
ception  of  ions  of  a  certain  mobility  is  determined  by  the  voltage  which  is  suffic¬ 
ient  to  move  the  ion  w*  hin  one  period  of  the  sine  wave  towards  the  receiving 
electrode. 

The  resolution  of  the  measurement  can  be  made  very  good  with  the  first-order 
mode,  but  the  accuracy  is  much  smaller  than  with  the  second-order  application. 
The  resolution  with  the  second-order  mode  depends  on  the  number  of  receiving 
electrodes,  and  on  the  overall  range  of  mobility  to  be  comprehended  in  one  meas¬ 
urement,  since  this  in  turn  determines  the  dU/dx,  or  the  increase  of  driving  volt¬ 
age  along  the  axis  of  the  chamber. 

The  fact  that  the  mobility  spectrum  does  not  depend  on  M0  is  advantageous  for 
measurements  in  higher  altitudes.  For  several  reasons,  it  is  difficult  to  meas¬ 
ure  wind  speed  in  an  atmosphere  of  low  density.  Any  errors  made  in  such  a 
measurement  influence  the  mobility  spectrum  only  linearly,  and  not  in  its  square, 
as  in  the  case  with  the  dc -operated  chamber. 

Resolution  and  accuracy  of  the  measurement  depend  on  many  factors.  The  more 
essential  ones  will  be  discussed  in  the  following  sections,  and  a  particular  dis¬ 
cussion  of  the  limitations  will  be  presented  in  one  of  the  following  chapters. 

1.  4.  3  MEASURING  ION  NUMBERS  BY  MEASURING  CURRENT 

In  the  original  conductivity  apparatus  (RUTHERFORD,  1899;  GERDIEN,  1905; 
and  others)  the  "receiving  electrode"  was  fully  insulated.  Either  the  "charging" 
or  the  "discharging"  methods  have  been  applied.  In  the  first  case,  the  receiv¬ 
ing  electrode  had  no  charge  at  the  beginning,  and  the  rate  of  charging  was  de¬ 
termined  by  measuring  electrometrically  the  potential  difference  between  this 
electrode  and  ground.  After  a  while,  the  electrode  was  artificially  discharged 
and  the  process  began  again.  Correspondingly,  with  the  discharging  method, 
the  receiving  electrode  was  charged  at  the  beginning  of  a  measuring  cycle  and 
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the  velocity  of  the  air  flow  in  the  chamber.  Equations  apply  to  plate  condenser.) 


then  discharged  by  the  arriving  ions  of  the  opposite  sign.  In  these  two  cases 
(which  are  discussed  extensively  in  the  literature)  the  capture  of  the  ion  deter¬ 
mined  the  signal  which  was  read  in  the  electrometer.  Later  on,  the  applica¬ 
tion  of  the  GERDIEN  method  was  frequently  changed  by  the  introduction  of  the 
so-called  current  method.  Here,  the  receiving  electrode  is  connected  to  ground 
by  a  high-ohmic  resistor.  The  current  through  this  resistor  is  measured,  and 
it  is  supposed  that  this  current  represents  the  number  of  impacting  ions  times 
their  charges.  Some  discussion  on  this  principle  and  on  the  problem  of  whether 
or  not  this  supposition  is  correct  has  been  conducted  in  recent  years.  For  the 
integral  GERDIEN  chamber  this  supposition  will  hold;  deviations  will  be  small 
and  hardly  noticeable. 


With  the  modified  GERDIEN  chamber,  the  situation  is  different.  During  a  dis¬ 
cussion  at  a  meeting  in  1963  (DOLEZALEK,  1963),  it  was  stated  that  transients 
caused  at  a  receiving  electrode  by  ions  which  pass  by,  but  do  not  hit  it,  do  not 
disturb  the  measurement  if  only  the  actual  charges  of  impacting  ions  are  meas¬ 
ured  by  applying  an  integration  ever  a  sufficient  time.  A  more  accurate  ex¬ 
pression  was  derived  later  on  by  KADANOFF  (1964),  who  calculated  the  current 
caused  by  the  moving  ion  on  the  basis  tha  cylindrical  symmetry  can  be  supposed. 


The  charges  induced  on  the  inner  and  outer  electrodes  of  a  cyiindrical  conden¬ 
ser  by  an  electric  charge  moving  between  them  are  (not  regarding  edge  effects) 
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where  Q  is  the  moving  charge,  Q1  is  the  induced  charge,  R  and  r  are  the  radii 
of  the  outer  and  inner  electrodes,  res'pectively,  and  pis  the  radius  of  the  loca¬ 
tion  of  Q  at  a  given  time.  The  current  I  flowing  from  the  outer  and  inner  elec¬ 
trode,  respectively,  can  be  calculated  by  differentiation  over  time  (t)  as  follows: 
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With  these  formulas,  and  supposing  certain  mobilities  and  values  for  the  en¬ 
trance  point  of  the  ions  into  the  chamber  (location  of  ion  gate)  as  well  as  radii 
of  the  electrodes,  the  location  of  the  ion  in  the  assumed  chamber,  the  velocity 
of  this  ion  in  the  chamber,  and  the  current  induced  in  the  high-ohmic  resistor 
between  the  outer  electrode  and  ground  have  been  calculated.  They  are  entered 
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in  figures  12,  13,  and  14.  When  interpreting  these  figures,  we  must  keep  in 
mind,  that  the  driving  voltage  has  a  sinusoidal  shape.  This  determines  the  form 
of  these  curves.  A  quality  .ive  interpretation  of  the  figures  and  the  equations  leads 
to  the  following  result. 

A  current  is  caused  by  an  incoming  ion  as  soon,  as  the  respective  receiving  e- 
lectrode  "sees,!  that  ion.  For  the  first  one  of  the  receiving  electrodes,  this  oc¬ 
curs  when  the  ion  enters  through  the  rings  23  and  24  into  the  chamber.  If  this 
ion  is  not  to  hit  the  receiving  electrode,  the  current  caused  by  it  flows  in  one 
direction  as  long  as  the  ion  approaches  the  receiving  electrode,  and  in  the  oppo¬ 
site  direction  when  the  ion  moves  away  from  it.  It  ceases  to  flow,  when  the  e- 
lectrode  cannot  "see"  the  ion  anymore  whether  it  be  that  the  ion  is  intercepted 
at  another  receiving  electrode,  or  that  it  moved  too  far  downstream  the  chamber. 
In  both  cases  the  time  integral  of  the  current  extended  over  the  full  period  is 
zero.  However,  if  the  ion  hits  that  particular  receiving  electrode,  the  meas¬ 
urement  is  done  at  this  very  moment,  the  time  integral  is  not  zero,  but  has  the 
magnitude  of  the  charge  of  that  ion.  In  figures  12,  13,  and  14,  the  parameters 
have  been  chosen  in  such  a  way  that  the  ion  almost  hits  the  receiving  electrode. 

It  can  be  seen,  that  these  three  curves  would  be  very  nearly  the  same  if  the  ion 
would  hit.  That  means  that  during  the  passing  of  an  ion,  an  electrometer  with 
a  very  short  time  constant  attached  to  this  electrode  would  not  be  able  to  dis¬ 
tinguish  between  ions  which  would  hit  and  others  which  would  not.  In  ether  words 
we  again  approach  the  charge  measurement  method  as  used  with  GERDIEN  cham¬ 
bers  in  ea.  iy  days  -  and  it  can  be  the  "charging"  of  the  "discharging"  mode.  This 
depends  on  the  immediate  history  of  the  receiving  electrode  under  consideration. 

If  one  would  allow  ions  of  both  signs  to  be  measured  (i,  e.  ,  open  the  ion  gate 
twice  per  period  of  driving  voltage),  one  would  operate  essentially  with  a  dis¬ 
charging  mode.  This  would  require  a  rather  precise  setting  of  the  integrating 
time  for  each  individual  measurement.  On  the  other  hand,  of  one  opens  the  gate 
only  once  per  period  (compare  figure  9),  only  ions  of  one  sign  enter  and  the 
charging  mode  applies.  In  this  case, 'the  individual  measurement  can  be  extend¬ 
ed  over  several  or  even  many  periods  of  the  driving  voltage,  and  it  becomes  in¬ 
dependent  of  the  wind  velocity  in  the  chamber.  This  is  preferable.  The  dis¬ 
charging  of  the  receiving  electrode  may  then  be  done  by  grounding  it,  or  by 
switching  the  ion  gate  to  let  the  ions  pass  with  the  other  sign.  Again,  the  ground¬ 
ing  would  be  preferable,  since  then  one  always  operates  with  the  charging  mode. 

In  conclusion, the  time  constant  of  the  electrometer  attached  to  the  receiving 
electrode  must  be  great  enough  to  obtain  an  integral  over  a  time  which  is  long¬ 
er  than  the  time  in  which  an  ion  travels  from  the  exit  of  the  ion  gate  to  the  re¬ 
ceiving  electrode  under  consideration.  The  electrometer  time  constant  must 
be  shorter  than  the  period  of  time  between  the  switching  over  from  one  sign  to 
the  other.  If  these  conditions  are  fulfilled,  the  electrometer  should  show  how 
the  capacity  of  the  receiving  electrode  versus  ground  is  charged  up  by  the  in¬ 
tercepted  ions.  Because  of  the  high  ohmic  resistance  across  this  capacity,  this 
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Figure  12  LOCAT ION  OF  MO VI NG  ION  IN  MODI F IED  GERDI E N  CHAMBER 
(Outer  electrode  is  receiving  electrode.) 


charging -up  process  occurs  according  to  an  exponential  function  which,  at  a 
closer  sight,  dissolves  into  a  number  of  stepwise  increases  and  decreases, 
each  pair  of  which  has  the  duration  of  one  period  of  the  driving  voltage.  If  the 
high  ohmic  resistor  is  too  small,  however,  there  may  be  no  increase  at  all.  If 
one  should  be  chosen  too  high,  one  would  risk  getting  an  unkown  leak  resistance 
between  electrode  and  ground. 

1.  4.  4  DETERMINATION  OF  DRIVING  VOLTAGE  PARAMETERS 

The  discussion  of  the  preceding  section  encompasses  essentially  the  general 
theory  for  the  ion  motion  in  the  chamber.  The  formula  for  the  location  of  an 
ion  which  entered  the  chamber  through  the  ion  gate  located  at  the  radius  (r  +  yj  ) 
is: 


2k 
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(cos  (a>t 


6)  >•  cos 


(r  +  yj)- 


(82) 


ar.d,  for  an  ion  entering  when  the  driving  voltage  passes  the  value  zero(i.e.,  6^0) 
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If  p  is  replaced  by  either  R  or  r,  the  equations  provide  the  U/oi  values  for  a 
tangential  impace  on  either  the  outer  or  inner  electrode,  respectively. 


It  is  obvious  from  figure  8  (in  section  1.4.  1)  that  tangential  impact.  Is  the  de¬ 
sirable  method  for  measuring  the  ion  spectrum.  This  fact  is  confirmed  by  the 
quantitative  discussion  of  the  preceding  section. 


However,  it  has  been  shown  by  DOLEZALEK  and  OSTEIl  (1964)  that  even  with 
ions  impinging  at  a  greater  angle,  information  on  the  ion  spectrum  may  be  de¬ 
rived.  This  fact  may  add  valuable  information  in  cases  with  a  swept  or  stepped 
driving  voltage. 

For  the  case  of  tangential  impact  ,  U/oi  values  are  shown  in  figure  15,  covering 
the  mobility  range  which  corresponds  to  altitudes  between  about  30  and  80  km. 

The  four  curves  shown  are  valid  for  (a)  ions  entering  at  zero  driving  voltage 
exactly  (6=0)  and  in  the  middle  between  rings  23  and  24  of  the  laboratory  model 
of  modified  GERDIEN  chamber  (modelA),  and  for  (b)  ions  entering  at  <f>  ~  +  io  degrees, 
again  in  the  middle  between  the  rings,  and,  finally,  (c)  and  (d)  for  ions  entering 
at  6  =  0  but  very  close  to  rings  23  and  24,  respectively.  After  the  selection  of 
the  adequate  U,'«  values,  the  individual  parameters  u  and  <„  must  be  chosen. 

The  amplitude  of  the  driving  voltage,  U  (peak  to  zero)  must  be  greater  than 
disturbing  dc  voltages  which  may  occur  in  the  chamber  (compare  figure  6, 
lower  part),  but  its  size  is  limited  also  by  the  necessity  of  avoiding  some  pos- 
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Figure  15  DRIVING  VOLTAGE  PARAMETERS  UA.F0R  TANGENTIAL  IMPACT  OF  IONS 
(Model  A  of  modified  GERDIEN  chamber:  (a)  ions  entering  in  the  middle  between  rings 
23  and  24  and  at  zero  driving  voltage  (i.e,;  =  0).  (b)  ions  entering  in  the  middle 

between  rings  23  and  24,  at  a  phase  of  10  degrees  prior  or  after  zero  driving 
voltage;  (c)  ions  entering  very  close  to  ring  23  at  zero  driving  voltage;  and 
(d)  ions  entering  very  close  to  ring  24  at  zero  driving  voltage. 
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sible  disturbances  to  be  discussed  in  chapter  1.7  .  Furthermore,  U  must  be 
smali  enough  to  avoid  ionization  by  collision  in  the  chamber,  to  avoid  corona 
discharges  between  the  rings  20  to  27,  and  between  them  and  the  adjacent  rings 
of  the  set  10  to  17.  On  the  other  hand,  the  greater  U,  the  smaller  the  disturb¬ 
ing  influence  of  diffusion. 

The  frequency  of  the  driving  voltage,  «  =  2  u  v  is  determined  by  several 
factors  as  well.  Since  we  want  tangential  impact  at  the  receiving  electrode, 
the  distance  between  the  impact  point  of  the  first  electrode  and  rings  23/24  in 
the  axial  direction  should  be: 

v  2  v  v 

Ax  =  Tv  =  —  =  -  .  (84) 

V  OJ 

Thus,  for  a  given  wind  speed  in  the  chamber  v  and  given  geometrical  dimen¬ 
sions  of  the  chamber,  a  minimum  «  is  determined.  The  same  consideration 
applies  for  the  distances  between  the  first  and  subsequent  receiving  electrodes, 
or  simply  for  the  electrode  length.  This  is  plotted  in  figure  16  for  two  differ¬ 
ent  wind  speeds.  On  the  other  hand,  the  frequency  of  the  driving  voltage  must 
be  kept  in  a  range  in  which  it  is  possible  to  construct  a  properly  working  ion 
gate  impulse  generator  (compare  figure  9  and  section  1.  4.  7  below);  it  should 
be  kept  small  enough  to  allow  a  sufficient  cancellation  of  the  main  displacement 
current  (see  section  1.  4.  5  below).  Furthermore,  the  frequency  of  driving 
voltage  determines  the  opening  period  of  the  ion  gate  (to  a  certain  extent),  and 
since  the  barring  voltage  in  the  gate  must  stay  well  below  the  onset  of  corona 
discharge,  there  is  a  dependence  between  frequency  and  length  of  gate  rings 
(rings  13  and  14  in  model  A). 

1.  4.  5  ION  CURRENT  AND  DISPLACEMENT  CURRENT 

If  the  number  density  of  ions  of  one  mobility  range  is  ,  and  if  all  these  ions 
have  but  one  elementary  charge,  the  electric  current  as  carried  by  these  ions 
to  the  receiving  electrode  is  given  by: 


!k  =  ^Gnke  “  v  +  r)  e 


(85) 


As  pointed  out  in  the  preceding  section,  the  cancellation  of  the  displacement 
current  induced  in  the  receiving  electrode  by  "passing -by"  ions  requires  the 
summing  up  of  the  current  measurement  for  a  time  at  least  as  long  as  that  re¬ 
quired  for  one  ion  to  pass  downstream  the  length  of  one  receiving  electrode.  In 
this  period  of  time,  Ax/v  ,  the  charge  accumulated  on  an  insulated  receiv¬ 
ing  electrode  of  length  Ax  will  be: 
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CIRCULAR  FREQUENCY  OF  DRIVING  VOLTAGE  FOR  A  GIVEN  LENGTH  OF  RECEIVING  ELECTRODE  ACCORDING  TO  u;  =  |- 

Figure  16  RELATION  BETWEEN  FREQUENCY  OF  DRIVING  VOLTAGE  AND 
LENGTH  OF  RECEIVING  ELECTRODES 
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Q2  =  Ax  Ay  2n  (yj  +  r)  nk 


If  the  current  is  collected  over  one  period  of  the  driving  voltage,  one  samples 
all  ions  entered  in  the  time  r  .  This  quantity  may  be  derived  from  figure  9 
and  amounts  to: 

2  v  Ay 

Q3  =  nk  e  v  [Ayz  +  2 Ay  (r  +  yj)]  rv  +  — - —  -  / 

L  k  EG 

where  the  letters  d  and  r  of  figure  9  have  been  replaced  by  the  symbols  Ay 
and  (r  +  yj  )  .  Here,  EG  is  the  voltage  applied  to  the  ion  gate  for  closing  it, 
and  l  is  the  length  of  the  ion  gate  (ring  14  or  13). 


If  the  receiving  electrode  under  consideration  has,  versus  ground,  a  capacity  of 
Ck  and  if  the  leak  resistance  is  sufficiently  great  so  that  the  product  of  Ck 
with  it  is  much  greater  than  the  duration  of  one  period  of  the  driving  voltage, 
this  capacity  will  get  a  voltage  of: 


U3  = 


—  [Ay2  +  2  Ay  (r  +  y^] 


from  the  ions  of  one  mobility  range  entering  during  one  period  of  driving  vol¬ 
tage. 


In  addition  to  the  displacement  current  induced  by  "passing-bv"  ions,  we  also 
have  the  displacement  current  induced  directly  by  the  driving  voltage.  Since 
the  driving  electrodes  are  supplied  by  an  ac  voltage,  we  will  get  another  ac 
voltage  with  a  90  degree  phase  shift  on  the  receiving  electrodes  as  well,  oven 
if  there  were  no  ions  in  the  chamber  at  all.  If  we  were  to  measure  the  current 
per  receiving  electrode  exactly  for  one  period  or  for  a  full  number  of  periods, 
this  "main  displacement  current"  would  just  cancel.  However,  there  are  two 
reasons  why  we  should  not  count  on  this  automatic  cancellation:  at  first,  it 
would  be  very  difficult  to  arrange  the  measurement  to  a  full  number  of  periods 
exactly,  and  since  this  main  displacement  current  will  be  much  greater  than 
the  current  provided  by  the  ions,  any  small  error  in  the  timing  would  cause 
great  disturbance.  Second,  the  main  displacement  is  assumed  to  be  so  great 
that  it  would  quickly  saturate  any  sensitive  electrometric  amplifier.  It  is 
known  that  the  desaturation  of  these  amplifiers  occurs  with  an  extremely  long 
time  constant.  Thus,  the  amplifiers  would  just  jump  from  positive  to  negative 
saturation  and  vice  versa,  and  no  measurement  would  be  possible. 
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For  these  reasons,  a  more  detailed  discussion  on  the  main  displacement  cur¬ 
rent  is  necessary. 

To  obtain  some  numerical  values,  we  take  the  more  disadvantageous  case,  -- 
that  the  admittance  between  receiving  electrode  and  ground: 


is,  due  to  a  relatively  high  frequency,  determined  by  the  capacitance.  For 
i'  =  100  Hz  (i.  e.  ,  co  =  628  sec  _1  ),  this  is  the  case,  if  the  numerical 

value  of  Rk  in  Ohm  is  more  than  1.  6  x  10“^/  Ck>  where  C  k  is  measured  in 
Farads;  e.g.  ,  if  the  capacity  is  10  pF,  the  resistance  must  be  greater  than  1.  6 
Cigaohm  to  have  a  contribution  of  less  than  1  per  cent  to  the  admittance.  For 
v  -  1  kHz,  the  resistance  must  only  be  greater  than  160  megohms.  In  this 
case,  the  displacement  current  may  be  calculated  by  the  series  connection  of  Ck  an<^ 
Cpj,  where  Cp  is  the  effective  capacity  between  one  receiving  electrode  and  the 
driving  electrodes: 
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this  is  the  amplitude  (peak  to  zero)  of  the  main  displacement  current. 


(90) 


For  example,  if  the  modified  GERDIEN  chamber  is  moving  through  the  atmos¬ 
phere  at,  say,  an  approximate  altitude  of  50  km  with  half  the  velocity  of  sound, 
and  the  number  density  of  all  ions  is  5  x  10^  m“^,  and  we  want  to  measure  in  a 

4  Q  *3 

mobility  range  where  5x10°  m"°  ions  exist,  we  derive,  from  eq.  (85),  an  ion 
current  to  the  adequate  receiving  electrode  of  about  (using  the  dimensions  of 
figure  10): 


1.13  x  HT11  A  . 

3  —  | 

On  tlis  other  hctndj  with  to  =  6.  28  x  10°  sec."  and  U  =40V  (using  the  U/  co  value 
from  figure  15),  and  with  Ck  =  50  pF  and  CD  =  15  pF  ,  we  get  for  the  peak 
value  of  the  main  displacement  current  to  the  same  receiving  electrode: 

2.88  x  10-6  A. 

The  intended  accumumulation  of  charge  would  reduce  the  difference,  and  if  one 
measurement  is  made  by  integrating  over  a  great  number  of  periods,  the  signal 
can  be  increased  by  several  orders  of  10,  while  the  noise  due  to  displacement 
current  never  will  become  greater  than  the  charge  accumulation  during  one  half 
period. 
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However,  for  a  safe  measurement,  the  displacement  current  must  be  elimina¬ 
ted.  Figure  17  shows  a  workable  principle.  It  requires  a  duplication  of  the 
driving  voltage  amplitude.  The  critical  point  of  the  arrangement  of  figure  17 
in  the  potentiometer  Pq,  which  must  allow  a  very  sensitive  regulation.  In  fact 
this  potentiometer  was  made  up  by  an  arrangement  of  eleven  sensitive  single 
potentiometers.  This  done,  the  displacement  current  could  be  cancelled  with  ?. 
remaining  part  of  slightly  less  than  lO-^A.  An  even  better  cancellation  would 
probably  be  possible  by  a  better  shielding  and  a  more  sophisticated  application 
of  the  condenser  of  about  8  pF.  The  whole  arrangement  of  figure  17  must  be 
carefully  kept  free  of  phase  shifts. 

1.  4.  6  THE  AMPLIFICATION  OF  THE  SIGNAL 

According  to  eq.  (85)  the  current  to  be  measured,  will  be  smaller  than  10“^  A, 
and  it  should  be  tried  to  amplify  a  current  of  about  5  x  10*^  A  in  order  to  pro¬ 
vide  ample  flexibility.  This  current  will  be  ac  and  it  would  be  desirable,  if  a 
frequency  of  several  kHz  still  could  be  handled.  A  discussion  of  the  problems 
involved  has  been  presented  by  DOLFZALEK  and  OSTER  (1964)  and  will  not  be 
repeated  here. 

A  promising  beginning  to  solving  this  problem  is  given  by  the  so-called"  ad¬ 
mittance  neutralizer"  (VOLKERS,  1963/4).  Here,  an  input  stage  with  a  rel¬ 
atively  great  impedance  (better  than  10®  H )  is  provided.  In  addition,  artificial 
"negative  capacitances"  and  "negative  conductances"  are  fed  into  the  input,  with 
the  goal  of  raising  the  input  resistance  to  infinity  and  to  lower  the  input  capaci¬ 
tance  to  zero.  In  principle,  this  goal  could  be  reached,  and  we  succeeded  in 
measuring  alternating  currents  of  about  3  x  10"^  A  at  2  kHz.  It  turned  out, 
however,  that  a  further  decrease  was  impossible  because  of  the  existing  noise 
and  it  was  already  then  difficult  to  get  the  proper  settings  of  the  regulation  po¬ 
tentiometers.  It  seems  promising,  to  increase  the  "natural"  input  impedance 
to  about  1014  by  applying  a  good  electrometer  tube  or  a  modern  field  effect  tran¬ 
sistor  in  the  input  stage  of  the  admittance  neutralizer.  In  addition,  a  filter  for 
frequencies  greater  than  about  10  kHz  may  decrease  much  of  the  thermal  noise. 
It  has  not  as  yet  been  possible  to  test  this  possibility  by  an  experiment. 

Another  possibility  is  given  by  the  combination  of  an  electrometric  input  stage 
and  an  operational  amplifier.  A  practical  arrangement,  following  a  proposal  of 
SAPPUPO  (1964)  is  shown  in  a  later  chapter  of  this  report.  Here,  the  critical 
point  is  the  proper  balancing  of  the  zero  signal  put  into  one  of  the  two  inputs  of 
the  operational  amplifier  by  an  artificial  voltage  put  into  the  second  one.  Ex¬ 
perimentally,  we  had  good  results  down  to  current  of  10-1®  A  and  even  less, 
with  frequencies  up  to  several  hundred  Hertz.  This  arrangement  could  be  made 
more  sophisticated  by  applying  the  bridge  electrometer  of  DOLEZALEK  (1956, 
1961)  and  feeding  the  outputs  of  the  two  members  of  the  bridge  into  two  opera¬ 
tional  amplifiers,  which  then  could  be  combined  by  a  third  one.  The  main  ad¬ 
vantage  of  the  bridge  principle  is  the  provision  for  the  cancellation  of  disturb¬ 
ing  elements;  in  the  case  of  the  modified  GERDIEN  chamber,  it  could  be  for  a 
cancellation  of  the  main  displacement  current. 
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Figure  17  PRINCIPLE  FOR  THE  CANCELLATION  OF  THE  MAIN  DISPLACEMENT  CURRENT 

IN  THE  MODIFIED  GERDIEN  CHAMBER 
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1.  4.  7  ION  FENCE  AND  ION  GATE 


As  discussed  in  section  1-.  4.  1,  air  is  flowing  through  the  full  cross  section  of 
the  modified  GERDIEN  chamber,  but  ions  are  entering  only  through  a  circular 
slit.  From  the  rest  of  the  air,  the  ions  are  removed  by  the  application  of  strong 
dc  fields. 

We  call  this  the  ion  filter  or  the  ion  fence,  and  in  Model  A  (Laboratory  Model) 
it  was  made  up  by  rings  10  to  13  and  14  to  17  (figures  9  and  10).  The  field  to  be 
applied  must  be  stronger  than: 

.  (91) 

vd 


where  v  is  the  wind  velocity  through  these  rings,  d  is  the  distance  between  two 
adjacent  rings,  k  is  the  smallest  expected  mobility,  and  l  is  the  axial  length  of 
the  rings.  The  relation  between  uG  =  Edand  length  /  is  plotted  in  figure  18.  Here, 
the  mobility  values  are  converted  into  altitude  values  according  to  table  202  in 
serHon  1.  6.  4,  and  velocity  values  have  been  taken  from  PEDERSEN  (1964).  It 
should  be  mentioned,  however,  that  the  d-values  applied  in  figure  18,  being  1 
cm,  may  well  be  lessened  for  the  measurement. 

On  the  other  hand,  the  a.pplied  voltage  must  be  so  small  that  no  corona  discharge 
occurs.  Because  of  the  considerable  length  of  the  cathode  fall  under  low  pres¬ 
sures,  the  PASCHEN  curve  will  have  a  minimum  at  rather  high  voltages  (in  the 
order  of  a  few  hundred  Volt).  This  has  been  proven  experimentally  in  a  bell 
jar  (see  figure  43). 

In  Model  A,  the  "ion  gate"  is  situated  between  rings  13  and  14f  Again,  the  bar¬ 
ring  voltage  during  the  time  in  which  the  gate  has  to  be  closed,  must  suffice  eq. 

(91) ,  while  the  shortest  possible  opening  period  is  given  by: 

l  Ay  (92) 

~  ~v  Ti 

The  number  of  ions  entering  into  the  chamber  during  this  time  has  already 
been  given  by  eq.  (88)  in  section  1.  4.  5.  All  of  these  equations --(88),  (91),  and 

(92)  --are  easily  derived  from  figure  9.  The  letter  d  of  figure  9  and  eq.  (91)  is 
replaced  by  Ay  in  eq,  (92),  since  Ay  was  used  by  us. for  the  denomination  of  the 
slit-width  of  the  ion  orifice. 

The  practical  verification  of  a  generator  to  supply  the  pulses  for  ion  gate  (trig¬ 
gered  by  the  driving  voltage),  riding  on  that  portion  of  driving  voltage  which  is 
adequate  to  rings  23  and  24,  poses  many  difficulties.  In  fact,  this  generator 
and  the  amplification  of  the  signal  have  made  necessary  many  experiments  and 
required  much  effort. 
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1.  4.  8  THE  LENGTH  OF  THE  RECEIVING  ELECTRODES 


In  the  preceding  section  eq.  (92)  also  determines  practically  the  length  Ax  of  the 
receiving  electrodes.  This  can  be  seen  from  the  following  consideration: 

The  opening  period  for  the  ion  gate  r  ,  should  only  be  a  small  portion  of  the 
period  T  of  the  driving  voltage.  Otherwise,  the  resolution  of  the  measurement 
would  go  down  and  the  ion  number  densities  in  the  lower  mobility  ranges  would 
be  falsified  (see  section  1.  4.  1).  To  keep  r  small  enough,  the  wind  velocity  v 
must  be  great  enough. 

On  the  other  hand,  we  want  the  ions  to  be  moved  by  a  full  period  while  passing 
one  receiving  electrode, 

2  7rv  (93) 

Ax  >  - 

(o 


(compare  figure  16  and  eq.  (84)).  If  we  assume  a  value  T/r  of  10,  the  com¬ 
bination  of  eqs.  (92)  and  (93)  gives  approximately: 


2  77  V 

Ax  >  -  >  10  l 

co 


(94) 


here  the  second  member  on  the  right  hand  side  of  eq,  (92)  has  been  neglected. 

Eq.  (94)  is  essential  for  the  determination  of  the  overall  length  of  the  modified 
GERDIEN  chamber.  Lastly,  the  limits  are  set  by  the  smallest  practical  length 
of  ion  gate  l  and  the  greatest  allowable  barring  field  in  the  gate  (see  eq.  (91)). 
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1.  5  AIR  FLOW  IN  THE  GERDIEN  CHAMBER 


Two  characteristics  of  the  air  flow  inside  the  GERDIEN  chamber  may  influence 
the  measurement  of  mobility  and  conductivity:  turbulence  and  rounded  wind  pro¬ 
files.  For  turbulence,  refer  to  section  1.  3.4.  The  fact,  that  the  air  inside  the 
GERDIEN  chamber  will  not  have  the  same  velocity  at  all  places  in  the  cross  sec¬ 
tion  will  be  discussed  in  this  chapter. 

In  a  GERDIEN  chamber  operated  by  a  dc  driving  voltage,  the  accurate  knowledge 
of  the  wind  velocities  for  all  radii  which  are  crossed  by  the  moving  ion  is  essen¬ 
tial.  In  the  "modified  GERDIEN  chamber,  "  operated  by  ac  driving  voltage,  it 
is  sufficient,  if  a  consideration  of  wind  speeds  leads  to  the  conclusion,  that 
they  do  not  interfere  with  the  measurement.  This  is  the  case  if  the  ion  during 
its  vertical  motion  is  not  carried  past  the  adequate  receiving  electrode  by  too 
great  a  wind  velocity  for  some  portion  of  its  path. 

M.  BERNSTEIN  (1965)  has  conducted  theoretical  investigations  on  the  wind  vel¬ 
ocity  profile  for  the  modified  GERDIEN  chamber,  Model  Bl.  This  model  will 
be  used,  at  first,  for  some  dc  measurements.  It  is  described  in  Part  3  of  this 
report.  The  following  presentation  reflects  the  result  of  BERNSTEIN'S  investi¬ 
gations. 

Neglecting  the  influence  of  the  "spider  legs"  (i.  e. ,  the  four  legs  at  the  upstream 
end  of  chamber,  supporting  the  sets  of  rings),  the  velocity  distribution  is  given 
by  the  equation 


(95) 

v  (r)  =  V  [B  —  k(r)2  +  <f>  log  r] 

where  the  significance  of  the  letters  is  different  from  that  given  in  the  list  of  symbol 

v(r)  is  the  limiting  flow  velocity  in  the  region  from 
r  to  r  +  dr;* 

r  is  the  distance  from  the  axis  of  the  chamber; 

v  is  the  mean  flow  velocity  in  the  region; 
a2  In  a  —  (a2  —  1)  In  a  Rj 

3=  2  - 

(cr2  +  1)  In  a  —  (a  +  1) 


k  =  2 

5  =  2 


In  a 

Rj2  (a2  +  1)  In  a  -  (a2  +  1* 


(a2  4  1)  In  a  -  (a2  +  1) 


‘"Limiting’’  flow  velocity  refers  to  the  "final”  velocity  distribution  which  is  asymptotically  approached  when  x  »  d 
(refer  to  cq.  (96)). 
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a  -  r2/Ri  ,  the  ratio  of  the  outer  to  the  inner 

radius  of  a  "region,  11  and  a  "region" 
is  a  cylindrical  layer  within  the  cy¬ 
lindrical  condenser,  determined,  in 
a  particular  case,  by  two  adjacent 
coaxial  rings  at  the  upstream  end  of 
the  chamber. 


Aerodynamic  assumptions  behind  this  expression,  eq,  (95),  are  that  the  flow  is 
steady,  viscous,  and  incompressible,  and  that  the  entering  wind  is  at  a  uniform 
velocity  throughout  a  region. 

Comparing  the  flow  in  the  "regions"  as  defined  above  with  the  simpler  case  of 
flow  between  parallel  planes  (making  the  approximation  that  the  flow  in  the  in¬ 
dividual  regions  of  the  GERDIEN  chamber  may  be  considered  as  a  flow  between 
parallel  planes  (R2  -  Rj  «Rj)),  limiting  velocities  are  shown  in  figure  19  for  sev¬ 
en  different  regions,  defined  in  caption  to  figure  19. 


From  figure  19  it  is  seen  that  the  flow  in  the  parallel  plane  case  is  a  good  ap¬ 
proximation  of  the  flow  in  the  coaxial-cylinder-enclosed  regions  of  the  chamber, 
using  the  limiting  velocity  distribution  as  a  criterion.  Thus,  an  analysis  of  the 
flow  in  the  parallel-plane  case  will  yield  useful  information. 


To  determine  the  shape  of  the  flow  velocity  distribution  in  the  parallel-plane 
case,  a  parameter  $  is  introduced  by 


4> 


4  x  v 

77 


(96) 


where 


x  is  the  distance  parallel  to  the  velocity 
vectors,  counted  from  the  upstream  end 
of  the  chamber 

v  is  the  kinematic  viscosity  of  the  air, 
d  is  the  separation  between  the  parallel  planes. 

4 

In  figure  20  velocity  distributions  corresponding  to  several  selected  values  of 
d>  are  given.  For  $  >0.16  ,  the  flow  profile  becomes  indistinguishable  from  the 
limiting  velocity  distribution. 

For  region  6,  however,  information  on  <£is  not  sufficient  to  determine  the  actual 
velocity  distribution;  four  other  factors  must  also  be  considered.  They  are: 
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between  parallel  planes. 

is  the  normalized  radial  separation  between  the  region's  boundaries; 


* 
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Figure  20  THE  RATIOv/v  FOR  DIFFERENT  PARAMETERS  <1> 


(1)  that  there  is  actually  a  difference  between  the  flow  velocity  distribution  in 
region  6  and  the  studied  parallel-plane  case  flow  velocity  distribution;  (2)  that 
the  velocity  profile  is  decidedly  non-uniform  at  the  entrance  of  region  6;  (3)  that 
there  are  turbulent  values  from  the  rings  and  their  supports  ("spider  legs");  and 
(4)  that  the  flow  may  not  be  laminar  in  region  6.  For  a  quantitative  evaluation 
of  these  four  factors,  knowledge  of  the  actual  mean  wind  velocity  v  in  the  cham¬ 
ber  is  necessary.  This  wind  velocity  may  be  found  upon  considering  the  require¬ 
ment  that  laminar  flow  must  be  maintained  in  region  2  (ion  gate)  for  obvious 
reasons,  together  with  the  desire  to  increase  the  mean  wind  velocity  as  much  as 
possible.  As  it  is  well  known,  undamped  turbulence  does  not  occur,  if  the 
REYNOLD'S  number 


v  d 

Re  =  -  <  1400 

v 


(97) 


For  region  2,  we  may  derive,  by  using  v2  =  1400  —p 

and  V  =  v  (A  /A  )  =  1.41  v 
2  oo  oo  ent  00 

where  vM  is  the  falling  velocity  of  the  chamber 

Am  is  the  area  of  entering  flow  far  ahead  of  chamber 
Aentis  the  available  area  inside  the  chamber 


that  the  maximum  allowable  falling  velocity  is 


(Vlmax  =  1-00  x  105  x  v  . 

Kinematic  viscosity  is  given  for  the  altitude  h  in  km  by 
v  =  iq(0-0694  h  —  5.17)  m2/sec  _  , 

Thus  we  arrive  at  a  value  of 


(98) 

(99) 


(100) 


(101) 


Voo  =  10(0.0336  h  +  0.211)  m/sec  *  Ouc, 

From  this  value  for  the  parameter  v  and  the  parameter  <£can  be  calculated, 

for  each  region,  and  for  any  altitude  between  10  and  40  km.  This  has  been  done, 
and  the  result  is  communicated  in  the  following  table; 


*  It  should  be  remembered,  however,  that  laminar  flow  may  still  exist,  if  Re  >  1400.  This  should  be  tested  experimentally 
in  the  wind  tunnel.  For  our  theoretical  consideration,  we  stay  at  values  of  Re  <  1400. 
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PARAMETER  FOR  REGIONS  J  TO  5  IN  DIFFERENT  HEIGHTS 


(Calculated  for  the  downstream  end  of  each  region  between  adja¬ 
cent  rings  (exit  from  rings  into  chamber  proper))' 


h  in  km 

Parameter  cSfor  Region  No. 

'Values  to 
be  multiplied 
with 

1 

2 

3 

4 

5 

10 

4.  7 

6.  0 

6.  3 

6.  0 

3.  8 

x  10~2 

15 

6.  1 

8.  6 

9.  0 

8.  6 

3.4 

x  10“2 

20 

9.9 

13.  0 

13.  0 

13.  0 

7.9 

x  10"2 

25 

1.  5 

1.6 

1.  6 

1.  6 

1.  2 

x  10”1 

30 

1.  6 

1.  6 

1.  6 

1.6 

1.2 

x  10"1 

35 

1.  6 

1.  6 

1.  6 

1.  6 

1.  2 

x  10"  1 

40 

1.  6 

1.  6 

1.  6 

1.  6 

1.2 

1 

o 

r— < 

X 

At  the  exits  from  the  rings  into  the  chamber  proper,  the  velocity  distributions 
in'region  1  to  5  begin  to  influence  the  ion  trajectories  in  region  6.  Figure  21  is 
a  plot  of  the  parameter  ^  for  different  altitudes;  here  x  is  the  distance  measured 
from  the  end  of  regions  1  to  5  or  the  beginning  of  region  6. 

From  the  knowledge  of  parameter  <£,  v.elocity  distributions  in  regions  1  to  5  can 
be  derived,  and  a  first  approximation  to  the  velocity  distribution  in  region  6  can 
be  found,  applying  the  information  from  figures  19  and  20.  The  average  velo¬ 
city  in  region  6  is  given  by 

\6)  =  U5  v~  • 

For  a  better  approximation  to  the  velocity  distribution  in  region  6,  the  other 
three  factors,  as  listed  above,  may  be  studied.  It  can  be  shown  that  the  non¬ 
uniformity  of  the  entrance  profile  at  the  upstream  end  of  region  6  is  absorbed  in 
about  0.  1  meter,  and  that  the  difference  between  flow  velocity  inside  the  wake 
and  outside  of  it  is  not  more  than  0.  5  m/sec.  Furthermore  it  can  be  shown,  by 
application  of  the  REYNOLD'S  number  criterion,  that  a  fully  developed  turbu¬ 
lence  in  region  6  cannot  occur  if  h  >  26  km  .  Since  no  perturbation  mechanism 
seems  to  be  present  in  region  6,  it  is  entirely  possible  that  the  flow  will  never 
become  a  fully  turbulent  flow.  However,  even  if  it  does,  velocity  distributions 
could  be  calculated. 
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For  the  genuine  modified  GERDIEN  chamber  with  ac  operation,  this  discussion 
on  the  air  flow  inside  the  chamber  for  a  given  falling  speed  is  sufficient.  How¬ 
ever,  if  the  chamber  will  be  used  in  a  dc  operation,  a  more  detailed  discussion 
for  region  6  may  be  required.  This  discussion  can  at  best  be  conducted,  if  the 
motion  of  the  ions  under  the  influence  of  the  electric  field  is  considered.  This 
would  then  be  a  special  case  for  the  dc  operation,  and  will  not  be  reported  here. 


1.  6  EXPECTED  ENVIRONMENTAL  CONDITIONS 
FOR  ION  MEASUREMENTS  IN  THE  TERRESTRIAL  STRATOSPHERE 


AND  MESOSPHERE 


1.6.1  INTRODUCTION 

It  has  often  been  said,  that  the  range  between  30  and  80  km  is  the  least  known 
one  in  the  terrestrial  atmosphere.  Nevertheless,  for  the  design  of  any  instru¬ 
ment  which  is  to  make  measurements  in  this  range,  some  theoretical  or  obser¬ 
ved  values  must  be  assumed.  To  create  a  basis  for  out  discussion,  we  have 
compiled  such  values  from  other  publications.  Since  this  was  done  a  few  years 
ago,  some  more  reliable  or  accurate  values  have  been  published.  In  this  re¬ 
spect,  we  refer,  above  all,  to  the  paper  of  COLE  and  PIERCE  (1965).  While 
the  more  accurate  values  are  presently  being  used  for  the  improvement  of  our 
design  and  will  be  used  later  on,  the  earlier  tables  still  are  valuable  for  any 
first  approximation.  For  this  reason,  they  are  printed  unchanged  in  the  follow¬ 
ing  sections  1.  6.  3  and  1.  6.  4. 

1.  6.  2  REMARK  ON  THE  MEAN  FREE  PATH 

Table  104  in  section  1.  6.  3  communicates  two  sets  of  values  for  the  mean  free 
path  of  neutral  particles.  For  some  of  our  later  discussions,  approximate  da¬ 
ta  on  the  actual  free  paths--in  particular,  data  on  the  occurrences  of  free  paths 
longer  than  th'e  mean  one--is  important.  Therefore,  in  figures  22  and  23,  the 
percentage  of  free  paths  longer  than  the  mean  ones  are  presented,  calculated  by 
use  of  MAXWELL'S  velocity  distribution  function,  for  the  two  sets  of  values  in 
table  104. 

Of  course,  the  free  path  for  ions  is  different  from  the  path  for  neutral  particles. 
For  a  recent  discussion  refer  to  SIKSNA  1964. 
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Figi™  22  PERCENTAGE  OF  NEUTRAL  PARTICLES  IN  THE  ATMOSPH ERE , 
WHOSE  FREE  PATH  IS  LONGER  THAN  l 
(The  results  ore  ootained  by  using  the  mean  free  path  values  from  the 
first  column  of  tcbltr  104.) 
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iguro  23  PERCENTAGE  OF  NEUTRAL  PARTICLES  IN  THE  ATMOSPHERE, 
WHOSE  FREE  PATH  IS  LONGER  THAN  l 

(The  results  oro  obtained  by  using  the  mean  (ree  path  values  from  the 
second  column  of  table  104.) 


1.6.3  TABLES  FOR  NON-ELECTRIC  PARAMETERS 


The  T- values  are  taken  fro. n:  The  f -values  are  take  from: 

?.s.  JOKESO::,  Satellite  Environment  Handbook  F.S.  JOiL'SO!’,  S-tellit.  Enviro.’icent  Jancibo'n 

Stanford  University  Press  Stan  fora  Uuvcrmty  Press 

Stanford,  California,  1961  Stanford.  California,  . -'e,i 
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.RELATIVE  TEMPERATU 


S."iATIVE  PRESSURE 


h 

[km] 
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2.96x10° 

4.875X10'1 
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l.lOxlO3 

9 .333x10*  14 
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3.36/.103 

-4 

3.696x10 

70 

1* 

1.35x10 

-4 

1.475x10 

80 

0.9  xio'1 

-5 

2.317x10 

These  values  are  calculated  from 

the  values  of  tables  101  and  103 


1 


1.225 

.  -1 
6.091x10 

4.41*3x10  1 

2. 199x10* 1 

1. 270x10' l 

,  -1 
1.216x10 

1.651X10"1 

9.853xlo"2 


1.6.4  TABLES  FORELECTRIC  AND  DIFFUSION  PARAMETERS 


WHIPPLE,  and.  CLARK  (Journ.  Geoph.  Res.  64,  1959,  1363-1370) 
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1.7  LIMITATIONS  FOR  THE  MEASUREMENT  OF  IONIC  MOBILITY 
IN  THE  UPPER  ATMOSPHERE 


1.7.1  INTRODUCTION 

When  discussing  the  characteristic*  of  atmospheric  ions  and  the  definition  of 
ionic  mobility  in  the  earlier  chapters,  we  have  also  laid  the  foundation  for  a  dis¬ 
cussion  of  some  of  the  facts  which  determine  the  limits  of  the  applicability  of 
these  definitions.  When  ascending  in  the  atmosphere,  the  free  paths  of  the  at¬ 
mospheric  constituents  increase,  the  number  of  collisions  between  the  particles 
decreases,  and  parameters  which  involve  a  great  number  of  collisions  become 
more  and  more  inapplicable. 

This  gradual  changeover  is  of  great  theoretical  and  practical  interest.  The  fo¬ 
cal  point  is  given  by  the  problem  of  ontology- -"how  far  up"  we  still  can  use  the 
denomination  mobility.  At  an  altitude  of  about  80  km,  the  mean  free  paths  be¬ 
come  so  long,  that  it  becomes  unfeasible  to  construct  instruments  which  could 
measure  mobility.  Thus,  we  cannot  measure  it  any  more;  but  for  a  larger  scale 
treatment  of  atmospheric  processes;  the  parameter  mobility  may  still  be  appli¬ 
cable.  Still  higher  up,  this  possibility  also  ceases.  In  this  transition  region, 
the  application  of  denomination,  such  as  mobility  and  conductivity,  must  be  done 
with  caution,  and  it  is  an  open  question  as  to  how  far  out  this  may  still  be  toler¬ 
ated. 

The  different  aspects  of  mobility  and  the  different  problems  of  measuring  it  get 
into  these  difficulties  at  different  heights.  For  this  reason,  we  have  to  discuss 
these  changeovers  if  we  want  to  establish  safe  limits  for  the  application  of  our 
measuring  method.  On  the  other  hand,  some  of  the  difficulties  of  the  measure¬ 
ment  of  ionic  mobilities  and  conductivities  existing  in  the  lower  atmosphere  (see 
section  1.  3.4)  become  less  difficult  in  the  uppef1  atmosphere.  This  is  so,  for 
example,  for  the  problem  of  laminary  and  turbulent  flow  inside  a  GERDIEN 
chamber.  In  a  more  rarefied  atmosphere,  the  REYNOLD'S  numbers  are  such 
that  laminary  flow  is  always  ascertained.  Other  difficulties  change  their  nature, 
such  as  the  problems  provided  by  the  diffusion. 

1.7.2  DIFFUSION 

In  the  modified  GERDIEN  chamber,  diffusion  effects  provide  disturbances  at  two 
locations:  in  the  "ion  gate"  and  in  the  chamber  itself.  In  the  ion  gate,  no  electric 
field  exists  while  the  ions  pass  through.  There  will  be  a  loss  of  ions  to  the  rather 
closely  situated  walls.  Inside  the  chamber,  diffusion  effects  will  tend  to  broaden 
the  ion  path  and  thus  create  a  loss  in  resolution,  and  some  ion  loss  to  the  elec¬ 
trodes  where  the  ion  path  approaches  them  closely. 

Three  different  types  of  diffusion  have  to  be  considered: 
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1.  Gas  kinetic  diffusion,  governed  by  the  "diffusion  coefficient,  D,  and 
the  gradient  of  ion  density,  grad  n 

2.  "Electric  diffusion"  caused  by  the  mirror-charge  effect  attracting  ions 
to  walls  and  electrodes 

3.  "Eddy  diffusion"  due  to  turbulent  flow  of  the  air  containing  ions  by  which 
a  much  greater  number  of  ions  is  put  into  locations,  where  great  density 
gradients  cause  them  to  follow  diffusion. 

The  diffusion  coefficient  is  related  to  the  mobility  of  the  ions  according  to 

EINSTEIN's  law,  which  we  write  in  the  form 

\ 

D  .  k  -HI  <103> 

e 


where 

D=  diffusion  coefficient  of  FICK's  Law  in  sec~* 
k  =  ion  mobility  in  m2  V sec-1 

k  =  BOLTZMANN  constant,  1.4x  10“^  Wsec  degree"^ 

T=  temperature  in  degrees  KELVIN 

e  =  elementary  charge,  1.  6  x  10  Coulomb. 

Values  for  D,  calculated  according  to  eq.  (103),  are  listed  in  table  208  in  sec¬ 
tion  1.  6.  4. 

In  the  same  section,  in  tables  209  and. 210,  some  other  parameters  are  deiived 
from  D,  namely,  flux  numbers  and  relative  flux  numbers.  These  relative  flux 
numbers  give  the  ratio  of  ions  removed  due  to  diffusion  out  of  an  ion  stream  un¬ 
der  certain  ion  number  density  gradients.  In  fact,  table  210  gives  a  removal 
velocity,  and  this  velocity  can  be  compared  with  the  wind  velocity  in  the  modi¬ 
fied  GERDIEN  chamber  for  a  first  approximation  to  the  measuring  errors  caused 
by  diffusion.  This  statement  is  based  on  the  fact  that  the  walls  represent  ab¬ 
solute  ion  sinks,  all  ions  touching  the  wall  are  lost,  they  deposit  their  electric 
charge  and  turn  to  neutral  molecules.  In  fact,  it  is  not  even  necessary  for  them 
to  touch  the  wall,  since  mirror-charge  become  very  great  when  the  ions  ap¬ 
proach  closely;  they  are  practically  lost  as  soon  as  they  achieve  a  certain  ap¬ 
proach.  While,  in  principle,  this  method  for  an  estimation  of  ion  losses  due  to 
diffusion  (to  walls  as  well  as  the  broadening  of  an  ion  trail  in  an  otherwise  neu¬ 
tral  gaseous  environment)  is  certainly  a  correct*  one,  it  is  difficult  to  apply  it 

*  Within  the  vaguely  defined  limits  of  FICK's  diffusion  laws. 
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for  a  more  accurate  consideration,  since  the  definition  of  ion  density  gradients 
becomes  difficult. 

Before  dealing  with  another  approach,  we  want  to  use  the  numbers  of  table  210 
of  section  1.  6.  4. ,  and  other  similarly  derived  numbers  for  a  first  approxima¬ 
tion. 

We  see  from  table  210  that  in  an  altitude  of  80  km  the  ion  removal  velocity  for 
ions  in  the  first  millimeter  is  about  half  the  speed  of  sound.  However,  at  the 
other  end  of  the  table,  the  removal  speeds  are  in  the  order  of  very  little  frac¬ 
tions  of  1  m/sec,  and  thus  negligible.  Disturbances,  according  to  this  con¬ 
sideration,  may  be  expected  in  altitudes  greater  than  60  to  70  km,  and  it  is  here 
where  loss  factors  should  be  determined  experimentally. 

For  another  approximation,  based  on  the  same  concept,  we  use  the  values  of 
table  204  (section  1.6.4.),  dealing  with  average  distances  between  ions.  From 
them,  we  determine  the  thickness  of  the  "first  layer"  of  ions  to  —  yT  of  this  dis¬ 
tance.  Thus,  v  e  obtain  ion  density  gradients  for  this  first  layer  in  the  order  of 
1011  to  1014  m-4.  The  great  discrepancy  stems  from  the  fact  that  the  average 
distance  is  calculated  according  to  two  different  assumptions,  both  hypothetical 
ones,  as  explained  in  section  1.6.4.  The  relative  ion  flux,  calculated  for  the 
first  layer,  then  comes  to  values  between  600  and  2000  m/sec  for  the  height  of 
80  km.  That  means  that  here  the  first  layer  disappears  completely  (for  a  wind 
velocity  of  300  m / sec,  we  get  a  traveling  distance  paralled  to  the  wall  of  only 
0.3  m).  For  the  "second  layer,  "  the  situation  is  somewhat  different.  Since 
these  ions  are  not  directly  confronted  with  an  absolute  ion  sink,  the  ion  density 
gradient  is  smaller. 


A  more  sophisticated  approach  may  be  taken  by  starting  from  equations  quoted 
by  L.  B,  LOEB  (1955,  p.  199).  Here  a  parcel  of  gas  is  considered,  which  con¬ 
tains  a  great  number  of  "particles"  (molecules),  which  move  at  random  due  to 
their  temperature.  The  density  of  the  considered  type  of  particles  in  the  par¬ 
cel  is  NQ  ;  outside  of  it  this  density  is  zero  at  the  time  -  =  0. 


The  parameter 


N„ 


J  Ai r  D i 


(104) 


indicates  the  number  of  particles  of  the  considered  type,  which  have,  after  a 
time  t  has  passed,  displacement  components  along  the  x-axis  in  the  magnitude 
between  x  and  x  +  dx  .  The  square  root  of  the  average  squared  displacement 
along  the  x-axis  derived  from  eq.  (103)  is: 


(105) 
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and  the  average  displacement  is 


x  =  \J  4Dt/7T  .  (106) 

By  using  the  diffusion  coefficients  from  table  208  of  section  1.  6.4,  we  calculate 
for  the  different  heights  the  average  displacement  per  second  (see  the  following 
table): 


AVERAGE  DISPLACEMENT  OF  IONS  PER  SECOND 
DUE  TO  GAS  KINETIC  DIFFUSION 


With  a  wind  velocity  equal  to  that  of  sound,  ions  travel  downstream  with  an  av¬ 
erage  velocity  of,  for  example,  300  m/sec.  That  means  that  they  pass  the  zone 
of  rings  at  the  upstream  end  of  the  modified  GERDIEN  chamber,  being  about 
3  cm  long,  in  10~^  sec,  and  J 10-4=  10~2  .  With  the  same  velocity,  which  would 
'apply  to  the  uppermost  ranges  of  falling  height,  the  ions  would  traverse  the  whole 
modified  GERDIEN  chamber  of  about  1  m  in  1/300  second,  the  square  root  of 
which  is  5.  8  x  10"^  sec’1'2.  The  values  of  the  table  for  the  range  of  80  km  must 
be  multiplied  by  these  square  roots  of  time  in  order  to  determine  the  effect  of 
diffusion  on  (a'  ion  losses  in  the  ring  zone,  and  (b)  ion  path  broadening  inside 
the  chamber 

For  an  application  of  the  GERDIEN  chamber  in  lower  heights,  at  10  km,  the 
rings  must  be  longer,  for  example,  12  cm,  and  the  wind  velocity  may  be  as 
small  as  5  m/ sec. 

This  leads  to  square  root  of  time  values  of  4.  9  x  10~^  sec^^or  the  ring  zone, 
and  4.  45  x  10"  for  the  chamber  proper.  Thus  we  get  the  actual  average  dis¬ 
placement  as  indicated  in  the  following  table: 
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EXTREME  CASES  OF  AVERAGE  DISPLACEMENT 


x  for  ring  zone 

x  for  chamber  proper 

[m] 

[m] 

at  10  km 

1.  6  x  10-4 

1.47  x  10“3 

at  80  km 

4.96  x  10“  3 

2.  28  x  10“2 

The  true  displacement  occurs  due  to  a  MAXWELLIAN  curve.  From  the  last 
table  we  see,  that  gas  kinetic  diffusion  creates  a  small  concern  for  the  height  of 
10  km,  but  must  be  taken  into  account  at  80  km.  In  fact,  it  prevents  another  li¬ 
mitation  fGr  the  application  of  the  modified  GERDIEN  chamber  in  altitudes  ap¬ 
proaching  80  km;  if  the  distance  between  the  rings  is  10  mm,  and  the  displace¬ 
ment  according  to  the  table  is  in  the  average  about  5  mm,  only  half  of  the  ions 
will  pass,  roughly  ^speaking,  and  the  resolution  of  mobilities  at  the  lower  end 
of  the  range  will  become  rather  poor  at  that  extreme  height. 

In  conclusion,  all  our  different  approaches  lead  to  the  fact,  that  for  80-km  di¬ 
ffusion  disturbances  are  severe,  and  that  post-calibration  must  investigate  the 
influence  of  diffusion  at  altitudes  of  70  and  60  km  or  less. 

The  "electric  diffusion"  caused  by  the  mirror-charge  effect  leads  to  the  follow¬ 
ing  facts.  If  we  again  take  a  distance  between  ion  and  wall  in  the  order  of  1mm, 
the  electric  field  generated  by  the  mirror-charge  effect,  is  in  the  order  of 
mV/m  With  the  mobility  values  for  altitudes  between  ground  and  80  km,  the  re¬ 
sulting  ion  velocity  towards  the  wall  is  in  the  order  of  a  few  mm  per  second  at 
80  km,  and  of  less  than  a  thousandth  of  that  at  10  km.  Thus,  the  effect  of  the 
electric  diffusion  can  be  neglected;  in  lower  altitudes  it  does  not  play  a  role  at 
all,  and  at  the  maximum  altitude,  where  it  may  begin  to  be  noticeable,  its  ef¬ 
fect  is  much  smaller  than  the  effect  of  gas-kinetic  diffusion. 

The  most  severe  influence  is  that  of  eddy  diffusion,  Eddy  diffusion  occurs  as 
soon  as  the  flow  of  the  air,  containing  ions,  is  turbulent.  For  the  lower  atmos¬ 
phere,  GUNN  (1964)  has  discussed  this  effect. 

Losses  of  about  39  percent  of  the  small  ions  occurred  when  the  CARNEGIE  ex¬ 
pedition  used  a  special  type  of  GERDIEN  chamber.  In  this  chamber,  eddy  dif¬ 
fusion  was  produced  at  the  upstream  end,  because  coaxial  rings,  40  cm  long, 
and  about  2  cm  spaced  were  inserted  here.  GUNN  does  not  discuss  to  what  ex¬ 
tent  this  changed  the  mobility  spectrum  of  the  incoming  ions.  While  in  gas  ki¬ 
netic  diffusion  the  loss  is  a  function  of  mobility,  the  losses  in  eddy  diffusion 
presumably  depend  on  mobility  to  a  much  less  degree.  However,  this  degree 
is  not  known  and  cannot  be  determined  theoretically.  For  this  reason,  it  is 
planned  to  avoid  turbulence  during  the  measurements  with  the  modified  GERDIEN 
chamber. 
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The  probable  diffusion  losses  will  be  determined  after  each  flight.  Since  these 
investigations  belong  to  the  velocity-density  dependent  complex,  they  will  be  in¬ 
cluded  in  the  post-calibration  periods. 


1.  7.  3  FIELD- STRENGTH  DEPENDENT  MOBILITY  UNDER  LOW  PRESSURE 

If  the  ions  can  gain  much  energy  between  collisions,  i.  e. ,  when  moving  in  a  gas 
of  low  density  under  the  influence  of  a  great  electric  field,  their  velocity  is  no 
more  a  linear  function  of  the  field  strength 


v;  ~  E  ; 

and  it  finally  becomes  dependent  on  the  square  root  of  the  field 


(107) 


vj  ~  E1/2  . 


(108) 


Eq.  (107)  applies  when  the  kinetic  energy  acquired  by  the  ions  from  the  elec¬ 
tric  field  is  small  in  comparison  to  their  thermal  energy: 


le  E  «  kB  T 


(109) 


where 

l  =  mean  free  path 

E  =  electric  field  strength 

kg  =  BOLTZMANN  constant 

T  =  temperature  in  degrees  KELVIN 

e  =  electric  charge  of  the  ion. 

Energy  transfer  between  neutral  particles  (molecules)  and  ions  is  small  when 
the  respective  masses  are  quite  different.  If  this  is  the  case,  the  ions  may 
store  energy  through  many  collisions: 


e  E  Z  «  kg  T 


(110) 


where  and  mj  are  the  masses  of  neutral  molecules  and  ions,  respectively. 
This  may  be  expressed  also  as  follows,  substituting  the  mean  free  path  by  the 
collision  cross  section: 


m2  ml  \ 

-  +  - J  e  E  «  pa 

ml  m2  / 


(HI) 
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where 


p  =  pressure 

a  =  collision  cross  section. 

For  mj »  and  introducing  the  magnitudes  for  the  other  parameters,  we  get 

E  V/m  (1121 

-  «  200  -  . 

p  torr 

The  limit  determined  in  eq.  (112)  is  never  surpassed  in  the  free  atmosphere, 
since  the  atmospheric  electric  field  strength  is  small  where  the  pressure  is  low. 
However,  in  the  GERDIEN  chamber,  stronger  fields  have  to  be  applied.  Since 
it  is  a  cylindrical  condenser,  the  limiting  field  with  decreasing  pressure  will  at 
first  be  reached  close  to  the  inner  (driving)  electrode.  Even  then,  the  influence 
on  the  measurement  would  be  small,  for  three  reasons:  (1)  the  transition  be¬ 
tween  the  regions  of  field-independent  and  field-dependent  mobilities  is  a  smooth 
one;  (2)  the  surpassing  of  the  critical  E/p  value  is  smooth  as  well,  and  (3)  ions 
move  mostly  close  to  the  outer  electrode  where  the  transition  occurs  at  last. 
However,  part  of  the  ion  path  at  altitudes  higher  than  60  to  70  km  will  certainly 
be  affected.  Thus,  the  eveluation  of  measurements  requires  a  consideration  of 
this  effect. 

Based  on  the  studies  of  G.  H.  WANNIER  and  L.  B.  LOEB,  the  following  formula, 
which  describes  the  movement  of  ions  in  the  cylindrical  condenser,  has  been 
derived: 


where 

=  arc  cos  \f  cos  2  cot 

p  =  distance  from  the  center  of  the  modified  Gerdien  chamber 
m2  =  molecular  mass 
(lq  =■  ionic  mass 
e  =  ionic  charge 
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U  =  voltage  applied  between  the  two  electrodes  of  the  modified  Gerdien 
chamber 

ii2  =  number  density  of  the  molecules 
a  =  collision  crose- section  of  ions  and  molecules 
R  =  radius  of  outer  electrode  of  the  Gerdien  chamber 

r  =  radius  of  inner  electrode  of  the  Gerdien  chainber 

E  =  elliptical  integral  of  second  order 
F  =  elliptical  integral  of  first  order 

Cj  =  integration  constant 

This  formula  becomes  more  and  more  inaccurate  above  a  certain  frequency  li- 
mit--in  the  Earth's  atmosphere  about  1  MHz  at  30  km  and  10  kHz  at  80  km. 

1.  7.  4  CHANGES  OF  THE  NATURE  OF  IONS  BY  DRIVING  VOLTAGE 

The  possibility  that  ions  exist  which  can  be  destroyed  by  the  application  of  a 
very  small  energy  has  been  indicated  in  section  1.  2.  2.  If  ions  could  really  ex¬ 
ist  with  binding  energies  as  low  as  0.  014  eV,  we  certainly  could  destroy  them 
in  the  GERDIEN  chamber.  Thus,  if  we  do  not  find  mobilities  which  may  be  at¬ 
tributed  to  such  ions,  '■he  existence  of  these  ions  is  not  yet  refuted  experimen¬ 
tally.  However,  the  probability  is  small. 

The  existence  of  ions,  on  the  other  hand,  is  very  probable.  Their  electron 
affinity  is  0.  9  ±  0.  1  eV  (NAWROCKI  and  PAPA,  1961),  i.  e. ,  it  is  smaller  than 
that  of  ail  other  probable  constituents. 

Here,  the  problem  exists,  whether  the  electric  field  in  the  GERDIEN  chamber 
could  destroy  the  ions  and  thus  change  the  mobility  values  prior  to  the  measure¬ 
ment  itself.  In  order  to  do  this,  the  electric  field  must  become  so  great  that 
it  convoys  more  energy  to  the  ions  than  they  have  due  to  their  own  thermal  mo¬ 
tion  (see  section  1.  7.  3.  ),  Thus,  there  is  a  common  lower  limit  for  the 
"WANNIER  r  _'ect"  and  the  possible  ion  destruction. 

For  a  closer  estimate,  we  considei’  the  occurrence  of  long  free  paths  as  plotted 
in  figures  22  and  23  of  section  1.  6.  2.  Supposing  a  maximum  field  strength  in¬ 
side  the  GERDIEN  chamber  of  110  V/m  exists  for  measurements  in  70  or  80  km, 
and  of  7.  6  kV /m  in  30  km,  we  can  easily  calculate  that: 
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at  80  km  about  10  percent 
at  70  km  about  10"^  Percent 
at  30  km  less  than  10" 4  Percent 
of  the  ions  may  acquire  energies  of  0.  8  eV  or  more. 

1.  7.  5.  SHIELDING  BY  A  SPACE  CHARGE  LAYER 

From  the  domain  of  electrochemistry,  the  concept  of  shielding  off  any  electric 
field  by  a  layer  of  space  charges  of  opposite  sign  in  front  of  any  electrode  has 
been  introduced  into  the  domain  of  atmospheric  physics.  The  general  applica¬ 
bility  of  this  concept  is  still  vague.  However,  the  so-called  DEBYE-length  (or 
DEBYE  shielding- distance): 


--with  ne  =  electron  number  density--is  often  quoted.  We  have  to  show,  that 
the  possibility  of  a  shielding  effect  of  this  kind  does  not  disturb  our  measure¬ 
ments. 

It  is  intended  to  keep  the  field  outside  the  Modified  GERDIEN  chamber  free  of 
any  disturbing  electric  fields  which  may  be  generated  from  the  instrument.  The 
first  set  of  rings  (Nos.  00  to  07  in  Model  A)  is  constructed  for  this  purpose. 

That  means  that  there  is  no  field  which  could  build  up  a  shielding  space  charge 
layer  by  the  electrode  effect. 

Inside  the  GERDIEN  chamber,  the  electron  density  is  practically  zero,  since 
all  electrons  which  are  not  removed  by  diffusion  between  the  rings  will  imme¬ 
diately  be  destroyed  by  the  electric  field  in  the  chamber.  A  shielding  effect  of 
the  ions  caused  by  the  electrode  effect,  however,  certainly  does  not  occur  over 
so  short  a  distance  as  given  in  the  chamber, 

1.  7.  6  ACQUISITION  OF  AN  OWN  CHARGE  BY  THE  INSTRUEMENT 

Basically  there  are  two  possibilities  that  the  whole  modified  GERDIEN  cham¬ 
ber,  including  everything  which  is  electrically  connected  to  it,  may  adopt  an 
electric  charge.  This  is  a  well  known  event  at  higher  altitudes.  If  this  were 
so  in  our  case,  the  intake  of  ions  would  be  a  selective  one,  depending  on  the  sign 
and  even  on  mobility  of  the  aeroions.  It  can  be  shown,  however,  that  the  pro¬ 
bability  of  such  a  charge-up  process  is  very  small. 

An  external  charge-up  would  occur  if  more  ions  of  one  sign  would  impact  on  the 
outside  of  the  instrument.  This  already  is  unlikely  to  have  a  great  effect,  since 
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it  is  intended  to  move  the  instrument  at  a  much  slower  rate  than  the  thermal 
velocities  of  both  ions  and  electrons.  And  even  if  this  occurs,  the  instrument, 
together  with  the  carrier,  is  so  much  greater  than  the  mean  free  path  that  the 
concept  of  conductivity  applies  to  its  full  extent.  This  has  the  consequence  that 
the  rather  great  surface  of  the  instrument  will  be  discharged  very  quickly,  once 
it  has  obtained  a  potential  difference  with  respect  to  the  environmental  atmos¬ 
phere. 

An  internal  charge-up  would  occur,  if  inside  the  chamber  not  all  ions  are  in¬ 
tercepted.  That  means,  if  ions  of  predominantly  one  sign  would  exit  at  the  cham¬ 
ber's  rear  end  with  the  airstream,  a  certain  charge-up  may  occur.  Again,  con¬ 
ductivity  is  so  great  that  the  potential  will  very  quickly  be  equalized.  If  so  wan¬ 
ted,  the  aft  rings  of  the  chamber  (30  to  37  in  Model  A)  may  be  subdivided  into 
two  sets.  Then,  there  are  several  easy  possibilities  to  establish  a  workable 
ion  fence  or  filter,  and  no  ions  would  exit  any  more. 

1.  7.  7  DISTURBANCES  BY  COSMIC  PARTICLES 

v 

There  is  basically  a  possibility  that  the  internal  space  of  the  modified  GERDIEN 
chamber  will  show  an  additional  ionization  due  to  cosmic  rays.  The  wails  of  the 
chamber  are  thick  and  partially  built  by  heavy  materials,  and  cosmic  rays  pen¬ 
etrating  them  may  give  rise  to  secondary  effects  which,  in  turn,  would  ionize 
many  neutral  particles  in  the  chamber. 


For  a  first  approximation,  we  compare  our  instrument  with  ionization  chambers 
used  for  the  measurement  of  cosmic  rays  in  the  upper  atmosphere.  These  ioni¬ 
zation  chambers  must  be  constructed  in  such  a  way  as  to  avoid  the  indicated  dis¬ 
turbances.  The  ionization  chamber  described  by  NEHER  (1953)  and  used  by  him 
still  in  1963  (NEHER  &  ANDERSON,  1964)  was  essentially  a  sphere  of  25  cm 
diameter,  made  of  steel,  with  a  thickness  of  500/^m,  This  gives  a  mass  thick¬ 
ness  of  3.  9  kg/m'-.  The  chamber  has  been  used  to  altitudes  of  about  30  km,  and 
an  air  mass  of  100  kg/m^  was  still  overhead  of  it. 


In  the  modified  GERDIEN  chamber,  Model  Bl,  the  receiving  electrodes  are 
made  of  steel  with  a  thickness  of  250  (im,  However,  there  are  aluminum  beams 
to  support  the  structure  against  the  expected  great  accelerations,  and  in  the  most 
unfavorable  case,  a  ray  may  penetrate  25  mm  of  A1  before  entering  the  chamber, 
which  corresponds  to  a  mass  of  67.  5  kg/m^.  If  secondary  particles  will  have 
but  small  angles  with  respect  to  the  primary  ray,  only  a  small  fraction  of  the 
chamber's  volume  would  be  affected.  To  obtain  still  another  number,  one  may 
divide  the  chamber's  whole  mass,  9.  1  kg,  by  its  surface,  5.  6  x  10"  m^,  and 

arrive  at  1.  63  kg/m^,  which  is  less  than  the  3.  9  kg/m^  of  NEHER's  ionization 
chamber. 


It  is  difficult  to  draw  a  conclusion  from  these  numbers.  In  addition,  it  must  be 
remembered  that  the  100  kg/m^  of  air  above  NEHER's  chamber  are  practically 
not  existing,  if  the  modified  GERDIEN  chamber  is  to  measure  at  a  70-km  alti- 
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tude.  Thus,  the  total  material  mass  would  be  slightly  smaller  or  much  smaller 
in  the  case  of  the  GERDIEN  than  in  the  case  of  NEHER's  ionization  chamber. 

1.  7.  8  OTHER  LIMITATIONS  AND  DISTURBANCES 

Limitations  of  a  more  technological  nature  are  given  by  the  fact,  that  a  velocity 
below  that  of  sound  could  hardly  be  maintained  in  altitudes  much  higher  than  80 
km.  Also,  the  rather  high  frequency  of  the  driving  voltage  which  then  v/ould  be 
necessary  would  pose  problems  with  respect  to  the  generation  of  the  ion  gate 
pulses  and  the  amplification.  The  driving  voltage  cannot  be  lowered  under  a  li¬ 
mit  given  by  VOLTA  potential  differences,  which  cannot  be  avoided. 

The  low  temperature  of  the  atmosphere  around  the  instrument,  which  comes  up 
from  lower  layers,  may  cause  difficulties  for  the  required  high  insulation.  Ice 
deposits  on  the  insulating  surfaces,  however,  may  be  precluded  by  heating  up 
the  atmosphere  at  the  surface  of  the  insulators  to  about  10°C  above  the  temperat¬ 
ure  of  the  atmosphere. 

The  most  obvious  limitation  of  the  application  of  the  instrument  is,  of  course, 
given  by  the  fact  that  the  requirement  of  many  collisions  between  the  ion  and 
neutral  particles  cannot  be  satisfied  more  where  the  mean  free  path  comes  into 
the  order  of  magnitude  of  the  instrument  dimensions.  This  has  been  treated  at 
different  locations  in  the  preceding  chapters. 

1.7.9  CONCLUSION 

A  general  upper  limit  for  the  application  of  the  method  of  the  modified  GERDIEN 
chamber--and  that  means  for  the  possibility  of  measuring  the  mobility  spectrum 
in  general--cannot  be  given.  Partially,  the  values  of  involved  parameters  are 
not  yet  known,  and  information  is  to  be  expected  from  our  intended  measurements. 
Partially,  the  upper  limit  depends  on  time  variable  atmospheric  features,  and 
no  general  number  can  be  quoted.  Anc},  in  part,  the  limit  depends  on  the  actual 
measurement  itself  inasmuch  as  a  certain  combination  of  wind  velocity  in  the 
chamber  and  density  of  pie  atmosphere  and  total  ion  number  density  determines 
whether  a  measurement  still  can  be  evaluated.  This  combination  will  not  be 
known  prior  to  the  measurement, 

To  give  a  rough  estimate,  we  may  say  that  no  reason  can  be  seen  why  measure¬ 
ments  up  to  about  60  km  should  not  give  reliable  results.  However,  the  accu¬ 
racy  at  60  km  may  already  be  slightly  smaller  than  at  55;  50,  40,  or  30  km. 
Furthermore,  there  is  hope,  that  reliable  information  may  be  obtained  up  to 
70  or  75  km,  while  values  from  a  height  of  80  km  probably  cannot  be  evaluated 
further  to  obtain  a  mobility  spectrum.  However,  some  information  still  may 
be  derived  from  these  measurements  ...t  the  peak. 
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PART  2 

LABORATORY  EXPERIMENTS 
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PRECEDING  PAGE  BLANK 


2.  1  LOW  PRESSURE  WIND  TUNNEL 


2.  1.  1  INTRODUCTION 

The  wind  tunnel  to  be  described  in  the  following  pages  has  been  constructed  for 
the  simulation  of  the  atmosphere  between  ground  and  about  80  km,  as  seen  from 
a  body  moving  with  velocities  between  zero  and  about  the  speed  of  sound.  These 
limits  have  been  chosen  for  the  following  reasons. 

For  many  in  situ  measurements  in  the  free  atmosphere  it  is  desirable  to  avoid 
the  shock  wave  which  is  generated  by  a  body  moving  with  a  speed  greater  than  that 
of  sound  (Mach  1).  Such  shock  waves  may  alter  the  atmospheric  parameters  to  be 
measured  -  at  higher  speeds,  even  create  additional  ionization  -  or  disturb  the 
measurement  otherwise.  A  measuring  platform  moving  through  the  atmosphere 
with  subsonic  speeds  may  be  provided  by  the  application  of  modern  forms  of  para¬ 
chutes  for  atmospheric  heights  below  about  80  km.  Since  the  mean  free  path  at  that 
height  comes  in  the  order  of  the  geometric  dimensions  of  many  practicable  instru¬ 
ments,  and  since  it  is  there  where  the  atmospheric  conditions  undergo  a  typical 
change  (the  mesopause),  ‘the  height  range  between  ground  and  about  80  km  may  be 
considered  a  more  or  less  uniform  environment  for  the  development  of  in  situ  in¬ 
strumentation.  The  possibility  of  providing  subsonic  falling  speeds  and  the  criti¬ 
cal  length  of  the  mean  free  path  occurs  in  the  atmospheres  of  other  planets  at  other 
altitudes,  but  again  they  constitute  the  upper  limit  of  environmental  conditions  fav¬ 
orable  for  the  application  of  many  instruments.  An  additional  requirement  for  the 
wind  tunnel  was  given  by  the  fact  that  the  simulation  of  the  environment  had  to  be 
maintained  for  times  comparable  to  measuring  times  in  the  free  atmosphere,  which 
often  will  count  by  thousands  of  seconds  or  more. 

The  actual  wind  tunnel  in  its  present  configuration  allows  for  the  fulfillment  of  these 
requirements  in  a  continuous  operation  which  may  last  for  hours:  wind  speeds  be¬ 
tween  zero  and  more  than  350  m/.sec  for  any  desired  pressure  between  about  10 
and  10-2  tori*,  and  to  a  certain  extent  for  higher  pressures  as  well.  It  is,  in  add¬ 
ition,  possible  to  exceed  these  velocities  and  to  lessen  the  pressure  to  a  certain 
dqgree.  For  a  short -time -measurement  the  wind  speed  can  be  increased  consider¬ 
ably:  the  tank  of  about  300  m^  content  pumped  down  to  about  10"^  torr  then  serves 
as  a  surge  action  vacuum  of  strong  capacity. 

Provisions  are  made  to  simulate  the  ion  content  of  the  atmosphere  as  well,  and 
temperature  simulation  may  be  added  if  so  desired. 

The  test  section  of  the  wind  tunnel  consists  of  a  pipe  with  a  variable  diameter.  For 
our  experiments  the  diameter  of  the  pipe  varies  between  10  and  20  cm.  However, 
this  diameter  as  well  as  a  number  of  other  features  can  easily  be  changed  if  the 
system  is  to  be  adapted  to  particular  experiments.  The  construction  of  this  wind 
tunnel  has  been  made  necessary  for  the  development  and  calibration  of  the  modified 
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GERDIEN  chamber,  intended  to  measure  in  altitudes  between  30  and  80  km,  because 
we  have  been  unable  to  find  another  facility  which  would  have  met  the  quoted  require¬ 
ments. 

It  is  not  surprising  that  no  other  facility  of  this  kind  exists.  There  are  three  main 
difficulties  to  overcome:  (1)  mechanical  pumps  have  but  a  small  effectiveness  be¬ 
low  a  pressure  of  about  1  torr  and  diffusion  pumps  cannot  be  used  for  pressures 
higher  than  a  few  tenths  of  a  torr.  Thus,  rather  great  mechanical  pumps  in  connec¬ 
tion  with  a  great  tank  are  the  only  means  of  providing  sufficient  wind  speed  in  the 
whole  range  of  pressure;  (2)  the  maintainance  of  small  wind  velocities,  i.  e.  ,  of  a 
small  constant  and  unidirectional  air  flow,  requires  a  high  degree  of  tightness  of 
the  whole  system,  higher  than  is  usually  required  for  this  pressure  range.  This  is 
made  difficult  by  the  necessity  of  providing  many  ports  and  feedtroughs,  and  to  com¬ 
bine  the  tunnel  of  a  rather  great  number  of  different  sections.  This  problem  could 
be  solved.  The  remaining  leak  rate  corresponds  to  a  wind  velocity  of  about  10 
cm/sec,  which  is  much  below  the  limit  of  all  measuring  instruments;  and  (3)  it  is 
difficult  to  measure  small  wind  velocities  under  the  required  low  pressures. 

The  wind  tunnel  is  briefly  described  below,  the  principle  of  its  operations  is  indi¬ 
cated,  the  attached  instrumentation  for  the  determination  of  pressure  and  wind 
speed  is  discussed,  and  information  on  actual  performance  is  given. 

2.  1.2  DESCRIPTION  OF  THE  SYSTEM 

Figure  24  shows  an  overall  drawing  of  the  wind  tunnel  and  Figure  25  is  a  general 
sketch  of  the  overall  wind  tunnel.  The  wind  tunnel  (figure  26)  js  attached  directly 
to  a  big  tank  (figure  2/)  whereas  the  pump  system  (figure  ^8)  connected  at  the 
opposite  end  of  this  tank.  Figure  29  shows  a  cross  section  of  the  wind  tunnel.  Pres¬ 
surized  air  or  gas  is  inserted  in  the  wind  tunnel  through  valves.  The  inside  of  the 
wind  tunnel  is  nickel -plated  so  that  during  the  operation  of  the  wind  tunnel  no  dust 
can  be  created  by  the  joining  air  stream.  The  jet  of  air  is  smoothed  out  by  the 
bafflings  in  the  first  section  before  entering  the  test  section.  The  aging  section 
can  be  used  for  a  temperature  control  of  the  air.  Presently  it  is  used  for  the  gen¬ 
eration  of  ions.  The  test  section  as  used  in  figure  29  is  60  cm  long  and  has  a  di¬ 
ameter  of  20  cm.  It  is  accessible  through  a  rectangular  port.  The  smoothing  sec¬ 
tion  and  the  baffling  section  (II)  prevent  the  air  from  being  disturbed  in  the  test 
section  when  valves  (C)  are  opened. 

By  changing  the  openings  of  valves  A  and  valves  C  and  by  using  different  valves  B 
in  the  bleeding  section  ,  the  whole  pressure -velocity  range  can  be  obtained.  Three 
sits  of  pumps  are  used  (figures  24,  25,  and  28).  The  first  set  of  pumps  consists 
ox  three  11  kw  pumps  in  parallel,  the  second  of  the  three  8  kw  pumps  in  parallel. 

For  the  third  stage,  two  37  kw  pumps  are  used.  The  air  removal  of  the  pumps 
amounts  to  about  15m^/sec  at  10”*  torr.  A  short  pipe  with  a  60-cm  diameter 
connects  the  pipes  to  the  tank. 
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AUXILIARY  SYSTEM 


PRESSURE  WIND  TUNNEL  FOR  CONTINUOUS  OPERATION 
(WITH  PUMPING  SYSTEM) 


SUBSONIC  LOW 
PRESSURE  WIND 
20ci<>  TUNNEL 
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Figure  25  OVERALL  SYSTEM  OF  THE  SUBSONIC  LOW. 
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Figure  29  CROSS  SECTION  THROUGH  LOW-PRESSURE  WIND  TUNNEL  FOR 
CONTINUOUS  OPERATION  SHOWING  SECTIONS 
AND  VALVES 


The  tank  to  which  the  pumps  are  connected  has  a  volume  of  283  m3.  One  purpose 
of  this  big  tank  is  to  smooth  out  all  pulsations  of  the  pump  system.  In  fact,  a 


steady  air  stream  is  obtained.  The  pumps  are  operative  to  below  10 
a  blank  flange,  they  will  pump  the  big  tank  down  to  about  10“3  torr. 
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torr.  Against 


A  smaller  pump  system  may  be  used  in  addition  (figure  25).  Here,  the  air  removal 
amounts  to  about  1.  5  m3/sec  only,  and  the  lowest  obtainable  pressure  is  about 
6  x  10"2  torr. 


2.1.3  OPERATION  OF  THE  SYSTEM 

In  a  fully  developed  flow  in  a  pipe,  either  laminar  or  turbulent,  the  velocity  of  a 
gas  is  given  by  a  function  $  of  the  following  form  (ROHSENOW  and  CHOI,  1961): 

v  =  4>  (pj>  ?2,  P >  R> 


where 

v  =velocity  of  the  gas  inside  the  tube 
Pj  =pressure  at  the  front  end  of  the  pipe 
P2  =pressure  at  the  rear  end  of  the  pipe 
R  resistance  of  the  pipe 

p  =density  of  the  gas 

v  =kinematic  viscosity  of  the  gas 
e  internal  energy  of  the  gas  per  unit  mass. 

Assuming  laminary  flow  in  the  pipe,  the  velocity  distribution  in  a  cross  section 
is  given  by 


Pl  “  P2  7  7 

v  (y)  =  — - - —  (yr2  -  y2); 

4  v  p  l 


(116) 


where  .y  =distance  from  the  axis  of  the  pipe 
yr  =the  radius  of  the  pipe 

The  density  of  the  gas  (air)  in  the  pipe  in  the  middle  between  the  two  ends  is  as 
follows: 

Pl  -  P2 


where  m  =the  mass  of  a  gas  molecule  (117) 

kB  =the  BOLTZMANN  constant,  1.  38  x  10"^  Ws/degree 
T  =the  temperature  of  the  gas. 

By  eqs.  116,  and  117,  it  is  shown  that  with  independent  control  of  the  pressures 
pl  and  p2  ,  i.  e.  ,  the  pressures  at  both  ends  of  the  pipe,  gas  density  and  gas  vel¬ 
ocity  in  the  pipe  can  both  be  controlled  independently  within  certain  limits. 
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This  is  the  basic  idea  for  the  operation  of  the  wind  tunnel:  by  maintaining  a  cer¬ 
tain  pressure  in  a  vessel  at  the  upstream  end  of  the  tunnel  (i.  e.  ,  baffling  section), 
and  maintaining  another  pressure  at  the  downstream  end  (bleeding  section),  the 
desired  values  for  density  and  wind  velocity  can  be  set.  There  are  five  practical 
means  of  obtaining  this  result:  (1)  the  pressure  in  the  gas  container  (P^  in  figure  29); 
(2)  the  set  of  valves  A;  (3)  the  set  of  valves  B;  (4)  the  set  of  valves  C;  and  (5)  the  gas 
removal  speed  of  the  pumps.  Because  of  technological  reasons,  the  gas  removal 
speed  of  the  pumps  cannot  be  controlled  at  will,  but  the  remaining  four  regulation 
possibilities  are  sufficient. 

In  fact,  the  stability  of  the  system  is  such,  that  a  point  in  the  density -velocity  dia¬ 
gram,  set  on  Friday  afternoon,  is  again  obtained  when  the  system  is  put  to  opera¬ 
tion  on  Monday  morning. 

To  obtain  the  gas  resistance  between  the  pumps  and  the  wind  tunnel  proper,  the 
pressure  has  been  measured  at  two  points  of  the  system  with  the  following  result: 


p ressure  at 

connection 

tank  -  wind  tunnel 

6 

4 

3.  5 

11 

15 

19 

millitorr 

pressure  at  pump 

4.  2 

3.  8 

3.  4 

6.  6 

10 

1 5 

millitorr 

quotient  of  both 

1.  4 

1.  08 

1.  02 

I.  67 

1.  5 

1.  26 

That  means  that  the  flow  resistance  betweer.  '.he  wind  tunnel  and  the  pumps  is  small. 

During  the  experiments,  the  gas  pressure  is  monitored  continuously  at  different 
points  of  the  wind  tunnel,  as  indicated  in  figure  8. 

2.  1.4  DETERMINATION  OF  PRESSURE  AND  VELOCITY 

For  the  determination  of  the  wind  velocity  in  the  tunnel,  Pitot  tubes  are  used. 
Depending  on  the  range,  Pitot  tubes  for  subsonic  flow  with  a  very  short  time  con¬ 
stant,  or  Pitot  tubes  for  supersonic  speeds  with  a  somewhat  larger  time  constant 
(shown  in  figure  30)  are  applied.  In  the  smoothing  section,  a  Pitot  tube  for  static 
pressure  is  fixed,  and  a  Pitot  tube  for  the  total  pressure  can  be  moved  by  a  motor 
across  a  diameter  of  the  tunnel.  Another  set  of  Pitot  tubes  has  been  applied  in¬ 
side  the  test  section  itself  during  some  of  our  experiments.  The  pressure  differ¬ 
ences  between  the  two  tubes  of  each  set  of  Pitot  tubes  is  measured  by  a  "Baratron" 
differential  pressure  gauge  type  77,  equipped  with  several  heads  of  different  sen¬ 
sitivity.  The  smallest  pressure  difference  to  be  read  is  about  3  x  10"~’  torr.  Its 
operation  principle  is  a  variable  condenser  (capacitive  sensor)  in  a  bridge  circuit 
and  in  a  heated  atmosphere;  phase  sensitive  ac  voltage  and  polarized  dc  voltage  are 
provided  for  recording  and  controlling. 
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The  static  pressure  is  measured  simultaneously  at  six  different  locations  along 
the  wind  tunnel  (figure  31)  mostly  by  thermocouple  gauges.  For  a  more  accurate 
determination  of  static  pressure,  an  .  Alphatron  absolute  pressure  gauge,  model 
530,  is  attached  to  the  smoothing  section.  It  allows  the  measurement  of  pressures 
from  760  to  less  than  10“^  torr  with  an  accuracy  of  ±  5  per  cent  of  full  scale  in 
the  0.  . .  1  millitorr  range  and  with  linear  response.  A  STOKES'  gauge  is  attach  d 
to  the  smoothing  section  as  well.  Figure  32  shows  the  instrument  panel  with  the 
pressure  gauges,  valves  A  and  C,  and  part  of  the  electronic  instrumentation  for 
a  particular  experiment.  The  measuring  head  of  the  Alphatron  and  one  measur¬ 
ing  head  of  the  Baratron  a^e  visible  in  figure  26,  sitting  on  the  smoothing  section 
of  the  tunnel. 

2.1.5  MEASUREMENT  OF  STATIC  PRESSURE 

For  the  measurement  of  static  pressure  (and  wind  speed)  in  the  axis  of  the  tunnel, 
a  cone  of  internal  diameter  of  5mm  was  fixed  close  to  the  axis  in  the  smoothing 
section  and,  in  some  of  the  experiments,  in  the  test  section  as  well.  This  rela¬ 
tively  large  diameter  was  chosen  for  measurements  in  the  subsonic  range  on  the 
basis  of  experiments  to  obtain  a  short  response  time.  The  central  positioning  of 
the  cone  helped  to  avoid  errors  due  to  radial  changes  of  static  pressure  in  the 
range  of  about  10~^  torr  or  below. 

The  exact  correction  value  for  the  applied  cone  are  not  known.  However,  from 
measurements  with  similar  cones  (TALBOT,  1954),  it  was  derived  that  the  remain¬ 
ing  error  at  about  10“^  torr  would  not  exceed  15  per  cent  (TALBOT,  1954;  ENKENHUS, 
1957).  Including  the  possible  error  of  the  Baratron  instrument,  the  total  error  for 
10-2  torr  and  about  Mach  1  comes  to  20  per  cent. 

It  decreases  rapidly  with  increasing  pressure;  e.  g.  ,  for  Mach  1  and  10~*  torr  it 
is  already  below  5  per  cent. 

2.  1.  6  MEASUREMENT  OF  IMPACT  PRESSURE 

As  mentioned  above,  the  open-end  impact  probes  can  be  moved  by  motors  in  radial 
direction  across  the  tunnel.  The  diameter  of  the  orifice  is  about  5  mm.  The  exter¬ 
nal  diameter  of  the  probe,  being  about  7  mm,  was  taken  as  a  characteristic  length 
for  the  determination  of  the  Reynolds  number.  With  Reynolds  numbers  greater  than 
100,  the  continuum  flow  theory  can  be  applied,  and  errors  introduced  by  slip  flow 
or  molecular  flow  can  be  neglected  (ENKENHUS,  1957).  Smaller  Reynolds  numbers 
occur  at  lower  pressures.  At  about  2  x  10""1  torr  (Reynolds  number  50),  the  total 
error  due  to  outgassing  and  viscosity  is  estimated  to  be  10  per  cent  for  a  wind  speed 
of  Mach  1.  At  still  smaller  pressures,  the  error  increases  rapidly  and  amounts  to 
about  40  per  cent  at  5  x  10~2  torr  and  Mach  1  (Reynolds  number  10).  Corrections 
with  respect  tc  viscosity  have  not  been  applied  because  available  correction  factors 
differ  very  much  and  do  not  seem  to  be  reliable  enough  for  our  particular  task. 

•Outgassing  from  the  material  in  the  pressure  j;auge  can  last  for  many  hours.  Under  rather  low  pressures,  either  it  must 
be  awaited  that  long,  or  the  error  must  be  estimated  and  accounted  for. 
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2.  1.  7  DETERMINATION  OF  FLOW  VELOCITY  (Wind  Speed) 


Two  different  equations  for  the  determination  of  the  flow  velocity  may  be  quoted, 
and  in  both  cases  the  velocity  is  given  by  Mach  numbers  and  not  in  a  temperature- 
independent  scale.  For  subsonic  flow  (M  <  1),  we  have 


A  ,  zri  ,A  y- 

l  2  ) 


(118) 


and  for  supersonic  flow  (M  >  1)  the  RAYLEIGH  formula  may  be 
applied: 


(119) 


In  both  formulas  are 

pt  =  the  total  or  impact  pressure 
p  =  the  static  pressure 

M  =  the  Mach  number  or  ratio  v/c  with  c  =  speed  of  sound 
y  =  the  ratio  of  specific  heating:  Cp/cv 
(for  air  at  20°C,  y  =  1.4). 


These  formulas  are  valid  under  certain  assumptions,  and  their  reliability  decreases 
when  the  pressure  gets  low.  A  study  of  these  assumptions  and  of  the  inaccuracy 
under  the  low-pressure  conditions  has  been  made,  and  it  was  decided,  that  for  the 
purpose  of  our  measurements,  eqs.  (118)  and  (119)  maybe  applied.  The  result¬ 
ing  error  for  the  determination  of  the  wind  speed  comes,  for  wind  speeds  of  Mach  1, 
to  about  10  per  cent  for  altitudes  of  about  60km,  and  to  about  20  per  cent  for  70  km, 
increasing  rapidly  if  altitudes  of  more  than  70  km  are  to  be  simulated.  Since  these 
accuracies  have  been  sufficient  for  our  purposes.no  attempt  was  made  to  improve 
the  measurement  by  the  application  different  types  of  probes  and  the  execution  of 
a  series  of  particular  experiments. 


2.  1.  8  SIMULATION  OF  OTHER  PARAMETERS  THAN  DENSITY  AND  WIND  SPEED 


For  the  simulation  of  tne  ionic  state  of  the  free  atmosphere,  there  are  provided 
two  possibilities  to  create  ions.  To  simulate  the  conditions  up  to  heights  of  about 
40  km,  polonium -foils  are  inserted  in  the  aging  section  of  the  wind  tunnel.  De¬ 
pending  on  the  desired  ion  density,  more  or  less  radioactive  material  is  used.  This 
ion  source  can  be  moved  by  a  motor  along  the  axis  of  the  aging  section.  This  pos¬ 
sibility  allows  the  provision  of  older  or  younger  ions  to  the  test  section;  furthermore 
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it  yields  a  rather  sensitive  regulation  of  the  ion  number  density  since  recombina¬ 
tion  and  diffusion  losses  depend  on  the  traveling  time. 

To  get  a  sufficient  ion  density  for  higher  altitudes,  a  box  containing  (at  present)  400 
millicurie  polonium  is  used.  This  box  is  connected  by  a  thin  tube  and  a  valve  with 
a  pressurized  air  bottle.  During  the  operation  of  the  wind  tunnel,  the  air  stream¬ 
ing  through  this  box  is  highly  ionized,  since  the  air  pressure  is  then  about  atmos¬ 
pheric.  The  ionized  air  is  released  through  a  nozzle  of  4-mm  diameter  into  the 
wind  tunnel  proper.  Iri  this  nozzle,  the  pressure  reduces  to  the  operating  pressure 
maintained  in  the  tunnel.  Many  ions  are  lost  in  the  nozzle  due  to  diffusion,  but  it 
is  still  possible  to  obtain  a  number  density  in  the  test  section,  which  corresponds 
to  a  height  of  70  km.  By  changing  the  pressure  in  the  box  and  the  diameter  of  the  noz 
zle,  a  still  greater  ion  density  could  be  obtained. 

The  wind  tunnel  described  here  is  not  a  recirculating  system.  That  means  that  fresh 
gas  supply  is  possible  throughout  the  operation.  Thus,  every  gas  or  gas  mixture 
which  can  be  provided  from  pressurized  bottles  and  which  does  not  disturb  or  en¬ 
danger  the  functions  of  the  wind  tunnel  and  the  pumps  can  be  used.  However,  it 
must  be  remembered  that  even  the  purest  gases  purchased  from  industry  or  stored 
in  bottles  will  be  contaminated  by  aerosol  particles.  This  fact  may  considerably 
disturb  the  simulation  of  the  conditions  in  the  free  atmosphere.  To  avoid  this  dis¬ 
advantage,  "absolute  filters"  can  be  inserted.  These  filter  out  aerosol  particles 
of  the  0.  5  /<m  size  to  better  than  99.  97  per  cent,  and  particles  of  greater  or  smaller 
size  to  100  per  cent,  practically. 

If  the  ionization  box,  as  described  in  the  preceding  section,  is  applied,  it  may  be 
fed  directly  with  a  gas  other  than  that  used  in  the  wind  tunnel.  In  this  way,  it  is 
possible,  to  a  certain  extent,  to  have  ions  of  one  chemical  nature  in  a  neut.ral  gas 
of  another  one. 

No  provisions  have  been  made,  so  far,,  to  control  the  temperature  or  to  simulate 
in  the  wind  tunnel  the  radiative  environment  given  in  the  free  atmosphere.  Simu¬ 
lation  of  temperature  is  possible  to  a  certain  extent,  and  facilities  for  radiation 
simulation  may  be  installed  within  certain  limits. 

2.  1.  9  EXPERIMENTALLY  OBTAINED  VELOCITY  PROFILES  AND 
MAXIMUM  VELOCITIES 

All  results  which  are  given  in  the  following  are  values  which  have  been  obtained 
during  continuous  operation  within  the  test  section,  which  during  these  experiments 
had  a  diameter  of  10  cm. 

In  figure  33,  as  an  example,  a  few  wind  profiles  are  given  which  have  been  obtained 
at  different  pressures.  Figure  34  shows  in  which  pressure-velocity  range  the  wind 
tunnel  can  be  operated  in  our  present  setup.  Each  point  in  this  figure  gives  the 
value  of  the  maximum  velocity  of  one  wind  profile  inside  the  test  section  and  the 


-112- 


65-12062 


Figure  33  EXAMPLES  FOR  VELOCITY  DISTRIBUTIONS  IN  PRESENT  SET-UP  OF  TEST 

SECTION 

(Each  curve  is  obtained  by  8  to  15  measurements  and  calculations.) 

(1)  static  pressure  0.24Torr,  valves  B21,  B22,  and  C  partially  open 

(2)  static  pressure  0.068  Torr,  valves  B21  and  B22  partially  open 

(3)  static  pressure  0.35Torr,  valves  B1  1  and  C  partially  open 

(4)  static  pressure  0.034  Torr,  maximum  obtainable  speed 

(5)  staric  pressure  0.16Torr,  maximum  obtainable  speed 

(6)  static  p.  sure  0.28  Torr,  maximum  obtainable  speed 


corresponding  static  pressure.  It  can  be  seen  that  pressure  and  velocity  can  be 
regulated  independently  from  each  other.  For  example,  at  a  static  pressure  of 
500  millitorr  (corresponding  to  a  height  of  55  km),  a  number  of  different  wind  speeds 
between  2  and  0,  1  Mach  have  been  simulated.  Even  much  slower  wind  speeds  as 
those  given  in  this  figure  can  be  obtained  at  this  static  pressure.  However,  the 
accuracy  of  the  velocity  measurement  becomes  more  and  more  inaccurate  the  lower 
the  velocity  of  the  streaming  air. 

The  maximum  wind  speed  which  can  be  obtained  with  a  given  test  section  can  be 
determined  from  the  pumping  speed  of  the  pumps.  The  lines  in  figure  34  give  the 
theoretical  obtainable  values.  As  can  be  seen  only  at  pressures  abo^/e  10  torr 
(corresponding  to  a  height  of  about  30  km),  the  obtainable  velocity  is  less  than  0.  3 
Mach.  Higher  velocities  in  this  pressure  range  can  be  obtained  by  operating  an 
additional  pump  system  (see  figure  26). 

As  can  be  seen  from  these  figures,  in  the  present  set-up  of  the  tunnel  and  with  the 
indicated  test  section,  velocities  up  to  about  Mach  I  can  be  obtained  for  a  continous 
operation  in  all  pressure  ranges.  In  the  range  between  10"^  and  10  torr,  about 
Mach  2  can  be  reached. 

Exact  measurements  of  the  higher  velocities  obtainable  in  a  surge  operation  have 
not  as  yet  been  made. 
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2.2  MODIFIED  GERDIEN  CHAMBER,  MODEL  A, 
LABORATORY  MODEL 


2.  2.  1  INTRODUCTION 

It  has  not  been  possible  to  find  a  method  for  measuring  the  wind  speed  in  the 
GERDIEN  chamber  without  introducing  probes  into  it.  The  production  of  visible 
smoke  balls  traveling  in  the  air  stream  ceased  to  be  feasible  for  pressures  be¬ 
low  0.4  torr.  The  introduction  of  probes  into  the  chamber  proper  whould  have 
disturbed  the  electric  field  in  it  and  thus  prohibited  the  intended  measurements. 
For  this  reason,  it  was  necessary  to  construct  a  "dummy  model"  with  exactly 
the  same  dimensions  but  without  the  electrical  parts.  Aerodynamic  probes 
(Pitot  tubes)  have  been  inserted  in  the  dummy  model.  While  a  velocity  profile 
was  obtained  in  the  dummy  model,  a  second  velocity  profile  was  measured  some 
50  cm  downstream  of  it.  Then  the  dummy  model  was  replaced  by  the  chamber 
proper,  and  the  wind  velocity  at  the  second,  downstream  location,  was  set  as 
before.  It  was  then  assumed  that  the  velocity  profile  in  the  chamber  was  the 
same  as  in  the  privious  dummy  model. 

2.2.2  MODEL  A 

The  dimensions  of  the  laboratory  model  (Model  A)  have  been  chosen  on  the  ba¬ 
sis  of  earlier  instruments.  The  receiving  as  well  as  the  driving  electrode  have 
each  been  divided  into  16  highly  insulated  and  gold-plated  parts.  Thus  it  is  pos¬ 
sible  to  exchange  their  functions. 

A  cross  section,  which  is  nearly  accurate  with  respect  to  sizes,  has  been  shown 
in  figure  10  (section  1.4.  1).  The  first  set  of  rings,  00  to  07,  was  added  after 
the  initial  experiments,  because  it  turned  out  that  the  edge  effect  across  the  ion 
gate  (rings  13-14),  created  by  the  potential  differences  between  12  and  14  as  well 
as  between  15  and  13,  increased  the  number  of  entering  ions  of  one  polarity  and 
decreased  that  of  the  other  polarity.  This  was  immediately  avoided,  after  the 
additional  rings  had  been  mounted,  together  with  a  "sleeving"  of  the  aging  sec¬ 
tion  of  the  wind  tunnel  in  which  the  ions  were  generated.  Thus,  the  ions  did  not 
"see"  any  electric  field  prior  to  entering  the  free  volume  in  the  chamber  pro¬ 
per. 

In  figure  10  (section  14.  1.  ),  data  on  the  material  used  to  construct  the  chamber 
are  given  as  well. 

Figure  35  shows  the  assembled  chamber*.  To  reduce  capacity,  the  outer  elec¬ 
trodes  have  been  used  as  receiving  electrodes;  they  could  be  connected  to  the 
feedthroughs  of  the  wind  tunnel  with  the  shortest  length.  Figure  36  shows  the 
chamber  mounted  to  the  cover  plate  of  the  test  section  of  the  wind  tunnel,  and 

•The  individual  pans  of  Model  A,  prior  to  assembly,  have  been  shown  in  figures  2a,  2b,  and  2c  of  Quarterly  Progress  Report 
No.  1  0963)  of  the  same  contract. 
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METERS  OUTSIDE  OF  THE  WIND  TUNNEL 


in  figure  37  one  .sees,  above,  the  extremely  highly  insulated  feed-throughs  for 
the  receiving  electrode  and  the  box  for  inserting  the  electrometers,  and,  below, 
the  inputs  to  the  driving  electrodes  and  the  rings  and  to  the  shields  of  some  cri¬ 
tical  connections. 

Figure  38  shows  the  dummy  model  mounted  on  its  cover  plate  for  the  wind  tun¬ 
nel,  with  one  set  of  Pitot  tubes  inserted.  Another  set  of  Pitot  tubes  has  been 
shown  in  figure  30  (section  2.  1.  4). 
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Figure  37  COVER  PLATE  OF  TEST  SECTION  OF  WIND  TUNNEL  WITH  THE 
MODEL  A  OF  THE  MODIFIED  GtRDIEN  CHAMBER 
MOUNTED  BEHIND  IT 

(above:  Box  for  insertion  of  electrometers,  below:  5  x  19  connections  to 
electrodes  and  rings  and  shields  of  chamber  below  ) 

P12363C 
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Figure  38  DUMMY  MODEL  OF  MODEL  A  MODIFIED  GERDIEN  CHAMBER 
MOUNTED  ON  A  COVER  PLATE  FOR  TEST  SECTION  OF 
WIND  TUNNEL 

(Two  Pitot  tubes  are  inserted,  one  of  which  can  be  moved  radially  from 
outside . ) 

P12147A 
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2.  3  SUPPLY  AND  AMPLIFICATION 


2.  3.  1  PROVISION  OF  SUPPLY  VOLTAGES 

The  modified  GERDIEN  chamber,  Model  A,  has  16  receiving  electrodes  (which 
may  be  used  in  any  combination,  down  to  just  1)  plus  48  other  electrodes  (in¬ 
cluding  rings)  to  which  a  maximum  of  29  different  potentials  plus  the  "ground" 
potential  must  be  provided.  Except  for  the  ground  potential,  all  these  potentials 
are  functions  of  time.  Of  these  potentials,  27  are  sine  wave  potentials  oscillat¬ 
ing  versus  ground  as  their  zero  line,  one  is  a  sine  wave  plus  a  superimposed 
constant  potential,  and  3,nother  one  is  a  sine  wave  plus  superimposed  constant 
potential  plus  superimposed  square  wave  pulses,  triggered  by  the  sine  wave. 

In  addition,  a  sine  wave  with  a  180-degree  phase  shift  must  be  provided  for  the 
cancellation  of  the  main  displacement  current  (see  section  1.4.  5  and  figure  17). 

The  capacitive  load  for  the  supply  units  is  considerable,  since  the  generator  for 
the  ion  gate  pulses,  including  its  power  supplies,  and  the  "sleeving"  of  the  aging 
section  of  the  wind  tunnel  must  follow  the  sine  wave  variations  without  delay. 

For  this  reason,  a  specific  powerful  supply  for  ring  24  was  provided. 

For  the  experiments  in  the  laboratory,  it  was  desirable  to  be  able  to  vary  the 
supply  to  every  electrode,  if  possible,  independent  of  the  potential  applied  to 
other  electrodes.  Also,  independent  measurements  with  voltmeters  and  oscillo¬ 
scopes  of  each  of  these  potentials  was  desired.  The  final  schematic  of  the  over¬ 
all  supply  system  is  given  in  figure  39,  which,  if  taken  together  with  figure  10, 
hardly  needs  explanation.  Figure  39  also  contains  the  schematic  for  serving  the 
electric  valves  of  the  wind  tunnel  (see  section  2.  1.  2  and  figure  29).  Figure  40 
shows  the  special  rack  for  the  generation  of  the  supply  voltages.  The  related 
power  supplies,  however,  have  been  mounted  on  the  racks  shown  in  figure  32 
(section  2.  1.4.  )  and  located  in  some  distance  to  avoid  electromagnetic  disturb¬ 
ances  from  them.  A  careful  grounding  system,  using  own  grounders  free  of  any 
noise  in  a  loop  connection,  and  grounding  leads  poor  in  capacitance  and  inductance 
as  well  as  resistance,  is  essential. 

2.  3.  2  AMPLIFICATION  OF  THE  SIGNAL  FROM  THE  RECEIVING  ELECTRODES 

The  experiments  conducted  to  find  the  best  method  for  the  amplification  of  the 
signal  coming  from  the  receiving  electrodes,  and  thus  for  the  measurement  of 
numbers  of  electric  charges,  are  not  reported  here.  In  part,  they  have  been 
described  by  DOLEZALEK  and  OSTER  (1964). 

For  the  measurement  of  dc  voltages  or  currents  the  following  methods  have  been 
applied: 

The  LINDEMANN- Electrometer  by  Spindler  &  Hoyer 
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SUPPLY  CIRCUIT  FOR  MODIFIED  GERDIEN  CHAMBER 
FEBRUARY  1965 


EXPLANATION  OF  SYMBOLS 
o  BANANA  JACKS  ON  FRONT  PANEL 
Q  8NC -CONNECTORS  ON  FRONT  PANEL 
y  REGULATION  OONE  ON  SUPPLY  SWITCH  BOARD  PANEL 
/  REGULATION  OONE  ON  INSTRUMENT  PROPER 

- SHIELD 

X  LIGHTED  SWITCH 

N _ ••'CONNECTIONS  MADE  BY  CABLES  WITH  BANANA  PLUGS 
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"MODIFIEO  GERDIEN  CHAMBER" 


Figure  39  TOTAL  SUPPLY  SYSTEM  FOR  MODEL  A  OF  THE  MODIFIED  GERDIEN  CHAMBER 
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Figure  40  RACK  FOR  THE  PROVISION  OF  SUPPLY  VOLTAGES  TO  THE  RINGS 
AND  ELECTRODES  FOR  THE  MODEL  A  MODIFIED  GERDIEN  CHAMBER 
(The  lower  part ,  recessed  for  shielding  purposes/  is  the  pulse  generator  for 
the  ion  gate.) 
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The  Electrostatic  Voltmeter  by  Trvib,  Tauber  &  Co. 

The  Cary  31  Vibrating  Reed  Electrometer  by  Applied  Physics  Co. 

A  Valve  Electrometer  of  the  straight-forward  bridge  circuit  type 
(after  DOLEZALEK  1956,  1961). 

For  the  ac  measurements,  the  following  methods  were  used: 

The  Bridge-Electrometer,  after  DOLEZALEK,  with  a  sensitive 
double  beam  oscilloscope. 

The  Admittance  Neutralizer  by  Micronia  Corp.  with  oscilloscope 

A  combination  of  electrometer  input  stage,  operational  amplifier 
by  Philbrick  Co. ,  and  a  10^'  Ohm  feedback,  following  the  schematic 
by  SAPUPPO  (1964). 

The  schematic  of  the  DOLELALEK  bridge  electrometer  has  been  shown  by 
DOLEZALEK  and  OSTER  (1964).  There  also  the  basic  schematic  of  the  SAPUPPO 
electrometer  is  presented.  The  actual  schematic  applied  in  our  measurements 
is  given  in  section  3,  2.  3.  (figure  55,  sheet  2)  below,  while  calibration  results 
are  presented  in  figures  41,  a.  and  b.  The  results  represented  in  figure  41b 
have  later  on  been  improved  toward  better  sensitivity. 
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Figure  4la  CALIBRATION  CURVE  OF  THE  ELECTROMETER  AFTER  SAPPUPO  (1964),  WITH 
A  FEEDBACK  RESISTOR  OF  IQ12  OHM,  AND  AN  IMPUT  SERIES  RESISTOR  OF 

!010  OHM  -  DC  OPERATION 


65-  12065 


SIGNAL -FREQUENCY  j>.2j 


Figure  41b  MINIMUM  INPUT  o.c.  CURRENT,  WHICH  STILL  GIVES  A  READABLE  OUTPUT 
OF  70  mV  FROM  THE  ELECTROMETER,  AFTER  SAPPUPO,  WITH  A  W-  OHM 

FEEDBACK  RESISTOR 


The  curves  present  the  minimum  measurable  current  as  a  function  of  frequency  and 
of  input  capacitance  be'.ween  zero  (i.e  . ,  only  the  natural  capacitance  of  the  in¬ 
put)  and  100  x  10"^  Farad  (i.e.,  in  addition  to  the  natural  input  capacitance). 


2.  4  LABORATORY  EXPERIMENTS 


2.4.1  INTRODUCTION 

In  this  chapter,  laboratory  experiments  which  have  been  conducted  to  determine 
the  characteristic  features  of  the  modified  GERDIEN  chamber  and  to  measure 
mobility  spectrum  of  ions  generated  artificially  in  air  (from  pressurized  bof- 
tle^  are  reported  .  A  report  on  experiments  for  the  design  and  improvement  of 
the  wind  tunnel  and  the  facilities  to  measure  pressure  and  air  velocity  in  it,  as 
well  as  on  other  supplementory  experiments,  such  as  generation  of  the  supply 
voltages  and  pulse  voltages,  cancellation  of  displacement  current,  elimination 
of  noise,  and  amplification  of  the  signal  is  not  given  here.  These  experiments 
have  been  reported  in  the  monthly  and  quarterly  reports  prepared  for  this  con¬ 
tract. 

2.  4.  2  INITIAL  EXPERIMENTS 

After  an  extended  measuring  period  used  to  determine  all  insulation  resistances 
and  all  capacitances  between  the  rings  and  electrodes  in  the  modified  GERDIEN 
chamber,  at  first,  the  so-called  "GERDIEN characteristic"  was  measured. 

That  is  the  current  output  of  the  chamber  as  a  function  of  the  applied  driving 
voltage,  measured  in  constant  wind  speed  and  with  a  constant  ion  supply.  This 
was  done  under  atmospheric  pressure,  in  a  FARADAY  cage,  with  dc,  and  for 
four  combinations  of  receiving  electrodes.  Here,  the  inner  electrodes  have  been 
used  as  receiving  electrodes.  Figure  42  shows  the  result.  For  these  experi¬ 
ments,  the  edge  effect  at  the  intake  of  the  chamber  was  not  cancelled  (rings  00 
to  07  had  not  yet  been  applied,  and  no  sleeving  was  provided  for  the  ion  genera¬ 
tion);  for  this  reason,  the  right  hand  parts  of  the  characteristics  show  the  typical 
decrease  of  output  with  increasing  driving  voltage  and  thus  increasing  edge  ef¬ 
fect. 

• 

Next,  an  experiment  was  conducted  to  determine  the  maximum  voltage  ar  ,’icable 
to  the  ion  fence  (or  ion  filter,  rings  10  to  17)  without  causing  corona  disenarge. 
The  result,  presented  in  figure  43,  corresponds  to  the  well-known  PASCHEN 
curve,  and  may  be  explained  by  the  fact,  that  the  length  of  the  cathode  fall  in¬ 
creases  with  decreasing  gas  density,  and  that  a  self- sustained  gas  discharge 
cannot  occur  without  the  provision  Ox  this  length.  In  the  experiments,  the  onset 
of  corona  discharge  was  determined  by  measuring  electrometrically  the  current 
between  the  two  involved  rings. 

Then,  preliminary  ion  spectra  have  been  measured  in  the  low-pressure  wind 
tunnel,  under  relatively  high  pressure  and  applying  a  dc  driving  voltage  to  all 
driving  electrodes,  which  were  interconnected.  Here,  the  outer  electrodes  have 
been  used  as  driving  electrode.  The  result  is  shown  in  figure  44.  Air  pressure 
in  the  wind  tunnel  was  about  2  torr,  co"  responding  to  40  km. 
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Figure  42  DC  GERDIEN  CHARACTERISTICS 

(AlnMspUiic  pressu re,  artlf icioV  ionization  of  loborotor>  oir.  RiCtiving  tltctrod*  of 
modified  GERDIEN  chombrr  divided  into  four  port*  of  equol  length;  abscissa: 
driving  voltage;  ordinate:  current  from  Receiving  Electrode.  Tne  moximo  of  current 
ore  shifted  to  smaller  driving  voltages  with  increasing  distance  of  receiving  elec* 
Irode  from  front-(upstrean)  end  of  chamber.) 


Figure  43  ONE  EXAMPLE  FOR  THE  ONSET  OF  CORONA  DISCHARGE  BETWEEN  TWO 
ADJACENT  COAXIAL  RINGS  AT  THE  INTAKE  OF  THE  MODIFIED  GERDIEN 
CHAMBER 


ELECTRODE  NO 

65-12068 

P.g*e  44  REI  ATIVE  ION  SPECTRUM  OBTAINED  Wll  H  DIFFERENTIAL  ION  COUNTER 
OF  THE  SECOND  OROER  AND  DC  DRIVING  VOLTAGE 
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The  driving  voltage,  equal  for  all  driving  electrodes,  was  +10  V.  The  ion  fence 
was  closed.  The  ion  gate  was  kept  open  all  the  time;  it  was  located  between  the 
fourth  and  fifth  of  the  coaxial  rings,  counted  from  the  inner  electrode,  including 
the  ring  which  is  flush  with  it.  The  wind  velocity  was  very  high,  about  the  speed 
of  sound.  Because  of  the  edge  effect,  the  numbers  at  the  ordinate  of  figure  44 
should  be  taken  only  as  relative  ones. 

The  abscissa  shows  the  receiving  electrodes  which  have  been  measured  one  by 
one--a  total  of  16  electrodes.  The  first  electrode  counted  from  the  upstream 
(front-)  end  has  number  60,  the  last  one  at  the  rear  end  of  the  chamber  has  num¬ 
ber  75.  Thus,  with  increasing  numbers  of  electrodes,  the  mobility  of  intercepted 
ion  decreases.  It  is  obvious  that  ions  with  a  still  smaller  mobility  existed  as 
well  but  have  not  been  measured,  since  the  driving  voltage  was  not  high  enough 
for  this  wind  speed.  An  increase  of  driving  voltage  was  not  attempted,  since 
this  would  automatically  increase  the  edge  effect. 

Experiments  with  ac  driving  voltages,  applied  in  a  great  variety,  have  been 
conducted  as  well  in  the  preparatory  stage,  and  some  of  them  have  been  used 
to  derive  information  on  the  ion  spectrum  (see  DOLEZALEK  and  OSTER,  1964). 
On  this  basis,  the  conditions  for  the  actual  laboratory  measurements  have  been 
determined. 

2.4.  3  PROCEDURE  FOR  LABORATORY  MEASUREMENTS  IN  THE 
WIND  TUNNEL 

By  describing  this  procedure,  the  significance  of  the  many  different  parameters 
discussed  in  this  and  in  our  earlier  reports  will  become  clearer. 

At  first  an  altitude  to  be  simulated  is  selected,  and  the  appropriate  air  pressure 
is  derived  frdm  a  table  or  curve.  Then,  a  reasonable  wind  velocity  is  chosen, 
based  on  parachute  flights  with  instruments  made  by  previous  workers  and  re¬ 
calling  that  it  is  possible  to  reduce  the  wind  velocity  inside  the  chamber,  but 
hard  to  increase  it  if  the  fall  velocity  must  be  kept  below  the  speed  of  sound. 

For  the  selected  altitude,  a  range  of  expected  ion  mobilities  is  estimated,  ex¬ 
tending  between  kj  and  k2*  From  these  values,  the  U/cu  values  are  determined. 
For  the  case  of  a  plate  condenser  with  distance  ^between  the  plates,  and  dis¬ 
tance  y^  between  the  ion  gate  and  receiving  electrode,  the  equation  2  U^/co  = 
dk/yj  is  easily  derived,  where  U^  is  the  driving  voltage  measured  from  zero  to 
peak.  For  the  present  model  of  modified  GERDIEN  chamber,  the  U/a>  value  can 
be  taken  directly  from  figure  15  in  section  1.4.4.  Thus,  two  U/<y  values  are  de¬ 
termined,  corresponding  to  maximum  and  a  minimum  expected  mobility  (kj  and 

k2). 

The  driving  voltage  frequency  will  be  the  same  for  all  electrodes,  while  the  am¬ 
plitude  U  will  increase  along  the  axis  of  the  chamber.  For  the  determination 
of  the  frequency,  we  use  the  equation  <u>2ffv/Ax,  where  v  is  the  wind  speed,  and 
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Ax  the  length  of  one  partial  receiving  electrode.  This  equation  is  obvious  from 
the  fact  that  at  least  one  full  cycle  of  driving  voltage  should  occur  while  an  ion 
is  passing  in  front  of  one  partial  receiving  electrode.  Only  then  are  we  sure 
that  it  will  reach  this  electrode,  and  the  ion  with  the  next  smaller  mobility  will 
reach  the  next  electrode.  Thu3,  we  get  a  minimum  frequency,  and  now  with  our 
two u/<u  values,  two  minimum  voltages  for  the  first  and  the  last  driving  electrodes. 

In  fact,  we  shall  apply  this  higher  voltage  to  the  driving  electrode  which  is  next 
to  last,  and  apply  a  much  higher  voltage  to  the  very  last  electrode  order  in  to 
ensure  that  all  ions  are  intercepted.  Then  the  intermediate  voltages  are  selected 
to  follow  either  an  arithmetical  or  geometrical  series  or  some  other  sequence. 

The  next  steps  aim  at  the  determination  of  opening  time  r  of  the  ion  gate.  Un¬ 
fortunately  we  are  not  completely  free  to  select  this  time,  if  the  effective  length 
of  the  ion  gate  l  and  the  length  of  the  partial  receiving  electrodes  Ax  are  of  the 
same  order.  To  understand  this,  see  equation  <o>2nv/Ax  and  equation  III  from 
figure  9  in  which  we  set  q  =  r/T,  and  take  the  case  that  the  second  right-hand 
side  member  is  small  (E  sufficiently  large):  2^vq 

®  <  l 

By  inserting  our  minimum  frequency  o>,  derived  from  wind  speed  and  length  of 
receiving  electrodes,  into  equation  III  we  determine  l,  the  fraction  of  opening 
time  of  ion  gate  to  the  full  cycle  of  driving  electrode,  and  can  set  our  ion  gate 
pulse  generator  accordingly  (consider,  however,  the  significance  of  length 
as  discussed  above). 

We  are  now  ready  to  perform  the  measurement.  At  first  we  determine  whether 
sufficient  ions  are  created  by  using  the  modified  GERDIEN  chamber  as  an  ion 
counter  with  a  large  dc  driving  voltage,  then  we  determine  the  actual  length  of 
ion  gate  as  indicated  above,  and  apply  the  calculated  driving  voltages  and  fre¬ 
quencies,  and  measure  the  current  coming  from  each  of  the  receiving  electrodes, 

2.  4.  4  DISCUSSION  OF  A  RESULT  FROM  THE  LABORATORY  MEASUREMENTS 

• 

As  discussed  in  section  1.  4.  2, ,  the  modified  GERDIEN  chamber  can  be  operated 
in  two  ac  modes:  as  a  differential  chamber  of  the  first  order,  or  as  a  differen¬ 
tial  chamber  of  the  second  order.  The  result  of  an  experiment,  made  under 
otherwise  equal  conditions  with  these  two  modes,  is  shwon  in  figures  45  and  46. 

For  the  measurement  according  to  the  differential  chamber  on  the  first  order, 
all  receiving  electrodes  have  been  connected  as  well  as  all  driving  electrodes. 

The  driving  voltage  was  swept,  while  the  u/cj  value  was  kept  constant,  and  the 
current  from  the  receiving  electrode  was  measured.  The  results  have  been 
plotted,  and  the  mobility  spectrum  has  been  derived  by  graphical  differentiation. 
Unless  the  curve  is  a  polygonic  trace,  it  may  be  difficult  to  decide  how  many 
tangents  can  be  put  to  it  and  where  to  place  them.  The  more  tangents  one  may 
reasonably  apply,  the  better  becomes  the  resolution  but  the  accuracy  suffers. 

In  particular,  it  may  be  difficult  to  determine  whether  any  smooth  transition 
between  straighter  parts  of  the  curve  presents  the  existence  of  ions  of  intermediate 
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Figure  46  MOBILITY  SPECTRUM  (Differential  chamber  of  the  second  order; 
pressure  =  0.8  Torr.) 


mobility  or  is  just  the  result  of  the  interception  of  ions  straying  because  of  dif¬ 
fusion.  We  decided  to  divide  the  curve  into  eight  different  slopes.  The  result 
is  shown  in  figure  45. 

For  the  measurement  with  the  differential  chamber  of  the  second  order,  a  smooth 
increase  of  driving  voltage  was  applied  by  bringing  different  voltages  in  small 
steps  to  each  of  the  16  driving  electrodes.  The  receiving  electrodes  have  been 
divided  into  only  three  portions-- 5,  5,  and  6  individual  electrodes.  Thus,  the 
resolution  is  not  better  than  given  by  this  number  of  three,  and  it  cannot  be  said 
where  the  upper  limit  of  mobility  was  situated.  On  the  other  side,  the  accuracy 
is  much  better.  The  reliability  of  the  measurement  could  be  demonstrated  by 
applying,  e.g.,  the  driving  voltage  responsible  for  one  receiving  electrode  to 
the  next  one:  no  more  ions  were  then  intercepted  there. 

The  number  of  ions  per  mobility  range,  of  course,  is  determined  to  the  degree 
of  artificial  ionization  maintained  in  the  wind  tunnel  during  the  experiments.  Thus, 
the  height  of  the  columns  in  figures  45  and  46  does  not  give  an  information  on 
natural  conditions.  However,  if  the  two  figures  are  compared  to  each  other,  it 
comes  out  that  for  the  same  mobility  range  the  same  number  of  ions  was  measured 
with  the  two  methods,  thus  again  confirming  the  reliability  of  the  measurement. 

If  we  compare  the  measured  mobilities  with  the  values  expected  for  the  height 
of  about  50  km  in  the  free  atmosphere  (see  figure  47,  which  is  plotted  from  table 
202  in  section  1.  6.  4,  and  table  103  in  section  1.  6.  3.  ),  we  find  that  the  coinci¬ 
dence  is  good.  No  attempt,  however,  was  made  to  determine  the  chemical  na¬ 
ture  of  the  ions  measured.  The  pollution  of  the  wind  tunnel  seemed  to  be  too 
great  to  allow  a  reliable  determination,  since  no  special  precaution  was  taken 
to  secure  a  completely  dust-free  air  for  these  experiments. 
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Figure  ESTIMATED  MOBILITY  VALUES  FOR  DIFFERENT  HEIGHTS 


PART  3 

DESIGN  OF  ION-SPECTROMETER 
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3.  1  DETERMINATION  OF  DIMENSIONS  FOR  MODEL  B1 


3.  1.  1  INTRODUCTION 

Based  on  the  theoretical  considerations  reported  in  Part  1  of  this  report,  and 
on  the  laboratory  experiments  described  in  Part  Z,  the  modified  GERDIEN  cham¬ 
ber  or  ion  spectrometer  for  the  measurements  in  the  free  atmosphere  has  been 
designed.  While  the  actual  measurements  in  the  mesophere  and  upper  stratos¬ 
phere  will  be  conducted  by  the  instrument,  falling  down  tinder  a  parachute  from 
a  height  of  about  80  km  (to  which  it  must  be  carried  by  a  rocket),  it  is  intended 
to  execute  at  first  a  series  of  measurements  in  lower  heights.  It  is  supposed 
that  the  instrument  will  be  carried  to  about  30  km  by  a  balloon  and  then  make  the 
measurements  while  falling  down  under  a  parachute.  Recovery  is  intended. 
Measurements  in  the  lower  stratosphere  will  be  made  by  applying  a  dc  driving 
voltage  and  using  the  instrument  as  a  differential  chamber  of  the  second  order. 

With  these  intentions,  it  was  decided  to  build  the  instrument  proper  (i.  e.  ,  the 
receiving  and  driving  electrodes,  the  aft  rings,  and  the  hull)  to  be  suitable  for 
the  ac  measurements  at  greater  heights,  but  to  attach  the  front  rings  which  are 
suitable  for  dc  measurements  at  lower  heights.  These  rings  can  be  replaced 
by  the  shorter  rings  for  ac  measurements,  if  the  instrument  is  recovered  in  a 
good  state  after  the  dc  measurements.  The  supply  unit  was  constructed  for  the 
dc  measurements,  while  the  amplifiers  will  be  the  same  ones  for  both  series. 

The  general  decision  of  the  size  of  the  chamber  v/as  derived  from  the  considera¬ 
tions  communicated  in  section  1.  4.  8  of  this  report.  There  it  has  been  stated 
that  the  length  of  the  individual  receiving  electrode  must  be  greater  than  the  length 
of  the  rings  constituing  the  ion  gate.  Since  we  did  not  want  to  reduce  the  length 
of  the  ion  gate  very  much,  we  had  to  make  the  receiving  electrodes  much  longer 
than  in  Model  A.  This  has  the  consequence  that  the  overall  length  of  the  chamber 
also  increased  considerably;  it  is  no^  about  90  cm. 

The  principal  technological  difficulty  was  given  by  the  fact  that  the  instrument 
should Withstand  accelerations  of  300  g,  which  may  occur  in  the  process  of  e- 
jecting  the  instrument  with  the  parachute  from  the  rocket.  Another  difficulty 
for  measurements  in  the  free  atmosphere  is  given  by  the  necessity  of  avoiding 
losses  of  insulation  caused  by  the  possible  condensation  of  water  vapor  on  the 
insulating  surfaces,  and  of  maintaining  a  constant  temperature  in  the  enviroment 
of  the  high-ohm  resistors  in  the  amplifiers  which  have  a  great  temperature  co¬ 
efficient.  To  overcome  these  difficulties,  a  heating  is  provided:  (a)  for  the  main¬ 
tenance  of  a  temperature  at  the  insulation  surfaces  slightly  higher  than  that  of 
the  enviroment  by  covering  the  whole  chamber  with  "heating  blankets"  (bifilar 
wiring  to  avoid  magnetic  fields),  and  (b)  in  the  preamplifiers  by  a  heater  and  a 
thermostat  and  a  temperature  probe.  In  addition,  ample  temperature  insulation 
has  been  applied. 
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3.  1.  2  DETERMINATION  OF  DESIGN  PARAMETERS 


The  inner  electrode  will  serve  as  the  driving  electrode.  It  is  divided  into  45  indi¬ 
vidual  rings,  numbered  201  to  245,  and  insulated  from  each  other.  All  rings 
are  connected  to  wires,  which  are  led  out  through  the  boring  in  the  inner  elec¬ 
trode  and  which  lead  to  two  connectors  mounted  on  the  outside  of  the  chamber 
{but  under  the  outer  layer  of  thermal  insulation).  By  plugging  other  connectors 
into  them,  which  contain  resistors,  any  distribution  of  driving  voltages  to  the 
45  individual  electrodes  can  be  supplied  (see  figure  55,  sheet  6).  For  the  dc 
experiment,  all  individual  driving  electrodes  are  interconnected. 

The  outer  electrodes  will  be  used  as  receiving  electrodes.  There  are  fifteen  of 
them,  but  for  the  dc  experiments  they  will  be  connected  in  groups,  resulting  in 
a  number  of  only  six  actual  electrodes.  These  receiving  electrodes  are  mounted 
on  very  small  teflon  insulating  stripes,  so  that  most  of  the  dielectricum  between 
them  and  the  outer  hull  is  air.  This  was  done  to  reduce  the  capacitance.  To 
avoid  conduction  through  the  air  between  them  and  the  outer  hull,  the  backside 
of  the  receiving  electrodes  is  covered  with  a  thin  layer  of  highly  insulating  ma¬ 
terial.  The  receiving  electrodes  are  made  of  stainless  steel.  The  preamplifiers 
will  be  mounted  directly  over  the  heating  blankets  which  serve  as  the  hull  around 
the  receiving  electrodes.  These  heating  blankets  are  covered  with  a  highly  con¬ 
ductive  layer  for  electrostatic  shielding. 

Since,  for  a  dc  operation,  no  ion  gate  pulsation  is  needed,  a  number  of  adjacent 
rings  of  the  three  different  sets  at  the  upstream  end  can  be  electrically  connf  cted. 
In  this  case,  they  have  been  made  simply  of  one  piece.  The  rings  of  the  ion 
fence  or  ion  filter,  numbered  1 1  to  16,  are  much  longer  than  in  Model  A  to  avoid 
the  otherwise  necessary  great  fence  voltages  for  the  measurements  in  the  lower 
stratosphere,  which  could  give  rise  to  corona  discharges  (see  sections  1.4.7 
and  2.4.  2).  The  first  set  of  rings,  serving  to  reduce  the  edge  effect,  is  numbered 
01  to  06,  while  the  third  set  at  the  upstream  end  and  the  set  at  the  downstream 
end  are  numbered  21  to  26  and  31  to  3 £>,  respectively  (see  figure  53),  all  numbers 
counting  from  inward  out. 

3.  1.  3  DETERMINATION  OF  EDECTRICAL  PARAMETERS 

To  derive  the  electrical  parameters  for  Model  Bl,  we  started  from  figure  48, 
which  presents  the  results  of  calculations,  made  on  the  basis  of  the  geometri¬ 
cal  dimensions  of  the  chamber  as  well  as  on  some  assumptions  on  the  air  velo¬ 
city  inside  the  chamber.  These  calculations  concerned  the  location  of  the  re¬ 
ceiving  electrodes  for  different  driving  dc  voltages,  under  the  assumption  of  cer¬ 
tain  mobilities  being  theoretically  representative  for  the  different  heights  between 
40  and  10  km.  Since  these  mobilities  have  been  theoretically  extrapolated  from 
measured  values  at  much  lower  altitudes,  they  are  rather  unexact.  For  this 
reason,  curves  have  been  added,  which  correspond  to  mobilities  being  by  factors 
of  3  either  greater  or  smaller  than  the  assumed  values. 
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To  arrive  at  the  results  of  figure  48,  at  first  certain  calculations  had  to  be  made 
with  respect  to  the  air  velocity  inside  the  chamber.  Based  on  the  actual  dimen¬ 
sions,  and  on  a  number  of  "pessimistic"  presumptions,  REYNOLDS'  numbers 
have  been  calculated,  and  from  these  the  maximum  velocities  may  be  determined 
by  experiments  in  the  wind  tunnel.  Table  1  gives  these  values. 

TABLE  1 

MAXIMUM  AIR  VELOCITIES  IN  A  MODIFIED  GERDIEN 
CHAMBER  FOR  LAMINARY  FLOW 
(Pessimistic  Values  for  the  Later  Theoretical  considerations) 


Height  (km) 

Maximum  Velocity  Inside  Chamber  (m/s) 

10 

5.  15 

20 

11.  08 

30 

24.  65 

40 

52.  07 

For  the  derivation  of  these  numbers,  a  flat  wind  profile  inside  the  chamber  was 
assumed.  A  correction  must  be  made  after  the  flight,  when  the  actual  average 
velocities  for  any  given  pressure  are  known  and  the  actual  wind  profile  can  be 
better  estimated. 

The  curves  of  figure  48  suggest  a  subdivision  of  the  15  receiving  electrodes  of 
the  chamber  to  the  six  planned  electrometer  amplifiers,  according  a  geometri¬ 
cal  series.  It  was  decided  to  connect  the  amplifiers  as  shown  in  table  2, 

TABLE  2 

CONNECTION  OF  ELECTROMETER  AMPLIFIERS  ■ 

TO  RECEIVING  ELECTRODES 


Amplifier  No. 

Is  Connected  to  Receiving  Electrode  No. 

A1 

101 

A2 

102 

A3 

103,  104 

A4 

105,  106,  107 

A  5 

108,  109,  HO,  111 

A6 

- - 

112,  113,  114,  115 
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It  had  been  decided  at  an  earlier  time  that  the  driving  voltages  for  the  dc  ex¬ 
periment  with  the  GERDIEN  chamber  would  be  stepped  through  the  whole  intended 
range  over  the  whole  altitude  range  of  the  flight.  According  to  the  curves  of  fig¬ 
ure  48,  10  V  would  be  a  good  value  for  the  driving  voltage  in  an  altitude  of  40  km, 
and  300  V  would  be  good  for  10  km.  It  was  decided,  however,  to  extend  the  range 
of  driving  voltages  from  5  to  300  V,  because  this  has  the  effect  of  increasing  the 
dynamic  range  for  the  measurement  of  current  coming  from  the  receiving  elec¬ 
trodes,  at  the  lower  end,  by  a  factor  of  2.  Under  these  circumstances,  and  a 
subdivision  of  driving  voltages  according  to  a  geometrical  series,  the  values  of 
table  3  have  been  chosen. 


TABLE  3 

DRIVING  VOLTAGES  FOR  MODIFIED  GERDIEN  CHAMBER 
(d.  c.  Experiments,  10,  25,  60,  140,  300  V) 

Every  voltage  is  applied  simultaneously  to  all  driving  electrodes,  Nos.  201  to 
245,  which  are  all  interconnected  for  this  particular  experiment. 

The  next  decision  had  to  be  made  with  respect  to  the  "programming,  "  the  se¬ 
quence  and  duration  of  switching  processes.  Two  requirements  had  to  be  con¬ 
sidered.  At  first,  the  commutator  of  the  telemetry  channel  allows  one  measure¬ 
ment  per  channel  each  0.  2  second.  Secondly,  if  one  full  cycle  is  accomplished 
in  the  time  of  1  minute,  we  shall  obtain  one  full  cycle  of  measurements  for  the 
following  falling  distances: 


TABLE  4 

ROUGH  VALUES  FOR  MEASURING  CYCLES  OF  D.  C,  EXPERIMENT 
(Starting  from  an  Assumed  Height  of  40  Km) 


« 

40  .  .  . 

37  km 

37  .  .  . 

35  km 

35  ... 

33  km 

33  ... 

31.5  km 

31.  5  . 

.  .  30.  0  km 

30.  0  . 

28.  5  km 

28.  5  . 

.  .  27.  0  km 

about  all  600  to  1200  m  from  27  to  20  km 

about  all  30  to  600  m  from  20  km  down  to  10  km  altitude. 
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These  numbers  are  very  rough,  because  they  do  not  take  into  account  the  differ¬ 
ence  between  air  speed  in  the  chamber  and  falling  speed  of  the  instrument  with 
respect  to  a  ground  fixed  frame.  However,  a  number  of  approximately  80  to 
150  full  cycles  could  be  measured  during  the  falling  from  40  to  10  km.  Thus, 
a  duration  of  a  full  cycle  of  1  minute  is  certainly  not  too  long. 

Each  cycle  consists  of  six  measured  curves.  This  may  be  explained  in  the  fol¬ 
lowing  way.  Before  a  particular  driving  voltage  is  applied,  all  receiving  elec¬ 
trodes  are  shortly  grounded.  After  they  are  disconnected  from  ground  and  the 
driving  voltage  is  applied,  ions  will  begin  to  sediment  on  the  receiving  electrodes 
and  will  raise  their  potential  in  the  form  of  a  curve  of  the  type  y  =  1  - 
exp  (-x).  Since  for  each  cycle  six  different  driving  electrodes  are  applied,  this 
will  happen  six  times  for  each  cycle,  i.  ev  for  every  curve  there  is  a  time  of 
roughly  10  seconds.  In  this  period,  "10"/0.  2 ^45  single  points  of  these  curves 
will  be  measured.  Unless  a  curve  has  reached  its  maximum  value  in  a  period 
less  than  0.4  second,  the  curve  itself  can  be  determined  from  these  points.  This 
remains  true,  even  if  the  time  to  reach  maximum  is  considerably  greater  than 
9  or  10  seconds.  Thus,  the  maximum  value  may  be  determined  by  calculation 
even  in  these  cases,  though  it  is  preferred  to  directly  measure  the  height  of  the 
curve.  There  is  a  great  chance  that  during  each  cycle  there  will  be  at  least  one 
case  in  which  this  is  possible  for  each  of  the  six  amplifiers.  This  is  a  conse¬ 
quence  of  the  stepping  of  driving  voltages  according  to  table  3. 

It  is  furthermore  planned  that  each  amplifier  be  checked  on  time  during  each 
cycle,  and  that  the  zero  point  be  tested  six  times  during  each  cycle.  This  "check¬ 
ing"  of  the  amplifier  is  not  a  full  calibration,  because  we  did  not  want  to  include 
the  high-ohmic  resistance  in  this  control  process  so  as  not  to  endanger  the  good 
insulation.  Also,  a  correct  calibration  of  an  electrometer  takes  a  considerable 
amount  of  time.  Always,  when  the  driving  voltage  is  changed,  the  amplifiers 
are  grounded.  Thus,  a  saturation  of  the  amplifiers  due  to  displacement  current 
is  avoided. 

• 

For  the  first  flight,  only  ions  of  one  sign  will  be  measured.  The  instrument  is 
built  in  such  a  way,  however,  that  by  desoldering  one  wire  it  will  be  able  to 
measure  alternatively  ions  of  both  signs,  Then,-  one  cycle  would  measure  pos¬ 
itive  ions,  the  next  negative  ions,  etc.  It  should  be  noted  that  the  values  of 
table  4  are  derived  from  maximum  falling  speeds.  If  measurements  of  less  dis¬ 
tant  regions  are  desired,  and  if  the  falling  speed  of  the  instrument  can  be  reduced 
by  a  greater  parachute,  the  spacing  could  be  made  much  closer. 

Finally,  it  had  to  be  determined  how  sensitive  the  electrometer-amplifiers  had 
to  be  made.  The  general  rule  applied  here  was  that  for  the  first  flight  it  should 
be  preferred  to  risk  saturation  and  the  risk  of  missing  measurements  due  to  in¬ 
sufficient  sensitivity  should  be  reduced. 

The  calculation  has  been  made  for  ion  number  densities  of  10®  and  10*®  m"^, 
since  the  expected  values  should  be  in  between  these  limits.  Considering  the 
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actual  area  of  the  ion  gate,  the  numbers  of  incoming  ions  has  been  calculated. 
From  these,  and  by  considering  figure  3,  the  following  limitations  for  the  cur¬ 
rent  measurements  have  been  derived  (table  5). 

TABLE  5 

LIMITS  OF  CURRENT  MEASUREMENT  BY  THE 
ELECTROMETER  AMPLIFIERS 
(Assuming  a  Mobility  Range  of  IQ 


g 

Currents  per  Amplifier  for  the  Driving  Voltages  of: 

300  Volt 

10  Volt 

5  Volt 

10 

average  per 
amplifier: 

0.  1  .  .  .  1  pA 

less  than 

0.  1  pA 

less  than 

0.  1  pA 

40 

amplifier  A1 
gets 

5  ...  50pA; 

A2  .  .  .  A6:  nothing 

average  per 
amplifier: 

0.  8  .  .  .  8  pA 

amplifiers 

Ai  . .  .  A3: 

nothing' 

amplifiers 

A4  . .  .  A6 
in  the 
average 

0.  8  . .  .  8  pA 

The  dynamic  range  of  the  electrometer-amplifier  can  be  made  rather  great, 
since  the  amplifier  is  able  to  measure  less  than  1  fA.  However,  the  dynamic 
range  of  the  telemetry  system  is  limited  to  about  25.  Remembering  our  general 
decision,  rather  being  too  sensitive  than  too  insensitive,  we  propose  a  dynamic 
range  of  the  electrometer  amplifiers  as  follows; 

TABLE  6 

DYNAMIC  RANGE  OF  ELECTROMETER  AMPLIFIERS 

0.  1  .  .  .  25  pA  . 


With  this  x  -\ge,  we  should  measure  at  least  something  for  all  cases  of  table  5, 
and  in  the  coses  where  the  expected  current  is  greater  than  25pA,  we  should  be 
able  to  extrapolate  the  measured  curve  in  order  to  derive  the  value  of  equilibrium 
between  charging  of  the  receiving  electrode  by  sedimented  ions  and  discharging 
it  by  amplifier  input  resistance.  We  should  remember,  however,  that,  due  to 
the  stepping  of  driving  voltages,  every  point  of  the  mobility  spectrum  is  measured 
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or  approximated  several  times.  All  these  values  will  be  checked  as  soon  as  the 
calculation  is  made  for  the  bent  wind  velocity  distribution  curves.  However,  we 
do  not  expect  essential  variations. 

It  was  decided  to  monitor  the  driving  voltages  all  the  time.  Since  the  dynamic 
range  of  driving  voltages  is  greater  than  the  dynamic  range  of  the  transmitter, 
two  transmitter  channels  have  to  be  used,  which  previously  were  called  Alpha 
and  Beta  here.  Table  7  shows  the  applied  voltages. 

TABLE  7 

MONITORING  OF  DRIVING  VOLTAGES 


Output  from 
Switches  81 
to  86  (Volt) 

Transmitting 

Channel 

Alpha  Gets 

Transmitting 

Channel  Beta 

Gets 

300 

2.  5  V 

2.  5  V  (zener) 

140 

1.  16  V 

2.  5  V  (zener) 

60 

500  mV 

2.  5  V  (zener) 

25 

(208  mV) 

2.  5  V 

10 

(  83  mV) 

1.  0  V 

5 

(  42  mV) 

500  mV 
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3.  2  CONSTRUCTION  AND  INTEGRATION  OF  CHAMBER 


3.  2.  1  OVERALL  DESIGN 

Mechanically,  the  whole  instrument  will  consist  of  four  parts.  When  falling  down, 
the  foremost  part  is  the  chamber  proper.  It  consists  of  the  inner  and  outer  elec¬ 
trodes,  the  insulators  for  the  outer  electrodes,  the  heating  blankets  as  a  hull 
around  the  outer  electrodes  with  some  distance  in  between,  the  preamplifier  di¬ 
rectly  attached  to  pins  reaching  from  the  outer  electrodes  through  the  blankets 
without  touching  them,  the  connector  plugs  for  the  subdivision  of  driving  volt¬ 
ages  and  the  connectors  for  the  supply  to  the  rings,  the  front  and  aft  sets  of 
rings,  the  thermal  insulation  layers,  the  thin  outer  hull,  and  an  intake  ring  in 
the  form  of  half  a  doughnut. 

This  chamber  is  made  of  a  frame  consisting  of  four  beams  in  axial  direction 
and  bracers  forming  a  periphery.  The  four  beams  are  extended  backwards  (up¬ 
wards  when  falling)  to  allow  the  uninhibited  exit  of  the  air  streaming  through  the 
chamber,  and  ending  at  a  conical  cover  for  the  other  parts. 

Behind  (above  when  falling)  this  cone  there  is  at  first  a  setup  of  three  "decks.  " 
Two  of  them  carry  the  second- stage  amplifiers,  while  the  third  one  is  carrying 
the  supply  units  and  the  programmer.  After  these  three  decks  there  come  the 
compartments  for  the  28-V  battery  (batteries)  and  for  the  transmitter.  At  last 
comes  the  parachute.  The  parachute,  the  transmitter,  and  the  batteries  will  be 
supplied  by  the  contracting  agency  and  are  not  described  here.  The  batteries 
deliver  28  Vdc,  unregulated. 

3.  2.  2  THE  CHAMBER  PROPER  AND  THE  "DECKS" 

Figures  49,  50,  51,  and  52  show  the  chamber  proper.  Figure  49  shows  the  cham¬ 
ber  without  the  preamplifiers  and  the  thermal  insulation;  one  looks  at  the  frame 
and  at  the  heating  blankets  between  the  parts  of  the  frame.  In  figure  50,  the  ther¬ 
mal  insulation  has  been  added  partially,  while  figure  51  shows  the  fully  assembled 
instrument  with  the  outer  hull.  In  figure  52,  the  same  is  shown  as  seen  from  the 
front  end.  The  ruler  is  scaled  in  inches. 

Figure  53  gives  an  overall  sketch  of  the  instrument,  partially  cut  to  show  in¬ 
ternal  parts.*  Photographs  of  the  internal  parts  have  been  submitted  with  Monthly 
Report  No.  16  of  this  contract. 

In  figure  54,  the  assembled  decks  can  be  seen.  The  lowest  deck  and  the  middle 
deck  carry  the  six  second- stage  amplifiers,  essentially  consisting  of  operational 
amplifiers  (type  Philbrick  PP  65a,  solid  state).  The  uppermost  deck  caries  the 
cylindrical  programmer,  and  below  it  at  the  left  hand  side,  the  four  converters 


•A  complete  set  of  detail  drawings  in  scale  1:1  has  been  made,  but  is  not  included  in  this  report. 
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for  the  supply  of  regulated  voltages,  and  at  the  right  hand  side  the  box  for  the 
polarity  inverter. 

The  programmer  was  developed  and  constructed  by  the  Geophysics  Section  of 
Avco  RAD.  Two  arms  which  carry  small  magnets,  driven  by  an  electric  motor 
by  the  way  of  a  great  set,  rotate  in  a  completely  sealed  cylindrical  box.  The  magnetic 
fields  penetrate  the  hull  of  this  cylinder  and  may  so  activate  magnetic  reed 
switches,  which  are  mounted  on  the  outer  side  of  the  cylinder,  inserted  in  rings 
of  plexiglass  (covered  by  the  black  tape).  These  reed  switches  operate  in  a  sealed- 
in  atmosphere  themselves.  Thus,  motor,  rotating  arms,  and  switching  action 
are  independent  of  the  outer  atmosphere  and  its  pressure. 

The  polarity  inverters  consist  essentially  of  a  relay,  triggered  by  the  program¬ 
mer  once  with  every  revolution.  The  time  for  one  revolution  of  the  programmer 
is  set  to  approximately  1  minute. 

3.  2.  3  THE  ELECTRICAL  SCHEMATIC 

Figure  55,  consisting  of  eight  sheets,  presents  the  whole  electrical  schematic 
of  the  system,  with  the  exception  of  the  transmitter. 

Sheet  1  shows  the  general  connections  between  the  electrodes  and  rings  and  pre¬ 
amplifiers  mounted  on  the  chamber  proper  on  one  side,  and  the  second- stage 
amplifiers  as  well  as  the  supply  units,  programmer,  and  connections  to  the  trans¬ 
mitter  on  the  other  side. 

The  other  sheets  give  the  details  of  individual  parts,  including  a  survey  on  the 
electric  program.  They  are  understandable  without  further  explanation. 

The  preamplifiers  are  small  units  of  about  6x4x2  cm  in  a  curved  shielding  box, 
mounted  on  a  printed  circuit  board.  Including  the  small  heater,  the  thermostat, 
and  the  heat  probe,  they  are  potted  with  a  heat-conducting  material,  but  the  in¬ 
put  connections,  which  have  to  maintain  an  extremely  high  insulation',  are  not 
included  in  the  potting.  The  second- stage  amplifiers  are  mounted,  three  by 
three,  on  printed  boards  as  well.  The  voltage  inverters  are  industry  items. 
Photographs  of  the  printed  boards  as  well  as  of  the  details  of  the  programmer 
have  been  submitted  with  Monthly  Report  No.  16.  Vacuum  tightness  is  obtained 
with  O- rings. 

Part  of  the  electrical  connections  are  fully  shielded  electrically.  This  has  been 
made  possible  by  the  application  of  special  connectors  which  allow  the  insertion 
of  unshielded  pins  as  well  as  shielded  ones.  Sheet  1  of  figure  55  indicates  the 
shielded  connections  as  well  as  the  placing  of  connectors. 

The  requirement  for  power  for  the  heating  blankets  is  considerable.  However, 
it  has  been  calculated  that  the  applied  28-V  batteries  will  be  sufficient.  If  nec¬ 
essary,  the  supply  for  the  heating  blankets  can  be  made  independent  of  the  sup¬ 
ply  for  the  other  parts. 
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Figure  49  THE  MODIFIED  GERDIEN  CHAMBER,  MODEL  Bl,  WITHOUT 
PREAMPLIFIERS  AND  THERMAL  INSULATION 
P13471A 
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Figure  51  MODIFIED  GERDIEN  CilAMBER,  MODEL  Bl,  FULLY  COVERED 

P13471C 
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Figure  55  ELECTRICAL  CONNECTIONS  FOR  MODIFIED  GERDIEN  CHAMBER 

MODEL  B-l 

Sheet  2:  Left:  Schematic  of  Preamolifier  (Electrometer) 

Right:  Schematic  of  Second-Stage  Amplifier  (Operational  Amplifier) 
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Figure  55  ELECTRICAL  CONNECTIONS  FOR  MODIFIED  GERDIEN  CHAMBER  MODEL  B- 
Sheet  4:  Time  Schedule  of  Programmer 
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Figure  55  ELECTRICAL  CONNECTIONS  FOR  MODIFIED  GERDIEN  CHAMBER  MODEL  B-l 
Sheet  5:  Additional  Resistor  Network  on  "Deck  3"  of  Payload 
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Sheet  8:  Wiring  of  Connector  Plug  for  Aft  Rings 


3,  3  TESTING  OF  MODEL  B1 


After  the  chamber  had  been  built  and  the  electronics  have  been  completed,  the 
wind  tunnel  was  enlarged  to  suit  the  new  instrument.  At  first,  the  new  chamber 
was  checked  under  atmospheric  pressure;  it  turned  out  to  be  extremely  sensitive 
to  the  slightest  amounts  of  ionization  or  pollution  in  the  laboratory  air. 

Then,  the  chamber  was  put  into  the  wind  tunnel.  A  dummy  model  for  the  pro¬ 
grammer  was  built  and  connected.  While  the  preamplifier  operated  under  the 
low  pressure  in  the  wind  tunnel,  the  second-stage  amplifiers,  the  supply, and  the 
programmer  were  installed  outside  the  tunnel,  the  connections  made  by  the  way 
of  feedthroughs.  A  pressure  of  a  few  torr  was  applied,  with  a  wind  speed  cor¬ 
responding  to  the  data  presented  in  section  3.  1.  3.  Instead  of  a  transmitter,  a 
sensitive  Sanborn  eight-channel  recorder  was  connected.  By  varying  the  driving 
voltage,  it  could  be  shown  immediately  that  the  new  chamber  was  sensitive  and 
delivered  the  data  for  the  mobility  spectra  in  the  required  amount. 

During  the  time  for  the  preparation  for  the  planned  experiments  in  the  free  atmos¬ 
phere,  the  chamber  will  be  used  in  the  laboratory  for  a  variety  of  experiments. 
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